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Preface 


This  is  the  fust  volume  of  a  two-volume  collection  of  papers  presented  at  the 
Symposium  on  the  Modification  of  Electromagnetic  Scattering  Cross  Se  ctions  in  the 
Resonant  Region  held  on  4  June  1963  at  zhe  Air  Force  Cambridge  Research 
Laboratories. 

Papers  in  the  present  volume  establish  the  need  for  techniques  which  signifi¬ 
cantly  reduce  the  electromagnetic  scattering  cross  section  of  bodies  in  the  fre¬ 
quency  range  where  wavelength  is  comparable  to  the  body  size.  One  technique 
which  holds  promise  for  use  in  the  resonance  region  is  impedance  loading  of  the 
body  surface.  This  technique  is  discussed  in  several  papers.  Other  papers  give 
results  demonstrating  the  effects  that  material  properties  have  in  changing  cross 
sections  in  the  resonant  region. 

Volume  II  of  the  collection  is  cla  .sified.  It  contains  the  papers  "Radar 
Absorbing  Materials  for  the  Resonance  Region  (U)"  by  K.AI,  Siegel  and  "Comm._nt 
on  Some  Radar  Camouflage  Problems  (U)"  by  W.F.  Bahret. 
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Foreword 

Objectives  of  tbe  ^wnposium  on  Modification  of  Electromagnetic  tvatlcrinji  f.ro>'. 

N?ction>  in  the  Hc-onant  llc^hm 


We  organized  this  conference  to  emphasize  resonant  region  scattering  because 
we  thought  there  were  significant  results  to  report.  Using  the  technical  ideas  anu 
techniques  discussed  June  4th,  we  believe  the  scattering  properties  of  important 
radar  target-  can  oe  significantly  modified  by  reactive  loading  recessed  into  Ok 
targets.  As  is  often  the  case  in  science,  the  experimental  findings  ore  out  in  ad¬ 
vance  of  theoretical  and  even  conceptual  explanations. 

We  were  gratified  by  the  high-quality  papers  and  audience  at  this  symposium. 

It  is  hoped  that  you  .will  now  look  at  this  record  critically  to  find  the  contradictions, 
mistakes,  and  anomalous  data.  A  sound  national  effort  on  this  theme  should  add 
another  tool  for  controlling  the  scattering  of  bodies. 

In  my  opinion  the  unsolved  problems  are  these:  What  other  general  shapes 
besides  the  electric  dipole  (rod)  or  string  of  dipoles  (multiply  loaded  rod)  w.ll  have 
a  minimum  total  scattering  cross  section  upon  application  of  suitable  localis  ed 
loading?  Montgomery1  hints  that  there  might  be  a  class  of  symmetrical  shapes 
with  such  properties.  What  is,  for  example,  the  suitable  loading  to  reduce  the 

total  cross  section  of  a  body  representable  by  both  electric:  and  magnetic  dipoles0 

*> 

What  is  the  significance  of  adding  loss  to  the  terminal  loacis  on  longer  rods"  (over 

\  / 2  in  length)  to  maximize  scatter  reduction?  How  does  one  make  the  (ransition 

from  a  frequency  range  where  specular  reflections  are  dominant  and  where  matched 
3 

dipoles’  or  ordinary  absorbing  materials  provide  minimum  back  scattering  to  the 
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resonant  region  where  slightly  lossy  and  lossless  reactive  loadings  appear  to  do  the 

4 

job?  Does  the  Conductron  paper' s  discussion  of  surface  impedance  loading  suitable 
for  resonant  shapes  answer  this  question,  or  is  the  situation  significantly  different 
for  terminal  loading  with  discrete  network  impedance  ports  on  the  body?  Can  the 
reactive  loading  method  be  used  to  redirect  scattering  in  the  forward  scatter  direc¬ 
tion,  for  example,  as  absorbers  tend  to  do,  rather  than  reduce  the  total  cross 
section? 

Why  do  the  values  of  inductive  reactive  loading  needed  to  minimize  scattering 
cross  section  decrease  appreciably  as  the  rod  is  made  thicker?  Why  is  the  band¬ 
width  better  with  dual  loading  than  with  single  (mid-point)  loading  on  rods  of  other¬ 
wise  similar  shapes?  Is  the  desired  reactive  loading  vs.  frequency  for  minimum 

g 

back  scatter  approximately  a  constant  value  as  the  data  show?  When  the  input 
impedances  vs.  frequency  characteristics  for  positions  of  loading  on  a  scatterer  are 
known,  can  they  be  synthesized  by  known  passive  networks  or  will  the  desired  fre¬ 
quency  response  violate  Foster's  reactance  theorem?  Do  the  external  modes  on  the 
object  (spherical  harmonics  in  the  case  of  a  sphere,  for  example)  dictate  the  proper 
positions  for  the  terminal  loadings  and  even  the  internal  cavity  structure  in  order 
to  realize  the  proper  phase  relations  of  the  surface  currents  to  produce  minimum, 
total,  or  differential  cross  section?  Why  does  the  bistatic  radar  cross  section  of  a 
rod  center  loaded  to  produce  minimum  monostatic  radar  scattering  for  all  attitudes 
increase  after  a  bistatic  angle  of  20  to  30  degrees?  Can  scattering  from  two  ortho¬ 
gonal  polarization  vectors  l  e  treated  independently  by  suitable  orthogonal  loading  on 
the  objects? 

When  the  volume  of  the  object  rather  than  its  surface  is  available,  what  values 

of  the  constitutive  parameters  will  augment,  reduce,  or  otherwise  modify  the 

7  Q 

scattering  of  resonant  objects?  The  work  of  Avco  (plasma  media)  and  Atlas  (di¬ 
electrics)  hints  at  unusual  possibilities  here. 

Can  reactive  loading  methods  be  successfully  combined  with  shaping  or  absorb¬ 
ing  techniques  to  cover  a  large  bandwidth  for  a  large  angular  region  of  observation? 

There  are  provisional  and  perhaps  correct  answers  to  all  these  questions. 
Probably  they  will  all  be  answerable  and  fall  like  a  row  of  dominoes  as  our  under¬ 
standing  of  the  resonant  region  scattering  and  the  influence  of  strategically  located 
terminal  loads  improves. 

In  spite  of  unanswered  questions,  some  very  useful  engineering  results  are 
available.  Thick  cargo-carrying  bodies  can  be  reduced  significantly  in  scatter 
cross  section  for  all  attitudes  over  the  peak  of  the  first  resonance  hump.  This 
loading  can  be  recessed  into  the  body  and  with  coaxial  or  lumped  network  impedance 
could  be  packaged  pretty  conveniently.  At  Ohio  State^  and  Boeing1  ®  good  strides  - 
are  being  made  to  combine  and  control  antenna  and  object  scattering  by  judicious 
use  cf  impedance  loading.  I  hope  the  prodding,  provocation,  and  discussions  oi'  this 
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conference  will  help  to  bolster  the  theoretical  picture  and  bring  useful  applications 
rather  soon. 

C.  J.  SLETTEN 
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THE  .MODIFICATION  OF  ELECTROMAGNETIC 
SCATTERING  CROSS  SECTIONS  IN  THE 

RESONANT  REGION 


I.  Some  Thoughts  on  Scattering  Cross  Sections  in  the 

Resonance  Region 

Richard  B.  Mack  and  Philipp  Blacksmith.  Jr. 

Microuave  IMi\-ics  Laboratory 
•\ir  Force  Cambridge  Research  Laboratories 
Bedford.  Mas-achusetts 


One  certain  way  to  start  an  argument  in  even  a  small  group  of  people  familiar 
with  electromagnetic  scattering  is  to  use  the  expression  "  resonance  region"  . 
Someone  will  immediately  object  that  for  most  shapes  there  is  no  resonance  in 
the  usual  circuit  sense  of  zero  susceptance  and  maximum  current.  Someone  else 
will  point  to  the  maxima  and  minima  in  curves  of  scattering  cross  section  vs. 
frequency,  and  then  the  discussion  really  begins.  Such  instances  simply  reaffirm 
what  most  of  us  readily  admit  -  that  when  the  characteristic  dimensions  of  a  tar¬ 
get  are  of  the  same  order  of  magnitude  as  the  wavelength,  the  electromagnetic 
properties  are  immersed  in  a  considerable  cloud  of  mystery. 

A  thin  dipole  is  by  far  the  simplest  of  such  targets,  for  it  has  essentially 
a  one -dimensional  line  current.  In  this  case,  the  measured  and  computed  back- 
scatter  cross  sections,  current  distributions,  and  driving-poim  odmittam  es  can 
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be  closely  correlated  to  give  a  clear  picture  of  its  behavior.  However,  as  soon 
as  currents  on  the  scatterer  are  permitted  to  have  additional  components,  such 
clarity  vanishes.  The  sphere,  for  example,  is  a  simple  three-dimensional  shape 
that  still  requires  four  to  seven  terms  of  the  Mie  solution  to  describe  its  proper¬ 
ties,  even  when  ka  is  only  1.  it  to  1.  5.  Already  there  is  some  difficulty  in  de¬ 
veloping  an  intuitive  feeling  for  the  results.  For  most  shapes,  including  the  very 
interesting  complex  ones,  even  this  kind  of  a  solution  is  not  presently  available 
when  the  body  dimensions  are  comparable  in  size  to  the  wavelength  of  the  imping¬ 
ing  wave.  If  the  scatterer  is  made  of  dielectric  or  lossy  material,  the  situation 
is  further  complicated.  The  result  is  that  our  ability  to  use  or  modify  target 
properties  in  the  resonance  region  is  severely  limited. 

Determining  the  scattering  cross  sections  from  an  interesting  shape  and 
learning  to  modify  them  are  the  first  steps  necessary  toward  an  ultimate  goal  of 
controlling  the  scattering  properties  -  developing  a  simple  positive  method  to 
make  the  cross  sections  appear  larger  or  negligibly  small  as  desired.  Alterna¬ 
tively,  if  energy  is  to  be  received  by  or  transmitted  from  the  scatterer,  the 
problem  is  to  place  antenna  structures  so  they  intercept  the  currents  with  maxi¬ 
mum  efficiency  or  excite  the  body  to  concentrate  radiated  energy  in  a  desired 
direction. 

Although  the  scattering  properties  of  a  target  are  usually  examined  by  them¬ 
selves  in  graphs  of  o  vs.  angle  or  ka,  their  principal  application  is  in  determin¬ 
ing  how  susceptible  a  target  may  be  to  detection  by  a  radar.  It  is  interesting  to 
examine  some  common  cross  sections  from  this  point  of  view.  The  radar  range 
equation  is 

O 

max 

where  RmaJJ,  P^,  G,  Smin,  \  are,  respectively,  the  maximum  range,  the  power 

transmitted,  antenna  gain,  minimum  detectable  signal,  and  the  wavelength,  tr  is 

the  backscatter  cross  section  of  the  target.  Note  that  P  ,  S  .  ,  G,  and  \  are 

t  min 

properties  of  the  radar  and  independent  of  the  target,  whereas  <he  backscatter 
cross  section  is  not  an  exclusive  property  of  the  target  because  it  is  a  function  of 
the  wavelength. 

The  obvious  and  most  important  conclusion  from  Eq.  (1)  is  that  large  changes 
in  care  necessary  to  produce  significant  changes  in  Rmax-  For  example,  a  12  db 
decrease  in  uonly  decreases  the  maximum  detectable  range  by  half;  a  40  db  de¬ 
crease  is  required  in  cr  to  reduce  R  by  00  percent.  This  is  illustrated  in  Fig- 

m  sx  2 

ure  1  for  a  radar  which  can  just  detect  a  target  of  1  in  cross  section  at  1000 
miles. 
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Figure  1.  Variation  of  Maximum  Range  With  Backscatter  Cross  Section.  Norm¬ 
alized  to  1000  miles  for  a  1  m-  target. 


In  Figure  2  the  backscatter  cross  sections  of  several  metal  targets  are  com¬ 
pared  as  a  function  of  k"a",  where  "a"  is  a  characteristic  dimension  for  each 
particular  target  and  is  listed  on  the  figure.  Data  for  the  sphere  is  theoretical; 
that  for  the  other  curves  is  experimental  1  “D  For  values  of  k"a"  greater  than 
7  or  8  ,  there  are  large  differences  between  the  values  of  <r  from  a  sphere  and 
those  from  a  flat  plate  at  normal  incidence  However,  when  k"a"  is  about  1.  1, 
the  cross  sections  from  all  but  the  dipole  are  within  a  -idb  range.  This  produces 
a  corresponding  difference  in  Rmax  of  about  ±12  percent  from  an  average  value. 

If  the  fiat  plate,  an  extreme  example,  is  excluded  there  is  a  considerable  range 
of  k"a"  over  which  the  remaining  curves  are  separated  by  not  more  than  12  db, 
which  corresponds  to  a  variation  of  2  to  1  in  Rmax  While  ibis  is  a  significant 
decrease  in  Rmax  it  could  be  compensated  by  existing  state-of-the-art  techniques 
For  example,  a  larger  antenna  producing  f>  db  additional  gain  would  be  sufficient 
to  restore  the  original  value  of  R 

max 

Also,  when  k"a"  is  near  1,  the  variation  of  backscatter  cross  section  with  tar¬ 
get  aspect  for  the  ditferent  shapes  is  very  similar  to  the  variation  observed  from  a 


dipole,  i  his  is  further  discussed  and  illustrated  in  the  paper  by  Sletten,  et  al. 

Consider  now  the  detection  of  a  class  of  targets  having  different  shapes  but 
approximately  the  same  size.  What  frequency  band  should  he  used  to  detect  the 
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targets  at  the  greatest  distances'’  Assume  Pt/Sm and  the  antenna  gain,  G,  to 
be  constant.  The  target  size  will  be  considered  constant  and  frequency  will  be 
the  variable.  Actually,  the  powers  and  sensitivities  which  are  available  make  the 
ratio  P  / S  -n  decrease  as  the  frequency  is  increased.  Assuming  the  antenna 
gain  to  be  constant  means  that  once  the  desired  frequency  is  determined,  the 
antenna  will  be  made  sufficiently  large  to  provide  the  required  gain.  At  lower 
frequencies  where  large  reflectors  become  impractical,  multiplate  techniques 
might,  be  used.  The  ideal  would  be  to  always  design  the  antenna  so  it  produces 
the  maximum  obtainable  gain  which,  from  considerations  of  mechanical  tolerances 
and  coherence,  is  currently  estimated  to  be  about  65  to  70  db.  Within  a  chosen 
frequency  band,  the  antenna  will  have  a  fixed  size  and  its  gain  will  exhibit  the 
usual  frequency  dependence.  This  is  not  considered  in  the  following. 

With  these  assumptions,  the  radar  range  equation  can  be  written  as 
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max 
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Figure  2.  Comparison  of  Scattering  Cross  Sections  for  Various  Shapes 


2  *> 

Thus,  unless  n  varies  at  least  as  i  or  1A~,  R„ .....  will  decrease  as  the  ire- 

nidA 

quency  is  increased.  In  the  Rayleigh  region  where  the  taroet  dimensions  are  mu-  h 
smaller  than  a  wavelength,  the  backscatter  cross  sections  ot  many  objects  ai  e 
known  to  vary  as  1  f\  .  With  such  targets  and  frequencies  R  changes  as  %  : 
where  f  is  the  frequency.  The  Rayleigh  backscatter  cross  section  oi  a  spheie, 
for  example,  is 


9a  ,6 
a  =  -jz~  (ka) 


R  =  K«  •—  — 
max  4  ^ 


5.  122 


where  D  =  sphere  diameter  =  2a,  k='Y:  ,  k^  = 


At  the  other  frequency  extreme,  the  optics  region  where  target  dimensions 
are  much  larger  than  a  wavelength,  the  backscatter  cross  sections  of  many  ob¬ 
jects  are  independent  of  frequency.  With  these  targets  and  frequencies,  R 

2  *  n*u\ 

changes  as  — ,  decreasing  as  the  frequency  is  increased.  The  optics  cross 
section  of  a  sphere,  for  example,  is 


a  =  ~  a 


/-  \A  A 
=  K(j)4D*  -s/I  or 


=  0.  S414  \'k. 


In  contrast,  a  flat  plate  which  is  much  larger  than  a  wavelength  gives 


< t  =  A  ■  ,  A  =  plate  area 

A 


K(4jtF  D 


square  plate  (D  =  edge  length) 
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R  =  K(4r)4  '•JTJT  D  circular  plate  {D  =  diameter)  (9) 

max  r 

so  that  Rmax  is  independent  of  frequency. 

When  k"a"  is  between  about  0.  2  and  10,  the  so-called  resonance  region, 
values  of  a/\“  can  be  taken  from  Figure  2.  Let  T  be  a  given  value  of  o/A“  in 
Figure  2.  Then 


R 


max 


KD 


T4  kx  ,  kx  =  \/D 


(10) 


and  the  usual  oscillations  of  o  vs.  k"  a"  are  modulated  by  the  dependence  on  X 
(or  1/f). 

Some  typical  curves  of  the  normalized  maximum  range,  »  vs-  k" a" 

are  given  in  Figure  3  for  the  backscatter  cross  sections  of  Figure  2.  Again,  "a" 
is  a  characteristic  length  for  the  target  as  given  in  the  table  of  Figure  2,  and  D 
may  he  considered  as  an  average  or  most  probable  target  size  equal  to  2  "  a" . 
Under  the  assumption  discussed  above  and  for  the  targets  considered.  Figure  3 
clearly  demonstrates  the  desirability  of  choosing  a  frequency  which  either  places 
k"  a"  between  about  0.  5  and  2  or  makes  k"  a"  as  low  as  is  compatible  with  the 
obtainable  antenna  sizes.  Since  for  many  targets  D  is  quite  large.  Figure  3  also 
emphasizes  the  need  for  physically  large  antennas  which  have  high  gain  at  the 
lower  frequencies. 

The  most  commonly  discussed  methods  of  modifying  the  cross  sections  of 
targets  are  shaping  and  the  use  of  absorbers.  Figures  2  and  3  indicate  that  shap¬ 
ing  is  effective  at  larger  values  of  k"  a"  but  is  not  very  effective  when  k"  a"  is 
near  1.  Likewise,  absorbers  are  effective  at  larger  k"ar  values,  but  at  smaller 
k"  a"  values  they  are  either  not  particularly  effective  or  are  bulky  at  lower  fre¬ 
quencies. 

Most  of  the  symposium  papers  are  concerned  with  a  third  method  of  modify¬ 
ing  the  cross  sections.  This  method  has  been  called  "reactive  loading"  and  is 
especially  applicable  to  targets  in  the  range  of  small  k"  a"  and  at  lower  frequen¬ 
cies  where  the  more  common  methods  lose  their  effectiveness.  Reactive  loading 
holds  the  further  promise  of  providing  a  simple  means  of  not  just  modifying  but 
of  actually  controlling  the  cross  sections  by  simply  closing  or  opening  a  switch. 

It  is  hoped  that  the  symposium  and  its  record  will  not  only  lead  to  improved 
techniques  for  controlling  scattering  cross  sections  in  the  resonance  region,  but 
will  also  lead  to  a  better  understanding  of  basic  scattering  properties  when  the 
object  dimensions  are  comparable  in  size  to  a  wavelength. 


1‘igure  3.  The  Normalized  Maximum  Range  for  Detection  of  Targets  of  Figure 
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II.  The  Minimization  of  the  Backscattering  of  a 
Cylinder  By  Central  Loading* 

Run-Mu  •  li*-n  urn!  \  l.n-ji,! 
Dll-  K.uiialinn  l  .ilmraliir\ 
Hu-  I  ni\i-r-u>  of  Mirliiaan 


Abstract 


A  theoretical  3nd  experimental  study  of  the  minimization  of  the  backscattering  of 
of  a  thin  cylinder  by  central  loading  is  presented.  The  induced  current  on  a  cen¬ 
trally  loaded  cylinder  illuminated  by  a  plane  wave  at  rormal  incidence  is  theoreti¬ 
cally  determined  and  experimentally  measured.  The  magnitude  and  phase  of  the 
induced  current  can  be  greatly  changed  bv  a  central  irr.pe  lame  The  optimum 
loading  to  achieve  zero  backscatter  in  the  broadside  direction  has  been  determined 
for  a  thin  cylinder  shorter  than  2  wavelengths.  The  optimum  central  impedance 
for  the  purpose  of  minimizing  the  broadside  hackscattering  from  a  thin  cylinder 
over  a  wide  range  of  frequencies  has  also  been  determined. 


‘‘This  paper  was  presented  at  the  Svmposium  on  Modification  of  Fleets  omagnetic 
Scattering  Cross  Section  in  the  Resonant  Region,  Air  force  Cambridge  Research 
Laboratories,  Bedford,  Massachusetts,  June  1,  1%’?.  The  research  in  this 
paper  was  supported  bv  Air  Force  Cambridge  Research  Laboratories  under 
Contract  AF  l«(f.28>-237*. 
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I.  1MH0DLCT10H 

The  first  known  use  of  reactive  loading  to  minimize  scattering  was  by  fams1 
who  applied  the  technique  to  metallic  posts  used  in  a  parallel  plate  pillbox  struc¬ 
ture.  The  idea  of  using  this  technique  to  decrease  the  radar  cross  section  of 

9 

objects  in  space  was  suggested  and  employed  by  Slettfn“  in  1950.  The  backscatter- 

ing  of  a  cylinder  with  and  without  a  central  load  has  be  en  studied  b\  several 
3  4  5 

authors.  ’  *  The  information  available  to  date  is  that  if  a  cylinder  of  near  a  half 
wavelength  is  loaded  with  a  high  reactive  impedance  at  its  center,  its  backscatter- 
ing  cross  section  can  be  significantly  reduced.  The  exact  way  in  which  the 
reactive  loading  behaves  and  the  optimum  method  of  loading  are  still  not  well 
understood.  It  is  the  purpose  of  this  paper  to  present  a  theoretical  and  an  experi¬ 
mental  study  of  the  subject  in  order  to  clarify  the  nature  of  the  loading  and  also  to 
present  an  optimum  loading  for  reducing  to  zero  the  broadside  backscattering  from 
a  cylinder. 

We  study  the  problem  from  the  viewpoint  of  the  induced  current  on  a  metallic 
object  illuminated  by  an  electromagnetic  wave,  since  the  induced  current  plays  the 
most  fundamental  role  in  the  scattering.  We  consider  a  perfectly  conducting  cylin¬ 
der  with  a  small  radius  and  a  length  shorter  than  2  wavelengths,  and  assume  that 
a  plane  wave  is  incident  broadside  on  the  cylinder.  A  current  is  induced  on  the 
cylinder  by  the  incident  plane  wave,  and  this  induced  curre  it  in  turn  produces  a 
scattered  electromagnetic  field.  If  an  impedance  is  added  at  the  center  of  the 
cylinder,  the  induced  current  and  the  scattered  field  are  modified.  There  are 
three  methods  by  which  the  addition  of  a  central  impedance  can  reduce  the  scattered 
field.  These  are:  (a)  by  reducing  the  magnitude  of  the  induced  current;  {b)  by 
reversing  the  phase  of  the  induced  current  over  some  part  of  the  cylinder;  and 
(c)  by  a  combination  of  (a)  and  (b).  The  third  method  is  the  most  effective  for 
reducing  the  scattering,  and  we  shall  show  that  with  central  loading  it  is  possible 
to  reduce  the  broadside  backscattering  to  zero. 

In  Section  2,  we  determine  theoretically  the  induced  current  on  the  cylinder 
with  central  loading  by  applying  an  integral  equation  method.  This  induced  current 
is  expressed  as  a  function  of  cylinder  dimensions  and  the  central  impedance.  In 
Section  3,  we  study  the  induced  current  on  a  cylinder  without  loading.  In  Section  4, 
the  induced  current  on  a  cylinder  loaded  with  an  infinite  impedance  is  found.  In 
Section  5,  we  study  the  induced  current  on  a  cylinder  with  a  resonant  length  for 
various  central  impedances.  In  Section  6,  the  induced  current  on  a  cylinder  with 
an  anti-resonant  length  is  studied  for  various  central  impedances.  Through  Sec¬ 
tions  3  to  6,  theoretical  and  experimental  results  are  compared.  In  Section  7,  the 
experimental  method  is  described.  In  Section  8,  we  obtain  the  optimum  impedance 
for  zero  broadside  backscattering  from  a  thin  cylinder.  The  optimum  impedance 
required  to  minimize  the  broadside  backscattering  from  a  thin  cylinder  over  a  wide 
range  of  frequencies  is  also  included. 
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In  this  paper,  attention  has  been  given  only  to  the  broadside  case,  since  the 
highest  return  is  in  this  direction.  We  have  limited  our  analysis  to  a  thin  cylinder 
in  the  interest  of  simplicity  in  developing  the  theory.  The  case  of  a  thicker  cylin¬ 
der  and  the  case  of  the  oblique  backscattering  will,  however,  be  studied  later.  The 
use  of  multiple  loading  on  a  cylinder  will  also  be  investigated  in  the  future. 


2.  INDUCED  CURRENT 

The  geometry  of  the  problem  is  as  shown  in  Figure  1.  A  cylinder  with  a  radius 
a  and  a  length  2h  is  assumed  to  be  perfectly  conducting.  A  plane  electromagnetic 
wave  with  the  E  field  parallel  to  the  axis  is  incident  normally  to  the  cylinder.  At 
the  center  of  the  cylinder  a  lumped  impedance  is  connected,  fhe  dimensions 
of  interest  are 


\  <  2h  <  2\ 


*  « 
o 


1 


where  X  is  the  wavelength  and  is  the  wave  number.  The  second  condition 
implies  that  the  cylinder  is  thin  and  that  we  can  assume  that  only  the  axia'>  current 
is  induced. 


6— *■ 

Incident 
EM  Wave 


l/4X<2h<2X 
PqO*  « I 


Figure  1. 


Cylinder  with  Central  Loading  Illuminated  by  an  Electromagnetic  Wave 
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2.1  Inteprel  Equation  for  the  Induced  Current  on  the  Cylinder 

In  order  to  determine  the  induced  current  on  the  cylinder,  we  apply  an  integral 
equation  method.  The  integral  equation  for  the  induced  current  is  derived  first. 
The  incident  tangential  electric  field  is  assumed  to  be 


(1) 


where  Eq  is  constant  along  the  cylinder. 

The  tangential  electric  field  at  the  cylinder  surface,  maintained  by  the  current 
and  charge  on  the  cylinder,  is 


Z 


ao 

3z 


-  jw  A 


z 


(2) 


where  0  is  the  scalar  potential  maintained  by  the  charge,  and  A  is  the  tangential 
component  of  the  vector  potential  maintained  by  the  current.  By  using  the  Lorentz 
condition 


V.i 

Ea  +  Ein  =  0  for  6  <z  <  h  and  -h  <  z  <  -6 .  <u 

Z  2  —  “  —  — 

At  the  gap,  the  electric  field  is  continuous.  That  is 

Ea  +  Ein  =  Eg  =  Z,  I  6(z)  ,  for  -6  <  z  <  6  •  (3) 

z  i  z  Lo  —  — 


Actually,  a  single  equation  can  be  obtained  by  combining  Eqs.  (7)  and  (C)  and  using 
the  result  with  Eq.  (4)  over  the  whole  length  of  the  cylinder  to  obtain 


A  +  £2 
z  o 


A 

z 


ZL*o 


for  -h  <  z  <  h 


(9) 


Equation  (9)  is  an  inhomogeneous  differential  equation  for  A  .  The  general  solu¬ 
tion  is  the  sum  of  the  complimentary  function  and  a  particular  integral  as  follows: 


a 

"z 


COS  + 


C9  sin 


/3qz  +  n  (z; 


■] 


(10) 


In  Eq.  (10)  v  is  l/\V0e0»  C1  and  C9  are  arbitrary  constants.  0(z)  is  a  particu¬ 
lar  integral  and  can  be  expressed  as 


Q  (z) 


y  [Eo<s)  - z 


T  1JS) 

L  o 


3 


sir.  ^o(z  -  s)  ds 


E 

=  -ip’  0-  -  cos  M) 
p  o 


1 

2 


sin 


(11) 


If  E  is  assumed  to  be  constant  along  the  cylinder,  C9  can  be  proved  to  be 
o  * 

zero  from  the  symmetry.  Equation  (10)  then  becomes 

E 

Az(z)  =  T~o  [C1  C0S  V  +  T0  (1  '  COS  Poz)  -  \  ZL  lo  sin  *0^] 

for  -  h  <  z  <  h  . 

From  Eq.  (12),  Cj  can  be  expressed  in  terms  of  A?(h)  as 


(12) 


Ka(z,  z')  = 


exp  VT-z'  )2  +  a2  J 


!(z-z')2  +  a2 


and  Iz(z’)  is  the  induced  current  on  the  cylinder. 

An  integral  equation  for  the  induced  current  on  the  cylinder  >s  then  obtained 
by  equating  Eqs.  (14)  and  (15)  as  follows: 


^  lz(z>)  Kd(z*  z)  dz'  =  T~  se"  ^oh  (j  vc  Az(h)  ~li)  (C°S  V  "  C0S  ^oh) 

-h  °  V  °  '  _ 


+  2  ZLTosin 


where  f  =  120jt  and  (18)  is  valid  for  -h>z>h.  It  is  noted  that  A  (h)  and  I  in 
o  z  o 

the  right-hand  side  of  Eq.  (18)  are  the  functions  of  Iz(z)  and  are  still  unknown. 


2.2  Solution  for  tin-  Induced  (.urr"nl  on  a  (  \Iindcr 

The  right-hand  side  of  Eq.  (18)  suggests  a  form  for  the  solution  V  fgil  as 


I  fa)  =  C  (cos  ,3  z  -  cos  3  h)  +  C  sin  3  (h  -  lz|  ) 
i.  CO  o  s  o 


(19) 


It  is  then  reasonable  to  divide  Eq.  (18)  into  tv.o  parts  as  follows: 


it 

{  (cos  fr'T.'  -  cos  0oh)  Kjfz,  z')  dz! 


c  ,»  'o 
-h 


-j4”  !  Ec\ 

—  sec  «oh  (i  vo  Az0>>  -  i;  (:OS  30z  -  cos  3Qh) 


(20) 
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Cg  ^  sin  0Q(h  -  jz'|)  Kd(z,  z)  dz' 


(21) 


-h 

-i2r  11% 

=  —  sec  In  ZT  I  sin  3rt(h  -  J z ; )  . 

i  .  O  i  o  u 

o 

Equations  (20)  and  (21)  are  valid  for  -h<  z<h  and  also  well  matched  at  the  end 
points,  z  =  ±h.  To  find  the  constants,  Cc  and  Cg,  we  can  match  both  sides  of 
Eqs.  (20)  and  (21)  at  the  center  of  the  cylinder,  z  =  0. 

By  setting  z  =  0  in  Eq.  (20),  C  is  determined  as 


Cc  =  rrr sec  V  I-1  vo  Az(h)  ~  -f~)  (1  - cos  ^0h) 

*o  cd  \  ^o 


(22) 


where 


u 

Tcd  =  1*  (cos  V'  "  C0S  ^oh)  Kd(0,  z  )  dz'  ■ 


(23) 


By  setting  z  =  0  in  Eq.  (21),  Cg  is  determined  as 


-J  2tt 

C  =  r— 7p —  sec  (3  h  Z,  I  sin  0  h 
s  ?o  Tsd  °  L  °  ° 


(24) 


where 


(25) 


h 

Tgd  =  \  sin  3Q(h  -  }z' ')  KjfO,  z  )  dz''  . 

-h 

Substitution  of  Eqs.  (22)  and  (24)  in  Eq.  (19)  gives 


Iz(z) 


-j4n 

r 

o 


Az(h)  -  j  (sec  3q h  -  1)  (cos  3^  -  cos  S^n) 


+  —  Zt  I  tan  3  h  sin  3  (h  -  |zi) 

~Asd  L  °  °  ° 


(26) 


Equation  (26)  is  not  the  final  form  of  the  solution  of  I  (z)  because  A  (h)  and  I 

z  z  o 

are  still  unknown. 

Since  IQ  =  Iz(z  =  0),  we  can  express  IQ  in  terms  of  Az(h)  by  letting  z  =  0  in 
(26),  and  Eq.  (26)  can  then  be  arranged  <±s 


Iz(z)  = 


~^T~  M  vo  Az(h,~3^)  j^1'  (cos  V  "  cos  5oh)  sin  ^ 


(27) 


where 


M'  =  =r—  (sec  3  h  -  1) 
cd  ° 


(28) 


N 


-Z^  tan  30b(sec  3Qh  +  cos  3Qh  -  2) 

T  ,  ZT  tan  3  h  sin  3h-j60T  ,T  , 
cd  L  o  o  J  cd  sd 


(29) 


Equation  (27)  still  has  one  unknown,  Az(h),  and  to  determine  this  we  do  the  follow¬ 
ing.  From  the  definition  of  vector  potential,  we  have 


Az(h) 


h 

V  L(z')  K  (h,  z')  dz' 
«./  Z  » 

-h 


(30) 


where  Kg(h,  z')  is  defined  in  Eq.  (17).  If  Eq.  (27)  is  substituted  in  Eq.  (30),  we 
obtain 


Vh>  =  Vo 


M'T  +  N'T 
_ ca  sa 

1  -  M'T  -  N'T 

ca  sa 


(31) 


where 


h 

Tca  =  \  (cos  *Q  z'  -  cos  JQh)  Ka'h,  z'  )  dz' 
-h 


(3  2) 


h 

Tsa  =  \  sin  *Vh  "  ^z  ^  KaQi,z')dz'  (33) 

-h 

Now  a  final  form  of  the  solution  of  Iz(z)  can  be  obtained  if  Eq.  (31)  is  substituted 
in  Eq.  (27).  After  rearrangement,  the  final  solution  for  Iz(z)  can  be  summarized 
as  follows: 


iE 


Iz(z) 


1 


30*o  VCOS^oh-MTca-XT 


sa / L- 


M(cos3  z 
o 


■  cos  3  h) 
o 


Xsin  .5  (h-  z 
o 


where 


M 


(1  -  cos  *0h) 


(35) 


-ZT  sin  *  h(l  -  cos  3  h)2 

M  =  _ L _ o - o -  (3fi) 

Tod  ZL  s”  0oh  -  i  60Tcd  Tsd  cos  3oh 

and  TC(j,  Ts(J,  Tca,  Tga  are  defined  in  Eq.  (23),  (25),  (32),  and  (33). 

Equation  (34)  gives  a  complete  expression  for  the  induced  current  on  a  cylinder 

with  a  central  load,  ZT  ,  when  illuminated  bv  a  constant  electric  field,  E  ,  at 

i-/  o 

broadside.  The  accuracy  of  Eq.  (34)  is  high  and  its  form  is  simple  and  suitable 
for  the  further  development  oi  theory.  It  is  checked  experimentally  in  later 
sections. 

As  a  matter  of  completeness  and  convenience,  integrals  TC(j,  Tg(j,  Tca,  and 
Tga  are  expressed  in  terms  of  better-known  integrals  as  follows: 

Tcd  =  Ca(h>  0)  “  Ca(!l'  h)  “  cos^oh  j\(h’ (,)  "  Ea(h'  h)]  (37) 

Tsd  =  sin  3oh  [Ca(h>  0)  "  Ca(h’h]  “  cos  'V1  [Sa(h>  0)  "  Vh,h^]  (38) 


Tca  =  Ca^h’h  *  “  cos  30h 


(30) 
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T  =  sin  5  h  C  (h,h)  -  cos  3_h  SJh,h) 
ScL  c  a.  O  c*. 

where 

h 

Ca(h,  0)  =  \  cos /^z7  Ka(0,  z' )  dz' 

-h 

h 

Ca(h,h)  =  ^  cos0oz  Ka(h,z/)dz/ 

-h 


h 

Ea(h,0)  =  ^Ka(0,z')dz' 

-h 

h 

Ea(h,h)  =  jKa(h,z')dz' 

-h 

h 

S  (h,  0)  =  \  sin  fl  z'  K  (0,  z')  dz7 

3  J  O  <x 

-h 

h 

S  (h,  h)  =  C  sin  8  z'  K  (h,  z')  dz7 
a  j  o  a 

-h 


(40) 


(41) 


(42) 


(43) 


(44) 


(45) 


(46) 


The  integrals  of  Eqs.  (41)  to  (46)  can  be  calculated  by  using  a  digital  computer. 


3.  INDUCED  CURRENT  ON  \  CYLINDER  WITHOUT  CENTRAL  LOADING 

The  first  and  the  simplest  case  to  be  studied  is  a  cylinder  without  loading. 
The  induced  current  on  the  cylinder  can  be  found  from  Eq.  (34)  by  letting  =  0. 
That  is 

jE„  /  1  -  cos  h 

r  \  -  O  I  — - — ^ - 

zw  30  (T  ,  +  T  )  cos  fl  h  -  T 

^o  \  cd  ca  ^o  ca 


(cos  SQz  -  cos  0oh) 


With  Eqs.  (37)  and  (39),  Mz)  can  be  expressed  as 

(1  -  cos  3  h)(eos  3  z  -  cos  3  h) 
_ o _ o _ o _ 

Ca(h,  0)  cos  3q h  -  Ea(h,  n)  cos2  3Q h  -  cyh.h)-  Ea(h,h)cos  5Qh 

(47) 

The  distribution  of  the  induced  current  along  the  cylinder  in  this  case  is  a  shifted 
cosine  curve.  The  maximum  induced  current  occurs  at 


iE 

I  (z)  =  — S- 
30  5 


z  =  0 
z  =  X/2 


for  3  h  <  — - 


for 


¥<  i>„h  <  \ 


(48) 


and  is  given  by 


I  (0) 


iE 
_ _ o 

30  3 


(1  -  cos  J  hi" 


C  (h,0)cos3  h-E  (h,0)cos  3  h-C  (h,  h)*-E  (h,h)cos3  h 

1—3  O  a  03  3  O  — * 


(49) 


or 


m  JEo 

r  .  2,  ,  “i 

-  sin  3  h 
o 

\2)  "3O0o 

Ca(h,0)cos£oh  -  Ea(h,0)  cos23Qh  -  Ca(h,h) Ea(h,h)  cos  3Q h 

Theoretical  and  experimental  results  of  1(0)  as  a  function  of  h  are  compared 
and  shown  graphically  in  Figure  2.  The  theoretical  results  are  calculated  from 
Eq.  (49)  with  a  computer.  The  experimental  results  are  obtained  by  measuring 
the  induced  current  at  the  center  of  a  3/16-in. -radius,  variable  length  cylinder; 
the  cylinder  length  can  be  varied  between  X/4  and  2X.  The  cylinder  is  illuminated 
by  a  plane  wave  with  a  frequency  of  1.038  kMc.  The  agreement  between  theory 
and  experiment  is  good  except  at  a  point  near  h  =  0.7  X  where  a  resonant  peak 
occurs.  The  discrepancy  at  this  point  may  be  due  to  theoretical  error  or  to  the 
fact  that  the  incident  electromagnetic  wave  in  the  experiment  is  not  constant  along 
the  cylinder  for  lengths  greater  than  X.  Since  the  general  behavior  between  X/4 
to  2X  is  quite  well  predicted  by  theory  and  confirmed  by  experiment,  we  have  not 
attempted  to  minimize  the  above  disagreement. 

The  current  distributions  along  the  cylinder  of  various  lengths  are  shown 
graphically  in  Figure  3.  Theoretical  and  experimental  results  are  in  good 
agreement. 
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Figure  2.  I  ( 0)  vs  h/A.  when  3Qa  =  0.11,  f  =  1.038  kAIc,  ZL  =  0. 

The  results  obtained  in  this  section  are  not  unexpected.  The  induced  current 
reaches  a  resonant  peak  on  a  cylinder  with  a  length  of  0.43X  or  a  length  of  1.4X 
for  0oa  =  0.11.  Obviously,  these  current  peaks  imply  large  backscattering  signals. 
We  propose  to  eliminate  these  current  peaks  by  the  use  of  suitable  impedance 
loading. 

4.  INDl  (  ED  Cl  RRENT  ON  A  C \  LINDER  WITH  AN  INFINITE  IMPEDANCE  AS  ITS  CENTRAL  LOAD 

The  second  case  to  be  studied  is  a  cylinder  with  an  infinite  impedance  loaded 
at  its  center.  Theoretically,  the  induced  current  can  be  obtained  from  Eq.  (34)  by 
setting  -*■  .  Experimentally,  an  infinite  impedance  is  approximated  by  a 
coaxial  cavity  tuned  at  its  anti -resonance  position.  This  coaxial  cavity  is  built 
inside  of  the  cylinder  as  described  m  a  later  section.  The  induced  current  on  the 
loaded  cylinder  is  then  measured  by  a  small  probe. 


THEORY 


h  =  0.258X 


Figure  3.  Distribution  of  Iz(z)  ’.then  =  o. 


When  Zjj-  *  ,  from  Eqs.  (35)  ana  (36) 
-(l-cos3oh) 


M 


sin  i)  h 
o 


for  3C(h  *  n"  . 


(51 ) 


The  substitution  of  Eq.  (51)  in  Eq.  (34)  gives 


iz(Z) 


-  i  E 
_ o 

30  3 


(1  -  cos  3oh)  { sin  3q  I  z  |  -  sin  ?Qh  x  sin  3o(h  -  !  z  | )] 
rsin0oh  (Ca(h,0)  -  (2 -cos  30n)  Cyh.h)  -  cos  ,3Qh  iyh.0)  +  ryh.h] 


■  (1  -  cos  3Qh)2  •c;a(h,  h) 


(52) 


The  induced  current  in  this  case  is  zero  at  the  center  of  the  eylindei,  and  is  dis¬ 
tributed  along  the  cylinder  as  a  combined  curve  of  a  sine  and  a  shifted  sine  curte. 
The  maximum  induced  current  occurs  at 


for  (3  h<  2tr 
o 


z  = 


h 

2  » 


(53) 
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and  is  given  by 


-jgc 

30  3 


3  h  ,  jS  n  . 

2  sin  -y  (l  -  cos  -2-J  (1  -  cot  3n) 


sin  3  h  [C  (h,0)  -  {2  -  cos  S  hi  C  (h,h)  -  cos  ,<?  h  E  (h,0)  +  E  Jh.hlf 


-U -cos 3  h)“  S  (h, h) 

O  a 


(54) 

The  theoretical  value  of  )  as  a  function  of  h/A.  is  shown  graphically  in 
Figure  4  in  comparison  with  an  experimental  curve.  The  experimental  results  are 
obtained  with  the  coaxial  cavity  set  near  its  anti-resonance  position  The  agree¬ 
ment  between  theory  and  experiment  is  not  very  good.  The  discrepancy  is  due  to 
the  fact  that  it  is  impossible  to  obtain  experimentally  an  infinite  impedance  from  a 
coaxial  cavity  structure.  This  reasoning  is  supported  by  the  fact  thi  t  there  is 
closer  agreement  between  theory  and  experiment  for  Z^  =  j2000  fi. 


Figure  4.  Iz(h/2)  vs  h/\  when  3Qa  =  0.11,  f  =  1.088  kMc,  ZL  =  large. 
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The  distribution  os'  the  induced  current  on  three  cylinders  of  different  lengths 
and  loaded  with  infinite  impedance*  is  shown  graphically  in  Figure  5.  The  agree¬ 
ment  between  theory  and  experiment  is  good  for  these  cases. 


THEORY 


1.0  - 

0.8  - 


•  •  • } 

a  a  a  ^EXPERIMENT 

o  o  oj 


2 1  =  6.34cm 
28  =  0.127cm 


Jz|/X 


Figure  5.  Distribution  of  I  (z)  for  h  =  0.2 13X, 
h  =  0.259A,  and  h  =  0.326X  when  =  °°  . 


The  important  information  obtained  with  infinite  impedance  at  the  center  of 
the  cylinder  is  as  follows:  The  induced  currents  at  resonant  lengths,  namely 
2h  =  0.43A.  and  2h  =  1.4X,  are  greatly  reduced  by  this  loading.  However,  the 
induced  current  appears  to  have  a  peak  when  2h  =  0.9A  with  this  loading.  This 
current  peak  should  be  suppressed  if  a  low  scattering  over  a  wide  frequency  range 
is  desired.  We  then  conclude  that  an  infinite  impedance,  or  a  very  high  impedance, 
for  the  central  loading  is  not  an  optimum  from  the  viewpoint  of  minimizing  the 
scattering  over  a  wide  range  of  frequencies.  We  shall  seek  an  optimum  loading  in 
later  sections. 


2*1 


3.  INDl  CEO  Cl  ROHM  0\  A  CYLINDER  OF  M  AR  RESONANT  LENGTH  SH  ii  \  ARIOl  CENTRAL 
IMPEDANCE^ 

Ln  this  section  the  effect  of  the  central  impedance  on  the  induced  current  of 
the  cylinder  with  a  resonant  length  is  studied.  For  experimental  convenience,  we 
choose  the  following  specific  case: 

a  =  0.0173A 

2h  =  0.43A 

ZL  =  JXL  • 

The  last  condition  restricts  the  central  impedance  to  be  purely  reactive,  since 
only  a  reactive  impedance  is  obtained  experimentally  from  a  coaxial  cavity. 

With  tne  above  conditions,  the  theoretical  value  of  the  induced  current  can  be 
expressed  as 


iE^  H (cos 3  z  -  0.  216)  +  sin  (77. 5*  -  |s  z|) 

I  (z)  =  ^rf-M  - 2 - ^ - - - 0 - 

2  Ju^o  j  0.215  -yj  (0.218  -  j  0.25) 

where 

.  -0.765  XL 

=  0.955  XL  -  24.6  * 


(55) 


(56) 


The  relative  magnitude  of  Iz(z)  is  calculated  and  shown  graphically  in  Figures 
6a  3nd  6b  for  the  following  values  of  Zj j  0,  »  ,  -jlOOOfl,  -jSOOft,  -j600£2, 
j  160012,  j 80012,  j  600S,  andj400f2.  In  these  curves,  we  observe  the  following 
facts: 


(1)  When  Z^  =  0  (no  loading),  the  induced  current  is  very  large  and  distributes 
along  the  cylinder  as  a  shifted  cosine  curve. 

(2)  When  Zj  =  00 ,  the  magnitude  of  the  induced  current  is  greatly  reduced 
from  the  value  for  Z^  =  0  and  the  distribution  of  the  induced  current  becomes 
double  humped  with  a  null  at  the  center. 

(3)  When  Z^  is  capacitive  and  finite,  the  induced  current  is  smaller  than  the 
case  of  Z^  =  0  but  larger  than  for  Z^  =  *>. 

(4)  When  is  inductive  and  finite,  the  magnitude  of  the  induced  current  is 
smaller  than  the  case  of  Z^  =  <=°  and  the  induced  current  starts  to  have  three  loops 
along  the  cylinder;  it  is  of  interest  to  note  that  the  phase  of  the  current  at  the 
center  loop  is  reversed. 
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The  most  important  and  significant  information  obtained  from  Figures  6a  and 
6b  is  that  when  Zp.  is  inductive  and  of  some  finite  value,  the  magnitude  of 
the  induced  current  is  reduced  to  a  value  smaller  than  the  case  of  Z^  =  * 
and  the  phase  of  the  induced  current  is  reversed  at  the  center  part  of  the 
cylinder.  It  is,  therefore,  possible  to  reduce  the  broadside  backscattering 
from  a  cylinder  to  zero  by  proper./  adjusting  the  value  of  Zj  .  In  effect, 
the  optimum  impedance  for  zero  broadside  backscattering  from  a  cylinder  of 
this  size  is  inductive  and  it  has  a  small  resistive  component,  as  we  shall 
see  in  a  later  section. 

To  obtain  an  experimental  verification  to  the  tneoretical  results  shown  in 
Figures  6a  and  Sb,  we  measured  the  induced  current  along  a  cylinder  of  this 
specified  dimension  and  with  various  cavity  lengths.  The  experimental  results 
are  shown  in  Figures  7a  and  7b.  We  obtained  a  family  of  curves  which  close¬ 
ly  resemble  the  theoretical  curves  shown  in  Figures  6a  and  6b.  When  the 
cavity  length  (total  length)  is  longer  than  6.  2  cm,  the  impedance  of  the  cavity 
is  capacitive;  it  becomes  inductive  for  cavity  lengths  shorter  than  6.  2  cm.  It 
It  should  be  noted  that  the  effective  cavity  length  is  greater  than  these  values 
since  it  is  loaded  with  a  dielectric  material  for  which  e  =  4.  0  .  The  approxi¬ 
mate  value  of  the  cavity  impedance  is  calculated  by  using  a  standard  impedance 
formula  for  a  transmission  line  and  assuming  that  a  capacitance  of  0.  4  p/if 
is  shunted  across  the  gap  at  the  center  of  the  cylinder. 

The  comparison  between  theory  and  experiment  is  made  in  Figure  8  v>here 
theoretical  curves  for  Z^  =  *  j  8001?  ,  *  ,  and  j  BOOQ  are  shown.  These 
curves  are  compared  with  experimental  results  for  £  =3.32  cm,  3.  10  cm, 
and  2.  01  cm  where  £  is  the  half  length  of  the  coaxial  cavity.  The  agree¬ 
ment  between  theory  and  experiment  is  very  good.  This  indicates  also  that 
the  calculated  value  of  the  cavity  impedance  is  quite  close  to  the  corresponding 
theoretical  impedance. 


6.  INDUCED  CURRENT  ON  \  CYLINDER  OF  NEAR  \NTI-RESON\NT  LENGTH  WITH  \  UUOl  > 
CENTRAL  IMPEDANCE." 

In  this  section  we  study  the  effect  of  the  central  impedance  on  the  induced 
current  of  a  cylinder  with  an  anti -resonant  length.  The  dimensions  of  the  cylin¬ 
der  and  the  central  impedance  are  chosen  to  be 

a  =  0.  0173A 
2h  =  0.  9  A 

H.  ’  J\- 


Relative  Amplitude  of  |Iz(z)j 


The  relative  magnitude  of  I^fz)  *s  calculated  and  shown  graphically  in  Figuse 
9a  and  9b  for  values  of  as  follows:  *,  -jl60Wi,  -j80(42f  -jSOOii,  j  lb'HUi, 

j  80UQ  and  j  600f2.  This  family  of  curves  is  quite  different  from  those  of  the  pre¬ 
ceding  section.  Although  a  purely  reactive  impedance  reduces  the  magnitude  of 
the  induced  current  and  tends  to  reverse  the  phase  of  the  induced  current,  U  is 
not  possible  to  reduce  the  broadside  backscattering  to  zero  because  current  nulls 
do  not  occur  in  this  case.  Actually,  an  optimum  impedance  for  zero  broadside 
backscattering  from  a  cylinder  of  this  size  should  have  a  large  resistive  compo¬ 
nent,  as  is  shown  in  a  later  section. 

The  experimental  results  for  the  induced  current  on  the  above  cylinder  with 
various  cavity  lengths  are  summarized  in  Figure  10.  The  general  shapes  of  the 
experimental  curves  are  similar  to  those  in  the  theoretical  curves.  The  compari¬ 
son  between  theory  and  experiment  is  made  in  Figure  11.  Three  typical  theoretical 
curves  are  compared  with  the  corresponding  experimental  results.  The  agreement 
between  theory  and  experiment  is  good  but  not  as  good  as  in  the  case  of  a  shoiter 
cylinder. 


7.  THE  EXPERIMENT 

In  Figure  12  is  shown  the  block  diagram  of  the  equipment  used  for  the  current 
measurement  on  the  reactivelv  loaded  cylinders.  In  the  experiment  the  cylindei 
is  illuminated  at  broadside  by  a  plane  wave  of  1.088 kMc  from  a  L-band  horn 
antenna  with  the  electric  field  vector  polarized  m  the  direction  parallel  to  the 
cylinder.  A  conventional  probing  method  with  a  small  current  probe  was  employed 
to  measure  the  induced  current  amplitude  on  the  cylinder. 


RELATIVE  AMPLITUDE  OF  |t2U)|  RELATIVE  AMPL 


Figure  9.  Curre.it  Distribution,  I7(?;)  ,  on  a  Cylinder  of  h  =  0.45 \  ,  a  =  0.  0173  \ 
for  Different  Central  Load  (Theoretical). 
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h  =  0.45AQ 

j  2a  =  0. 92  cm 

,  A  =  27. 57  cm 

o 

IP  4-  2S  =  0.127  cm 


_ theory  for  Z^  =  0 

&  &  experiment  for  Z  =  0  ; 

- theory  for  =  jSGOfJ  ; 

o  c  o  experiment  for  Z  =2.7  cm  \ 

- theory  for  Z^=  —  j  300^2  jg 

•  0  •  experiment  for  Z=  3.67  cm 
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Figure  11.  Current  Distribution  along  a  Cylinder  as  a  Function  of  Central  Load 
for  h  =  0.45  Aq  (Theoretical  and  Experimental). 


The  coaxial  line  leading  from  the  probe  was  covered  with  radar  absorbing 
material  (RAM)  and  oriented  perpendicular  to  the  E  field  to  minimize  the  influ¬ 
ence  of  its  presence-  The  measurement  area  was  lined  with  RAM  to  reduce 
unwanted  reflections. 

Figure  13  shows  partially  disassembled  components  of  the  loaded  cylinder. 
The  diameter  of  the  cylinder  is  about  0.95  cm  and  its  length  can  be  changed  from 
10em(h  =  0.182A)  to  51.29em(h  =  0.93A)  bv  the  combination  of  center  and  enci  pieces 
of  different  lengths.  The  center  sections  of  the  cylinder  contain  a  symmetric 
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Figure  12.  Block  Diagram  of  the  Experimental  Setup 


Figure  13.  Cylinder  with  Coaxial  Cavity 


coaxial  cavity  with  its  input  gap  at  the  center  of  the  cvlmdei.  R\  varving  the  cavity 
length,  various  input  leaetunces  that  represent  'he  centi.r  loads  of  the  cylinder  ate 
obtained.  The  coaxial  cavity  is  filled  with  a  dielectric  m  oidei  to  reduce  the  re¬ 
quired  cavity  length.  The  dielectric  used  is  ^  stvrac.ist  with  a  dielectric  constant 
of  4  and  a  loss  tangent  of  o.'mmsI. 
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8.  OPTIMl  M  IMPEDANCE  FOR  ZERO  R RO  ADSIDE  It \Ck?C ATTERIM,  FROM  A  THIN  C\  LINDER 


In  the  preceding  sections  we  have  concentrated  on  finding  the  induced  current 
on  a  cylinder  as  a  function  of  the  cylinder  dimension  and  the  central  impedance. 
We  now  seek  to  find  an  optimum  impedance  for  zero  broadside  baekscattering 
from  a  thin  cylinder. 

The  induced  current  has  been  found  in  Eqs.  (34)  and  (36).  The  scattered  tield 
due  to  this  induced  current  can  be  found  as  follows: 

The  vector  potential  maintained  by  the  induced  current  on  the  cylinder  in  the 
far  zone  of  the  cylinder  is 


A_  = 


u  j  E 
^o  J  o 


4-  30  0O  l  cos/?oh  -^Tca  “NTsa 
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\  (cos  - 


exp  (-j  3  R) 


cos  3  h)  D 

O  n 


■  dz 


-h 

h 


+  N  ^  sin  3  Ah  -  |  z  | ) 


exp  (-j^cR) 
R 


dz 


-h 


(59) 


where 

R  =  R  -zcos  0  -  distance  between  a  point  on  the  cylinder  and  the 
observation  point. 

The  scattered  electric  field  in  the  far  zone  of  the  cylinder  can  be  obtained  as 


Ej  =  -j «  Aa  =  j  us  A ^  sin  9  . 


(60) 


The  corresponding  Poynting  vector  is 


pS  ’  wo  M2 


(61) 


The  scattered  field  in  the  broadside  direction  can  be  obtained  by  letting 


o  ' 


9  =  90°  and  R  =  R 


(62) 
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We  then  obtain 


E*  (ft  =  90°)  =  -r-  E 


H30K0> 

R, 


M(sin3  h-j  hcos  i  h)-*-\(l  -cosd.h) 
_ o  o  o  o 


cos (3  h  -  MT  • 
o  ca 


XT 


sa 


(63) 


and 


PS(ft  =9C°) 


f  0  9 

°  KK 

O  O 


M(sin/3  h-3  hcosd  h)+X(l  -  cos  3  h) 
_ o  o _ o _  o 

cosa  h  -  MT  -XT, 
o  ca  sa 


2 


(64) 


The  Poynting  power  density  of  the  scattered  field  in  the  broadside  direction  is 
thus  expressed  as  a  function  of  the  central  impedance  because  N  in  Eq.  (64)  is  a 
function  of  Z^,  as  expressed  in  Eq.  (36). 

To  minimize  the  broadside  backscattering  to  zero,  we  can  simply  make 
Ps(ft  =  90°)  equal  to  zero.  That  is. 


X 

\f 


sin  £Qh  -  3Qh  cos  dQh 
1  -  cos  0Qh 


(65) 


Using  Eqs.  (35)  and  (36),  Eq.  (65)  can  be  rewritten  as 

Z,  sind  h(l  -cos  3  It)  sind  h  -  3  h  cos  3  h 

J- _ o _ _c _  _  _ _^o _ o  o 

ZL  sin23Qh  -  j  60Tsd  cos  3Qh  1  '  cos  (?oh 


(66) 


By  solving  for  Z^  in  Eq.  (66),  we  obtain  the  optimum  central  impedance  for  zero 
broadside  backscattering  as 


~j  60Tsd(1  -g0hcotftoh) 

2  cos  ah  -  2  +  ft  h  sin  3  h 
o  o  o 


(67) 


where 

Tsd  =  sin  ,3oh  (Ca(h'  0)  '  Ca(h' h)  1  ‘  cos  Soh  1  Sa(h'  0)  "  Sa(h'  h)  ! 
as  expressed  in  Eq.  (38). 

Equation  (67)  gi’'es  the  complete  expression  for  the  optimum  central  impedance 
for  zero  broadside  backscattering  from  a  thin  cylinder.  This  optimum  central 
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impedance  is  a  function  of  the  cylinder  dimension  and  its  numerical  value  can  be 
calculated  quite  easily.  The  expression  should  prove  useful  in  practical  design 
We  consider  a  cylinder  with  a  =  0.  0173  \  and  calculate  j~Zjl 

I—  J  ^ 

as  a  function  of  the  cylinder  length  h.  The  result  is  shown  granhically  in  Figure 
14.  From  this  figure  we  observe  the  following: 

(1)  In  general,  the  optimum  central  impedance  for  zero  broadside  backseat  - 
tering  should  have  both  the  resistive  and  the  reactive  component. 

(2)  For  a  cylinder  shorter  than  1  wavelength  (or  h<  0.5X),  the  optimum 
impedance  is  inductive  and  requires  a  resistive  component. 

(3)  For  a  cylinder  longer  than  1  wavelength  (or  n  >  0.5X),  the  optimum 
impedance  is  inductive  or  capacitive  but  requires  a  negative  resistive  component. 


Figure  14.  Optimum  Central  Impedance,  [  ZjJq  ,  for  Zero  Backscattering  from 
a  Cylinder  of  Radius  a=0.  0173  X  ,  as  a  Function  of  Cylinder  Length 

These  results  indicate  that  for  a  cylinder  shorter  than  1  wavelength,  a  passive 
impedance  loaded  at  the  center  of  the  cylinder  can  reduce  the  broadside  backscat¬ 
tering  to  zero.  However,  an  active  impedance  is  required  to  obtain  zero  broadside 
backscattering  from  a  cylinder  longer  than  1  wavelength.  Ot  course,  an  active 
impedance  may  not  be  needed  if  the  cylinder  is  loaded  at  two  points. 

We  also  consider  the  case  of  a  thicker  cylinder  with  a  =  0.0517X.  The  opti¬ 
mum  impedance  for  zero  broadside  backscattering  is  obtained  as  a  function  of  the 


cylinder  length  h  and  shown  graphical^  m  Figure  15;  -nd  j_-  seen  to  be  similai 
to  that  in  Figure  14.  The  resistive  component  remains  lek  tivelv  constant  •.*.  j»b 
respect  to  the  change  of  the  cylinder  thickness,  but  re  ?. -active  component  is 
reduced  almost  by  a  factor  of  two  for  the  thicker  cylinder. 


400t 


Figure  15.  Optimum  Central  Impedance,  [  Zl  lo  >  for  Zero  Backscattering  from 
a  Cylinder  of  Radius  a=0.  0517  A  ,  as  a  Function  of  Cylinder  Length 


To  determine  bandwidth  characteristics  of  this  technique,  we  calculate  the 
optimum  impedance  for  the  purpose  of  minimizing  the  broadside  backscattering 
from  a  fixed  cylinder  over  a  wide  range  of  frequencies.  We  considei  a  cylinder 
with  h  =  4  cm,  a  =  0.476  cm  to  cover  a  frequency  range  from  1  kMc  to  3  kMc.  The 
frequency  characteristic  of  the  optimum  impedance  is  shown  graphically  in  Figure 
16.  Within  this  frequency  range  the  optimum  impedance  is  inductive  and  requires 
a  resistive  component.  This  impedance  appears  to  be  obtainable  by  a  simple  net¬ 
work  synthesis.  It  is  also  noted  that  to  cover  a  wider  range  of  frequency,  an 
active  impedance  is  needed  to  reduce  the  broadside  backscattering  to  zero. 

A  final  remark  is  necessary  at  this  point.  Although  we  have  exclusively 
studied  the  broadside  illumination,  the  optimum  loading  obtained  in  this  study  is 
also  effective  for  the  case  of  the  oblique  illumination.  To  support  this  statement 
the  experimental  measurement  of  the  induced  curicnt  by  a  plane  wave  incident 
from  different  angles  are  presented  in  Figures  17  and  13.  In  these  figures  we 
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frequency  (Kmc) 


Figure  16.  Optimum  Central  Impedance,  l  Z^Jo  ,  for  Zero  Backscattering  from  a 
Cylinder  with  h=4  cm  ,  a=0.  476  cm  Over  the  Frequency  Range  of  IkMc  to  3kMc 
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Figure  17.  Current  Distributions  on  a  Center  Loaded  Cylinder  Due  to 
Plane  Waves  of  Different  Incident  Angle. 

h  =  0. 215A0 
26  =  0.050  in. 


The  current  scale  is  Relative  to  the  Maximum  Current  (*.0)  on  an 
Unloaded  Cylinder  with  h  =  0.215Ao,  <i>  -  0. 
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Figure  13.  Current  Distributions  on  a  Center  Load  Cylinder  Due  to 
Fiane  Waves  of  Different  Incident  Angies 

h  =  U.215.\0 
2o  =  0.050  in 
*r  =  * 

L  =  6.21  cm 

The  Current  Scale  is  Relative  to  the  Maximum  Current  ( ’  .0)  on  an 
Unloaded  Cylinder  with  h  =  0.215.\o,  o  =  C. 

observe  that  the  only  change  on  the  induced  current  for  the  oblique  illumination  .? 
the  reduction  in  its  magnitude.  The  distribution  of  the  induced  current  is  practi¬ 
cally  unchanged  from  the  case  of  the  broadside  illumination.  Additional  study  on 
the  oblique  illumination  case  will  be  made  in  the  future. 


9.  CONCI.I  SION 

We  have  studied  the  effect  of  central  loadings  on  the  induced  current  on  a  thin 
cylinder  that  is  illuminated  by  a  plane  wave  at  normal  incidence.  It  is  found  that 
both  the  induced  current  and  the  scattered  field  can  be  greatly  modified  by  a  cen¬ 
tral  load.  Through  Iheoi  etical  analysis  and  experimental  measurement,  we  have 
good  understanding  about  the  exact  nature  of  the  loading  for  the  reduction  ol  the 
scattering. 

We  have  also  obtained  the  optimum  central  loading  to  eliminate  ihe  broadside 
backscattering  from  a  thin  cylinder.  A  relatively  simple  formula  for  the  optimum 
impedance  is  presented,  and  it  should  prove  useful  in  practical  design  situations. 
Although  an  impedance  which  can  easily  be  obtained  by  a  simple  r.etwoik  synthesis 
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is  sufficient  to  reduce  to  zero  the  backscattering  f.  om  a  cylinder  shorter  than  1 
wavelength,  it  appears  that  an  active  impedance  is  required  to  eliminate  the  back- 
scattering  from  a  cylinder  longer  than  1  wavelength.  This  difficulty  can  probably 
be  overcome  by  loading  the  cylinder  at  two  points  with  two  passive  impedances. 

An  investigation  of  this  possibility  is  in  progress. 
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Abstract 


An  impedance-loading  technique  for  backscatter  reduction  of  long  thin  bodies  is 
discussed  and  experimental  verification  of  the  technique  presented.  It  is  shown  that 
a  significant  reduction  in  the  backscatter  is  possible  while  providing  a  receiving 
capability.  The  body  is  thus  converted  into  an  antenna  that  has  a  high  ratio  of 
receiving  cross  section  to  backscattering  cross  section  under  conditions  ox  low 
backscatter. 


1.  INTH0DTCT10N 

Backscattering  from  a  long  thin  body  in  the  region  of  longitudinal  incidence  has 
been  evaluated  by  Peters*’^  in  terms  of  a  traveling-wave  antenna  mode  excited  on 
the  body  by  an  incident  wave.  In  this  paper,  an  impedance-loading  technique  is 
described  which  effectively  results  in  cancelling  the  backscatter  field  produced  by 
such  a  mode.  The  body  is  considered  as  an  asymmetrical  dipole  antenna,  and  the 
value  o't  load  impedance  necessary  to  accomplish  the  cancellation  of  backscatter  is 
determined.  Using  the  Snnth-chart  representation  of  Kennaugh  and  Green'  for 


antenna  scattering,  the  backscattering  cross  section  is  evaluated  for  any  load  imped¬ 
ance  at  the  frequency  under  consideration. 

An  inciden.  wave  with  the  plane  of  polarization  parallel  to  the  axis  of  the  body 
excites  a  traveling  wave  on  the  body.  This  wave  is  reflected  from  the  rear  point, 
producing  a  secondary  traveling  wave  on  the  body.  It  is  this  backward  wave  that 
results  in  the  large  backscattering  cross  section.  Elimination  of  the  reflected  wave 
from  the  rear  point  would  affect  a  large  reduction  in  backscattering  over  a  range  of 
aspect  angles  in  which  the  travelling-wave  mode  is  the  primary  mode.  It  is  thus 
postulated  that  a  set  of  antenna  terminals  be  constructed  near  the  rear  point  so  as  to 
absorb  the  incident  wave  without  reflection,  or  to  reflect  part  of  the  received  signal 
with  'he  proper  uhase  and  amplitude  to  cancel  That  which  is  not  coupled  to  the  ter¬ 
minals.  In  this  paper,  such  a  technique  is  shown  to  be  possible. 


2.  IMPEDANCE  LOADING 

In  order  to  evaluate  the  backscattering  characteristics  of  a  long  thin  body  with 

an  impedance  load  near  the  rear  point,  the  general  concept  of  antenna  scattering 

should  be  considered.  The  electric  field  scattered  by  an  antenna  as  a  function  of 

3 

load  impedance  is  given  by 


E(Zl)  =  E(0) 


ZL2.4 

zl  +  za 


1(0)  et 


(1) 


where 

=  Rl  +  jXL  is  the  load  impedance, 

ZA  =  Ra  +  jXA  is  the  antenna  impedance, 

E(0)  is  the  scattered  field  of  the  antenna  with  Z^  =  0, 

—  T 

E  is  the  electric  field  transmitted  by  the  antenna  when  excited  by  a  unit 
voltage  source,  and 

1(0)  is  the  current  through  the  load  Z^  =  0, 

3 

This  equation  may  also  be  written  as 


n 


where 


Z  A  io  the  complex  conjugate  of  Z  A, 


E(Z^)  is  the  structural  scattered  field. 


—  4AI  — T 

E*  '  =  1(0)  E  is  the  antenna-mode  scattered  field, 

w  A 


i  -  Z 


r  = 


N 

L  . 


1  +  z 


tt  is  the  reflection  coefficient,  and 


N  ZL^XA 

Z^  =  - ^ - —  is  the  normalized  load  impedance. 

A 


If  the  condition 


r  EAM=  -e(z'a)  (3) 

is  fulfilled,  then  the  scattered  field  E(Z^ )  is  zero.  Ir.  order  to  accomplish  this 
over  some  range  of  aspect  angles,  one  should  have 


gAM  E(Za) 

I:™  |e(za)| 


f-1) 


over  that  range.  The  load  impedance  is  then  adjusted  to  provide  the  proper  phase 
and  amplitude  of  the  reflection  coefficient  so  as  to  fulfill  Eq.  (3). 

The  long  thin  body  with  a  set  of  terminals  near  the  rear  point  may  be  con¬ 
sidered  as  an  asymmetrical  dipole  antenna.  For  a  high  degree  of  asymmentry, 
the  primary  mode  in  the  region  of  longitudinal  incidence  is  the  traveling-wave 
mode.  When  short-circuited,  this  asymmetrical  dipole  reduces  to  the  long  thin 
body  considered  by  Peters1,  From  Eq.  (1)  the  structural  scattering  E (Z  A)  may 
be  written  as 


E(Z*  )  =  S(0)  -  EAM  e~''2° 


(5) 


where  o  is  the  phase  angle  of  the  antenna  impedance. 


One  can  also  write 


1 


E(0)  =  yEj 


(6) 
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.v here  >  is  the  reflection  coefficient  jt  the  end  of  the  both  and  E^  is  the  field 
radiated  b>  the  ‘rave  ling -wave  antenna  mode  excited  on  the  short-circuited  body. 
Thus  Eq.  <  4)  should  hold  for  the  lorn,  thin  body  in  the  region  of  longitudinal  inci¬ 
dence.  Considering  only  the  components  of  Eqs.  (3),  (5),  and  (6),  the  reflection 
coefficient  necessary  for  cancellation  then  becomes 


(7) 


AM  AM 

if  E^  =  E'  and  y  is  constant  as  a  function  of  aspect  angle,  then  it  should 
be  possible  to  find  a  load  impedance  that  would  result  in  cancelling  the  structural 
scattering.  If  Eq  W  t  E^^,  then  the  ioad  impedance  would  need  to  be  varied  as  a 
function  of  aspect  angle  to  affect  the  cancellation. 

For  the  case  considered,  the  terminals  are  some  finite  distance  from  the  end 
AM  AM 

of  the  bodv  so  that  Eq  and  E*  are  not  identical.  However,  if  the  distance  is 
sufficiently  small  they  are  nearly  the  same.  In  this  case  it  should  be  possible  to 
find  a  load  impedance  that  would  result  in  considerable  reduction  in  backscattering 
over  some  range  of  aspect  angles. 


3.  EXPERIMENTAL  EVALl  ATION 

In  order  to  verify  that  effective  cancellation  of  the  structural  scattering  is 
possible,  an  experimental  evaluation  was  made.  The  test  model,  shown  in  Figure 
1,  consists  of  a  30-degree  ogive  with  a  coaxial-fed  probe  extending  from  the  rear 
of  the  model.  Backscatter  and  receiving  cross  sections  were  measured  at  a  fre¬ 
quency  of  2000  Me  with  the  plane  of  polarization  parallel  to  the  axis  of  the  body. 

For  these  measurements,  the  model  was  supported  on  a  polyurethane  foam  tower. 

A  small  coaxial  cable  was  brought  out  from  the  model  and  connected  to  an  adjustable 
load  impedance,  a  power  meter,  and  to  other  equipment  located  at  the  base  of  the 
tower  For  the  polarization  used,  the  coaxial  cable  was  perpendicular  to  the  inci¬ 
dent  electric  field,  hence  its  effect  on  the  measurements  was  sufficiently  small. 

The  model  was  positioned  at  an  aspect  angle  corresponding  to  the  peak  of  the 
main  lobe,  and  the  load  impedance  was  adjusted  for  minimum  backscatter.  In  this 
manner,  the  peak  of  the  main  lobe  was  reduced  by  30  db.  The  entire  backscatter 
pattern  was  then  measured,  the  results  of  which  are  shown  in  Figure  2.  It  is  seen 
that  an  effective  backscatter  reduction  of  20  db  is  accomplished  over  the  main  lobe. 
The  30-db  reduction  of  the  peak  of  the  main  lobe  is  not  maintained  for  other  aspect 
angles  because  of  the  difference  between  the  two  antenna  modes  and  E^‘,  as 


10  LOGt 


Figure  2.  Backscatter  Reduc¬ 
tion  Characteristics 


Figure  1.  Backscatter  Model 
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Figure  3.  Receiving 
Characteristics 
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indicated  by  Eq.  (7).  In  order  to  maintain  nearly  complete  cancellation,  the  load 
impedance  would  need  to  be  adjusted  as  a  function  of  aspect  angle. 

The  receiving  characteristics  of  the  model  were  evaluated  under  conditions  of 
minimum  backscatter  and  maximum  received  power.  These  characteristics  are 
show'n  in  Figure  3,  and  indicate  that  the  received  power  under  conditions  of  mini¬ 
mum  backscatter  is  3  db  less  than  the  maximum  received  power  obtained.  Thus, 
one-half  of  the  received  power  is  reflected  from  the  terminals  with  the  proper 
phase  and  amplitude  to  effectively  cancel  'he  structural  scattering.  For  the  condi¬ 
tion  of  maximum  received  power,  a  backscatter  reduction  of  only  5  db  was  accom¬ 
plished.  The  measured  values  of  antenna  impedance  and  the  load  impedance  neces¬ 
sary  for  minimum  backscatter  were  found  to  be 

=  82  -  jl5  ohms 


and 


ZT  =  68  -  j98  ohms. 

3 

As  shown  by  Kennaugh  and  Green,  Eq.  (2)  may  be  represented  on  a  Smith-chart 

impedance  diagram.  This  allows  the  backscattered  field  to  be  evaluated  for  any  load 

impedance  by  experi mentally  evaluating  the  backscattering  cross  section  for  two  parti  - 

cular  values  of  load  impedance.  From  ZA  and  ZT  givenabove,  the  normalized  load 
N  ^ 

impedance  ZL  corresponding  to  minimum  backscatter  is  found  to  be  0.74-j  1.3  . 
Circles  about  this  point  represent  constant  backscatter,  as  shown  in  Figure  4.  Addi¬ 
tional  measurements  allow  calibration  of  these  circles  in  backscattering  cross  section. 

The  reflection  coefficient  j  F  |  is  found  to  be  0.  61  so  that  { F  |  2  =  0.  37.  Thus, 
only  about  one-third  of  the  received  power  should  be  required  to  cancel  the  struc¬ 
tural  scattering.  It  is  expected  that  if  various  losses  can  be  eliminated  in  the 
experimental  apparatus,  the  increased  power  corresponding  to  this  reflection 
coefficient  could  be  measured. 

4.  CONCLUSION 

It  has  been  shown  that  backscattering  due  to  the  traveling-wave  mode  on  a  long 
thin  body  can  effectively  be  cancelled  by  an  impedance-loading  technique.  Further¬ 
more,  the  cancellation  is  accomplished  with  only  a  small  reduction  in  received 
power.  Using  the  Smith-chart  representation  of  antenna  scattering,  the  backscat¬ 
tering  cross  section  is  evaluated  for  any  load  impedance  at  the  frequency  under 
consideration. 


Figure  4.  Backseattering  From  an  Asymmetric  Dipole 


Although  the  impedance-loading  technique  as  described  ior  the  long  thin  body 
is  rather  frequency  sensitive,  it  may  be  possible  to  extend  the  bandwidth  by  select  - 
ing  the  proper  thickness  and  asymmetry  of  the  scattering  body.  This  coupled  irh 
a  feedback  system  for  control  ot  the  load  impedance  may  result  in  a  large  back- 
scatter  reduction  over  a  wide  bandwidth  and  aide  range  of  aspect  angles.  Prelim¬ 
inary  experiments  in  this  area  indicate  hopeful  results. 
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Abstract 


Backscattering  from  a  monopole  (linear  scatterer)  can  be  considerably  dimin¬ 
ished  and  almost  suppressed  if  it  is  grounded  in  a  cavity  and  the  protruding  length 
of  the  monopole  does  not  exceed  0.  45  A  .  The  optimum  depth  of  the  cavity  is  about 
0.25  A. 

The  good  agreement  between  theory  and  experiment  justifies  the  rather  bold 
theoretical  assumptions  made  to  avoid  involved  mathematics  and  implies  their  use¬ 
fulness  in  similar  problems. 


1.  INTRODUCTION 

So  far  as  we  know,  backscattering  from  a  linear  scatterer  grounded  in  a  cavity 
has  not  previously  been  studied,  either  theoretically  or  experimentally.  The  only 
work  coming  close  to  ours  is  the  experimental  and  theoretical  work  on  unloaded 
antennas  cited  by  R.  W.  P.  Kang  in  The  Theory  of  Linear  Antennas,  Harvard  Univer¬ 
sity,  1956,  pp.  508,  509,516. 


In  our  initial  (1956)  study  of  backseatteiing  from  a  grounded  scatterer,  we 
restricted  ourselves  to  the  case  of  a  monopole  projecting  from  the  center  of  a 
cylindrical  cavity  in  ^  plane,  perfectly  conducting  screen  (Figure  1);  the  monopole 
was  excited  by  a  wave  traveling  along  the  surface  S  with  an  electric  vector  E 
parallel  to  the  scatterer.  The  restriction  to  cavity -loaded  monopoles  was  actually 
unnecessary,  since  the  formulas  obtained  led  easily  to  a  generalization  for  mono- 
poles  loaded  with  an  arbitrary  impedance.  Although  we  demonstrate-  this  possibility 
in  the  text,  we  did  not  change  our  study  to  the  general  case  because  it  had  been 
stimulated  by,  and  the  experiments  conducted  for,  the  specific  question  of  cavity- 
loaded  monopoles. 

A  further  restriction  we  set  was  to  consider  only  unmodulated,  harmonic, 
electromagnetic  waves.  In  a  sequel  to  this  report  we  will.deal  with  cases  of 
oblique -incidence  waves  and  modulated  or  pulsed  electromagnetic  waves. 


Figure  1.  Scheme  of  Scatterer  Arrangement 


As  in  almost  all  antenna  problems,  an  exact  solution  with  present  techniques  is 
difficult  if  not  impossible.;  the  more  or  less  dubious  initial  assumptions  of  a  theore¬ 
tical  treatment  must  be  justified  by  experimental  substantiation.  To  keepthe 


mathematical  treatment  relatively  simple,  we  made  rather  far-fetched  assumptions;5 
however,  we  carefully  checked  the  theory  by  conducting  very  reliable  experiments. 
These  experiments  were  conducted  by  Mr.  William  Kearns  after  some  exploratory 
measurements  by  Mr.  H.  Poehler.  The  experimental  setup  and  method  of  meas¬ 
urements  are  described  in  Appendix  A. 


2.  THEORY 


Our  theoretical  approach  to  the  pioblem  (see  Figure  1)  was  based  on  the  idea 
that  (1)  the  part  of  the  monopole  outside  the  cavity  could  be  regarded  as  a  lossy 
line  of  length  l ,  having  a  char.. cteristic  impedance  Z  with  an  additional  voltage 
dV  per  element  imposed  by  the  incoming  electric  field  E  ;  and  (2)  the  part  of  the 
monopole  within  the  cavity  could  be  regarded  as  a  lossless  line  of  length  6,  having  a 
characteristic  impedance  Zc,  grounded  at  ihe  bottom  of  the  cavity.  Strictly,  this  holds 
only  for  6  »  r ^  (r  =  radius  of  cavity),  but  to  simplify  the  mathematics  we  ex¬ 
tended  it  to  cases  where  6  «  r  . 

c 

Length  x  and  x  on  lines  l  and  6  are  measured  toward  the  screen  surface 
S  .  Their  currents  are  considered  to  be  positive  in  the  direction  measured.  The 
electrical  field  strength  E  of  the  incident  wave  is  taken  to  be  normal  to  the  sur¬ 
face  and  positive  in  the  outward  direction.  The  fieiJ  strength  of  the  incident  wave 
is  arbitrarily  consider  ed  to  be  zero  in  the  cavity. 

We  then  get  the  two  basic  relations: 


9V 

0x 


9L 

3x 


RI  .  L  |L  _  E(t,  x)  , 


GV-cf  ■ 


(1) 


where,  as  customary,  V  is  voltage,  R  is  resistance,  I  is  current,  G  is  shunt 
conductance,  L  is  inductance,  and  C  is  capacitance,  all  per  unit  length. 

If  we  suppose  that  the  incoming  plane  wave  is  incident  parallel  to  the  screen 
surface  S  ,  then  E  along  the  scatterer  is  a  function  of  only  the  time  t  . 

We  subject  Eq.  (1)  to  a  Laplace  transform  with  respect  to  time  and  then  to 
another  with  respect  to  length.  This  simplifies  the  rather  involved  analysis  and 


jJ { 

An  essential  disagreement  between  experimental  and  theoretical  results  would 
demand  a  "second-  or  third -order"  closer  theoretical  approach.  How  immensely 
such  a  closer  approach  complicates  the  formulas  is  seen  by  comparing  our  formu¬ 
las  given  in  King's  book  for  an  even  easier  problem. 
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gives  the  final  formula:* 


<2 


Tf  -  tanh  (r£t 


1+  ~  tanh  (rc6)  2  tanh 
1  +  tanh  (T  6)  tanh  (Ti) 

jL  C 


(2a) 


where  £  represents  the  backscattered  energy 

i 

r  =  [(R  +  uL)  (G+juC)]  2 


<7  _  TR+iwL, 
Z  " 


1 

l5 


and  i'  and  Z  are  the  analogous  quantities  for  the  cavity  parameters  R  ,  L  , 
c  c  *  c  c 

Gc,  and  Cc«  Specialized  for  6=0, 


~  |  I\f  -  tanh  (ri)J  ,  =  0  for  1=0  only. 


(2b) 


Equation  (2b)  describes  the  behavior  of  the  unloaded  monopole  directly  connected 
with  a  conducting  plane  and  also,  of  course,  the  backscattering  of  a  dipole  of 
length  21  in  freespace. 

If  t«X,  then  2 tanh  (Tf)  goes  to  1.  Therefore,  everything  after  the  minus 
sign  of  Ea.  (2a)  becomes  a  constant  that  will  be  small  compared  with  Ti .  Thus 
&  ~  |  | 2  ~  i  2.  The  backscattered  energy  increases  as  the  second  power  of  the 

ler^th  of  a  very  long  monopole. 

If  1«\,  then  2 tanh  ™  =  tan  (Tt )  and  £~Jri  -  tan  (Tl)\ 2  -  |(I\£ )3|  2~1  6. 

The  backscattered  energy  increases  as  the  sixth  power  of  the  length  of  a  very  short 
monopole. 

Unfortunately,  tne  evaluation  of  Eqs.  (2a)  and  (2b;  is  so  complicated  [even  the 
seemingly  simpler  formula  (2b)]  that  it  was  not  feasible  to  carry  out  calculations 
for  general- properties.  Equation  (2a)  was,  therefore,  calculated  with  a  computer 
for  certain  combinations  of  the  parameters  so  as  to  obtain  families  of  theoretical 
curves  for  comparison  with  the  experimental  ones.  The  parameters  Zc/Z,  r, 
and  r  had  to  be  chosen  as  close  as  possible  to  the  experimental  conditions,  but 
with  simplifications  necessary  to  keep  the  calculations  within  reasonable  limits. 

*The  complete  analysis  is  given  in  the  report  AFCRL-63-355  "Theoretical  and 
Experimental  Investigation  of  Backscattering  From  a  Cavity  loaded  Monopole" 
by  Werner  W.  Gerbes  and  William  J.  Kearns. 
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The  complex  propagation  constant  for  the  upper  part  of  the  monopole,  T=g-H/?  , 
where  a  is  the  attenuation  constant  and  p  is  the  phase  constant,  is  naturally  only 
approximately  a  constant.  But  wr-  can  set  p  equal  lo  that  of  frecspace,  so  that 
p  =  2 r/A  =  05 /c  .  The  attenuation  constant  a  is  related  to  the  radiation  resistance 
of  the  monopole.  To  avoid  complicated  calculations,  v/e  considered  a  a  parameter 
and  calculated  a  family  of  curves  for  every  depth  of  the  cavity,  setting  T  =  a  +  ia/c  . 

In  the  complex  propagation  constant  for  the  lower  part  of  the  monopole 
we  set  «c  =  0  because  the  radiation  of  this  part  can  be  neglected.  The  phase  con¬ 
stant  was  again  taken  as  approximately  equal  to  that  of  freespace.  Therefore, 


t  and  r  =i  it  . 
c  c  c 


For  Zc  we  used  the  formula  for  capacitance  of  a  coaxial  cable,  and  for  Z  we 
used  the  formula  for  capacitance  of  an  elongated  rotational  ellipsoid.  The  charac¬ 
teristic  resistance  is  inversely  proportional  to  the  capacitance  per  unit  length  C  . 

Thus  Z  /Z  =  C/C  .  Now 
c  c 


C 


c 


C  = 


Therefore,  with  f  /  2  r  =  p  , 

m  r 


A Ft 

arcosh  p 


(3) 
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where  r  is  the  radius  of  the  monopole,  r  is  the  radius  of  the  cavity,  and  £ 
m  c 

is  the  upper  length  of  the  monopole.  Evidently,  Z£/Z  is  a  function  of  £  and 
further  complicates  the  numerical  calculation. 

For  monopoles  of  small  length,  if  p  —  1  ,  Zc/Z  —  0/0  .  Ordinary  evaluation 
by  differentiation  gives 


a  constant  with  respect  to  I  .  For  monopoles  that  are  very  long,  if  p  »  1  ,  a 
short  calculation  gives  us 


showing  that  for  very  long  monopoles  Eq.  (2a)  changes  to  Eq.  (2b),  the  equation 
for  a  monopole  without  a  cavity.  This  means  the  cavity  has  no  effect  on  the  back- 
scattering  of  very  long  monopoles.  To  simplify  calculation  of  the  curves  in 
Figures  2  and  3  we  used  Zc/Z  =  1  as  a  convenient  mean  value,  which  automatically 
precludes  use  of  the  formula  for  long  monopoles. 

Basically,  there  are  two  ways  of  deriving  the  attenuation  constant  a  : 

(1)  directly,  from  the  radiation  resistance  of  a  monopole  divided  by  the  length; 

(2)  indirectly,  from  measurements  of  a  monopole  projecting  from  a  plane  surface 
without  a  cavity.  But  both  ways,  owing  to  the  assumptions  on  which  our  whole 
theory  had  to  be  based,  lead  to  more  involved  calculations  without  assuring  any 
greater  accuracy.  We  therefore  arbitrarily  chose  6  and  a  as  parameters  for 
the  calculation  of  the  families  of  theoretical  curves,  and  compared  them  with  the 
families  of  curves  measured  for  the  parameters  6  and  rm  .  All  measurements 
and  calculations  were  based  on  six  values  of  the  parameter  6  listed  as  follows, 
in  inches,  and  as  related  to  the  wavelength  used  (33  mm): 


^curve^ 

0 

1 

2 

3 

4 

5 

6  (in. ) 

0.  00 

0.  05 

0.  10 

0.  15 

0.  20 

0.  25 

a/x 

0.  00 

0.  04 

0.  08 

0.  12 

0.  16 

0.  20 
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It  may  be  appropriate  to  mention  at  this  time  that  the  theoretical  and  experi¬ 
mental  sets  not  only  match  substantially  in  shape,  but  even  coincide  quantitatively 
more  than  we  had  a  right  to  expect  from  the  somewhat  far-fetched  assumptions  of 
the  theory. 

The  theoretical  set  of  curves  of  backscattered  energy  as  a  function  of  the  pro¬ 
truding  length  l  of  the  monopole  for  a  =  0.  05  (  which  means  a  somewhat  low 
radiation  resistance  or  a  high  characteristic  impedance)  is  shown  in  Figure  2. 

The  curves  rise  with  the  sixth  power  to  a  first  sharp  maximum,  then  drop  to  show 
an  area  of  low  but  nearly  constant  backscattering  followed  by  an  additional  slight 
depression.  They  then  rise  to  a  second  maximum  and  drop  to  a  second  depression 
that  is  not  nearly  as  low  as  the  first.  This  behavior  repeats  itself  over  and  over.* 

Obviously,  the  first  depression  after  the  first  maximum  is  the  only  area  that 
exhibits  sufficient  backscattering  suppression.  With  a  recess  0.  20  X  deep,  the 
unsuppressed  backscattering  is  apparently  only  2  percent  of  the  maximum  scatter¬ 
ing  of  the  unrecessed  monopole. 

Figure  3  shows  the  theoretical  set  for  a  =  0.  10  ,  which  means  monopoles 
with  somewhat  higher  radiation  resistances  than  shown  in  Figure  2.  The  curves 
present  the  same  characteristics  as  the  curves  of  the  previous  set,  but  naturally 
not  so  pronounced  since  a  higher  radiation  resistance  means  more  strongly  damped 
circuits  with  smaller  resonance  effects.  In  this  case,  also,  a  recess  of  at  least 
0.  20  X  suppresses  about  98  percent  of  the  maximum  scattering  of  the  unrecessed 
monopole. 

Figures  2  and  3  indicate  that  for  satisfactory  scatter  suppression  the  protru¬ 
ding  length  of  the  recessed  monopole  must  be  less  than  0.  45  X  .  If  the  monopole 
protrudes  more  than  0.  45  X  ,  backscattering  cannot  be  sufficiently  suppi  essed  by 
means  of  a  cavity. 

These  theoretical  results  are  thoroughly  confirmed  by  the  experiments.  Fig¬ 
ure  4  shows  the  experimental  set  for  a  rod  1/16  in.  thick.  The  shape  of  the  curves, 
the  shift,  and  the  areas  of  possible  scatter  suppression  are  the  same  as  in  the 
theoretical  sets,  except  that  the  first  maxima  are  relatively  lower  with  respect  to 
the  second  maxima  than  in  the  theoretical  sets  (probably  due  to  the  choice  of 
Z  / Z  =  1).  Figure  5  shows  the  experimental  set  for  a  rod  1/32  in.  thick.  This 
set  also  shows  good  agreement  with  the  theory. 

Both  experimental  and  theorelicax  curves  show  a  defini.e  shift  of  the  maxima 
of  the  curves  with  increasing  depth  of  the  recess.  Although  it  appears  that  the 
positions  of  the  maxima  would  coincide  if  the  curves  were  plotted  versus  the  sum 
of  the  protruding  and  recessed  lengths,  closer  investigation  reveals  that  they  would 
still  not  exactly  coincide. 

*Our  simple  theory  obviously  yields  the  same  characteristics  as  those  produced  by 
the  more  complicated  formulas  cited  in  King's  book. 
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In  Eq,  (2a)  the  parameters  of  the  cavity  appear  only  in  the  grouped  form; 

Z£  tanh  (1^6),  which  represents  the  impedance  of  the  recessed  part.  Thus,  since 
the  recess  acts  as  an  impedance,  backscattering  could  be  suppressed  just  as  well 
by  grounding  the  monopole  over  a  specified  impedance.  Consequently,  if  in  Eq. 
(2a)  we  substitute  a  general  impedance  St  for  the  term  Zf  tanh  (T  o)  ,  v/e  obtain 
the  formula  for  the  general  case  of  backscattering  from  a  monopole  grounded  over 
an  impedance: 


e. 


r  £  -  tanh  (Ti ) 


1  +  Y  2  ta!lh  ^f 

1  +  ^  t jnh  (ri ) 


(2b) 


In  this  form  the  final  formula  also  covers  the  case  of  backscattering  from  a 
center -loaded  rod  in  space. 

Two  theoretical  curves  and  two  experimental  curves  for  each  of  the  basic 
parameters  are  compared  in  Figures  6  to  11.  The  shape,  and  even  the  size,  of 
the  curves  demonstrates  the  surprisingly  good  agreement  between  theory  and 
experiment. 

The  results  of  experiments  conducted  for  some  recesses  deeper  than  assumed 
for  theoretical  evaluation  are  given  in  Figures  12  and  13.  These  indicate  that  sup¬ 
pression  of  backscattering  can  be  achieved  with  recesses  deeper  than  0.  20  X  ,  but 
that  the  frequency  range  of  suppression  becomes  smaller  with  increased  depth. 

It  has  been  shown,  both  theoretically  and  experimentally,  that  backscattering 
from  a  monopole  can  be  considerably  diminished  and  almost  suppressed  if  the 
monopole  is  grounded  in  a  cavity  and  its  protruding  length  £  does  not  exceed 
0.  45  A.  ;  with  increasing  length,  cavity' -loading  becomes  increasingly  ineffective. 
Suppression  is  optimum  when  the  cavity  depth  is  about  0.  25  X  ;  deeper  recesses 
narrow  the  suppression  bandwidth  considerably. 
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ADJUSTABLE  CAVITY  BAR 


Figure  14.  Cross  Section  of  the  Mechanical  Setup 


Figure  15.  Scheme  of  the  Experimental  Arrangement 


Appendix  A 


1.  EXPERIMENTS 

The  backscattering  from  a  monopole  grounded  in  a  cavity  was  measured  for 
various  heights  of  the  monopole  above  the  ground  plane  and  for  various  depths  of 
the  cavity.  A  variable -depth  cavity  was  obtained  by  machining  a  hole  into  the 
ground  plane  and  precision-fitting  a  sliding  circular  brass  bar  into  the  hole.  A 
brass  rod  was  precision-fitted  into  a  hole  drilled  through  the  center  of  the  bar. 
This  was  the  monopole,  also  variable  in  height  and  independent  of  the  adjustments 
made  in  the  cavity  (Figure  14). 

In  the  block  diagram  (Figure  15)  showing  the  experimental  arrangement  of 
the  apparatus,  the  symbols  are  as  follows: 

T.  a  stabilized  frequency  source  operating  at  9080  Mcps,  modulated  at 
1000  cps 

PA  precision  waveguide  attenuator 
E  transmitter  arm 

HY  hybrid  junction  isolating  arms  H  and  E 

H  receiver  arm,  containing  a  crystal  detector  and  a  bolometer  amplifier 
for  recording  the  signal  level 

Z™  adjustable  impedance  transformer  attached  to  the  third  arm  of  HY  , 
offering  better  isolation  from  the  rest  of  the  circuit  when  used  with  L 
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L  precision  adjustable  load,  attached  to  the  third  arm  of  HY 
TR  fourth  arm,  carrying  the  transmitted  signal  and  reflected  signal 

2.  MEASUREMENT  PROCEDURES 

2. 1.  Without  a  scatterer  on  the  ground  plane,  Zj,  was  tuned  for  the  minimum 
signal  return,  which  was  found  to  be  at  about  85  db  below  the  transmitted  signal, 
--i^ce  this  was  approximately  the  same  as  the  inherent  noise  level  of  the  amplifier, 
the  reference  level  was  arbitrarily  chosen  to  be  10  db  above  noise.  For  all  meas¬ 
urements  the  reflected  signal  was  attenuated  to  the  reference  point,  and  the  amount 
of  attenuation  read  on  a  calibrated  db  scale  was  recorded  at  the  value  of  the  sig¬ 
nal  return. 

2.  2.  For  the  first  set  of  measurements,  the  top  of  the  circular  brass  bar  was 
left  on  a  level  with  the  ground  plane  so  that  no  cavity  existed.  The  top  of  the 
1/ IS -in. -diameter  monopole  was  raised  to  a  height  of  IX  above  'he  ground  plane, 
and  a  measurement  taken.  The  monopole  was  then  successively  depressed  so  that 
its  height  above  the  ground  plane  was  regularly  reduced  in  increments  of  0.  025  in. , 
and  reaomgs  taken  until  the  monopole  had  reached  zero  height  above  the  plane. 

2.  3.  The  monoDcle  was  then  again  extended  to  a  height  of  1 A  above  the  plane. 

This  time  the  circular  bar  was  depressed  into  the  ground  plane,  creating  a  cavity, 
and  the  measurements  described  in  Step  2.  2  were  repeated.  The  depth  of  the 
cavity  was  regularly  increased  in  increments  of  0.  05  inch.  In  each  case  the  mono¬ 
pole  was  raised  to  the  same  IX  height  above  the  ground  plane  and  incrementally 
reduced  until  the  height  was  again  zero. 

2.  4.  To  determine  the  influence  of  the  ratio  of  monopole  diameter  to  wavelength 
on  the  curves,  identical  measurements  v/ere  made  with  a  monopole  of  half  the 
diameter  (1/32  in. )  and  an  appropriate  cavity  bar.  A  comparison  cf  the  data  for 
both  monopoles  shows  that  the  curves  are  essentially  the  same,  except  that  the 
shape  is  more  pronounced  in  the  case  of  the  thin  rod. 

-.  5.  As  a  standard  reference,  a  hemisphere  0.  437  in.  in  diameter  was  placed  at 
the  same  point  at  which  the  monopole  had  stood.  The  return  from  the  hemisphere 
was  noted  to  be  13.  8  db  above  the  arbitrary  reference  level  for  the  monopolc. 
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V.  Scattering  From  Thick  Reactively  Loaded  Rods 

C.  J.  Sleiten,  P.  Blacksmith,  F.  S.  Holt,  and  B.  B.  Gorr 
Microwave  Physics  Laboratory 
•\ir  Force  Cambridpc  Research  Laboratories 
Bedford,  Massachusetts 


1.  INTRODUCTION 

For  objects  whose  dimensions  are  large  in  terms  of  wavelength,  considerable 
reduction  in  backscattering  cross  section,  o  ,  can  be  obtained  in  certain  specified 
directions  by  shaping  the  object.  Again,  for  objects  lax  ge  compared  to  wavelength, 
considerable  reduction  in  cr  for  all  target  attitudes  can  be  attained  by  coating  the 
object  with  absorbing  material.  For  objects  whose  dimensions  are  of  the  order  of 
a  wavelength,  it  was  conjectured  that  the  techniques  of  shaping  and  coating  with 
absorbing  material  would  lose  their  effectiveness  in  reducing  c  ,  and  that  addi¬ 
tional  techniques  would  be  required  for  effective  cont:  o'.  One  such  additional 
technique,  namely  passive  reactive  loading,  has  been  the  subject  of  considerable 
experimental  investigation  at  AFCRL. 

This  paper  is  principally  concerned  with  an  examination  of  the  effectiveness  of 
the  three  techniques  -  shaping,  coating  with  absorbing  material,  and  passive  re¬ 
active  loading  -  individually  and  in  combination,  in  reducing  o  for  objects  of 
resonant  dimensions. 
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2.  EXPERIMENTAL  SETUP  AND  DEFINITIONS 

AH  experimental  measurements  were  made  at  S  band  on  the  AFCRL,  freespace 
scattering  range  at  the  Ipswich  Field  Test  Site,  Ipswich,  Massachusetts.  All  tar¬ 
gets  were  placed  on  a  styrofoam  column  and  rotated  360°  about  the  vertical  axis  of 
the  column.  Elongated  objects  were  always  mounted  with  their  long  dimension 
horizontal;  for  such  objects,  horizontal  polarization  means  polarization  parallel 
to  the  long  dimension,  and  vertical  polarization  means  polarization  perpendicular 
to  the  long  dimension. 

The  terms  peak  and  peak  Oy  are  used  to  designate  the  maximum  a  un¬ 
der,  respectively,  horizontal  and  vertical  polarization  conditions  for  those  target 
aspects  attained  by  rotating  the  objects  360s  about  a  vertical  axis.  The  term 
overall  peak  a  denotes  the  larger  of  peak  cr^  and  peak  Cy.  For  some  of  the 
plots  the  peak  a's  are  normalized  for  each  target  to  the  overall  peak  a  for  that 
target  under  unloaded  er  shorted  conditions.  The  angle  of  rotation  of  the  target 
about  the  vertical  axis  is  designated  by  6  .  Unless  otherwise  explicitly  stated, 
all  elongated  objects  are  approximately  2  in.  long,  that  is,  A/ 2  at  S  band. 


3.  BACKSCATTER  REDUCTION  BY  MEANS  OF  ABSORBING  MATERIAL 

Conductron,  Inc.  absorbing  material.  Sample  No.  17  -  109,  has  been  tested 
on  several  shapes  in  the  resonance  region.  In  all  cases  the  targets  were  metal 
plates  coated  with  the  absorbing  material  on  one  side  only.  The  target  configura¬ 
tions  were  a  1.2  A  by  1.2  A  square,  a  A/2  by  A/2  square,  a  A/2  by  A/4  rec¬ 
tangle,  and  a  A/2  by  A/8  rectangle.  Figure  1  presents  the  results  for  vertical 
polarization,  and  Figure  2  the  results  for  horizontal  polarization. 

The  solid  curves  are  o  vs.  target  angle  0  for  the  metal  side  of  the  target, 
and  the  dashed  curves  are  a  vs.  0  for  the  absorber  side’  of  the  target.  The  a 
scale  is  in  db  with  an  arrow  denoting  a  reference  value  of  -  29.  2  sq  cm  (the 
backscattering  cross  section  of  a  1.  25-in.  -diameter  metal  sphere).  * 

For  vertical  polarization  (Figure  1),  the  absorbing  material  produces  a  large 
reduction  in  peak  a „  for  all  targets.  For  horizontal  polarization  (Figure  2),  the 
absorbing  material  is  very  effective  for  the  square  targets  but  relatively  ineffective 
for  the  elongated  targets  (Figures  2c  and  2d).  Note  that  for  the  metal  side  of  the 
A/2  by  A/8  target,  the  peak  <Ty  is  down  from  the  peak  cr^  about  8  db.  Thus  for 
elongated  targets  of  this  type,  the  condition  under  which  least  reduction  in  ^eak 
can  be  achieved  by  means  of  absorbing  material  (that  is,  horizontal  polaxization) 
is  just  the  condition  under  which  peak  a  is  greatest.  Conversely,  for  the  same 

* 

This  reference  value  ao  appears  in  many  of  the  subsequent  plots. 


type  of  target  the  condition  under  which  greatest  reduction  in  peak  c  can  be 
achieved  by  means  of  absorbing  material  (that  is,  vertical  polarization)  is  just 
the  condition  under  which  the  peak  a  is  least. 


Figure  1.  Backscattering  Cross  Section  a  vs  Target  Angle  0  for  Conductron,  Inc. 
Absorbing  Material  Sample  Number  17-109.  Polarization  Vertical. 
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Figure  2.  Backscattering  Cross  Section  cr  vs  Target  Angle  0  for  Conductron,  Inc. 
Absorbing  Material  Sample  Number  17-109.  Polarization  Horizontal. 

4.  flAC*C*.ATTEB  REDUCTION  BY  MEANS  OF  TARGET  SHAPING 


In  Figure  3  backscattering  cross  section  cr  vs.  target  angle  0  is  shown  for 
a  cylindrical  rod  and  a  truncated  cone  under  horizontal  polarization,  l'he  dimen¬ 
sions  of  the  targets  are  shown  in  the  figure.  Changing  from  the  cylindrical  shape 
to  the  conical  shape,  although  reducing  the  volume  by  a  factor  of  approximately  3, 
did  not  reduce  the  peak  cr^;  rather,  it  slightly  increased  it.  Also  note  that  cr  did 
not  decrease  for  the  nose -on  aspect  of  the  cone  compared  with  the  rod. 


Figure  3.  Backsoattering  Cross  Section  a  vs  Target  Angle  0  for  a  Rod  and  a  Cone.  Polarization 
Horizontal.  Frequency  r  3100.3  Mops. 
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'rhe  effect  of  rod  diameter  on  peak  cry  and  peak  Ojj  for  rods  approximately 
X/2  tong  is  shown  in  Figure  4.  The  curves  show  that  peak  is  essentially  inde¬ 
pendent  of  rod  diameter  over  the  range  of  values  considered  and  that  peak  &y 
increases  with  rod  diameter.  For  rods  approximately  X/2  long,  it  is  clear  that 
overall  peak  cr  occurs  for  horizontal  polarization.  It  is  also  clear  that  peak  Oy 
can  be  controlled  by  rod  diameter  but  peak  cannot. 

If  the  desired  level  of  reduction  in  overall  peak  a  is  arbitrarily  set  at  15  db 
below  the  overall  peak  cr  for  a  solid  X/2  rod,  then,  at  least  for  vertical  polariza¬ 
tion,  this  can  be  achieved  by  using  a  rod  whose  diameter  is  X/8  or  less.  Clearly 
then,  to  achieve  the  desired  level  of  reduction  in  overall  peak  cr  for  a  rod,  peak 
Ojj  must  be  reduced.  Neither  the  absorbing  material  nor  the  target  shapes  tested 
were  effective  in  reducing  pea„.  cr^.  This  leads  then  to  consideration  of  the  pas¬ 
sive  reactive  loading  technique. 

3.  BACKSCATTER  REDUCTION  BY  MEANS  OF  PASSIVE  REACTIVE  LOADING 

Early  work  in  the  technique  of  ’•educing  peak  ov,  appears  in  a  patent  filed  in 
1  H 

1946  by  Harley  A.  Iarr.s.  In  this  disclosure  lams  described  the  use  of  coaxial 

loading  together  with  dielectric  coating  to  effect  wide-band  scatter  reduction  for 

metal  spacer  posts  in  parallel  plate  construction.  lams  claimed  that  symmetrical 

loading  at  separated,  symmetrically  placed  points  (double  loading)  produced 

wider-band  performance  than  asymmetrical  loading  at  separated  asymmetrically 

placed  points  or  loading  at  one  point  (single  loading). 

In  1948  C.  J.  Sletten  at  AFCRL  successfully  designed  spacer  posts  for  the 

2 

parallel  plate  region  of  the  Volir  Antenna  ,  using  the  lams  idea  Sletten  also 
attempted  to  multiply  load  a  rod  to  produce  a  low  backscattering  freespace  object. 
This  viiempt  was  unsuccessful  but  Sletten1  s  ideas  motivated  both  a  theoretical  and 
experimental  investigation  of  passive  reactive  loading  of  monopoles  over  a  ground 

3 

plane  to  reduce  backscattering  by  Gerbes,  Poehler,  and  Kearns  in  1957. 

In  1958  AFCRL.  supported  further  research  in  passive  reactive  loading  of  thin 

monopoles  (about  X/45  in  diameter)  over  a  ground  plane  to  reduce  backscattering  by 
°  4 

As  and  Schmitt  at  Cruft  Laboratory,  Harvard  University.  The  Harvard  experimental 
results  agreed  well  with  their  theoretical  results  and  included  multiple  as  well  as 
single  loading. 

Independent  work  using  passive  reactive  loading  techniques  has  been  carried 
out  by  Andrew  Alford  Consulting  Engineers,  Inc.  in  the  reduction  of  scatter  from 
radio  towers.  Currently  this  company  is  investigating  the  possibilities  of  wide¬ 
band  backscatter  reduction  by  means  of  active  control  of  reactive  loading. 

In  the  last  few  years  under  Air  Force  support  the  Ohio  State  University 
Research  Foundation  and  the  University  of  Michigan  Radiation  Laboratory  have 
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Figure  4.  Peak  a  vs  Rod  Diameter  for  Rod  Length  \/2 

investigated  reactive  loading  of  various  configurations,  including  rods  and  slots, 
to  reduce  backscattering. 

Current  investigations  at  AFCRL  have  been  primarily  concerned  with  the  effec¬ 
tiveness  of  passive  reactive  loading  in  reducing  peak  <rH  for  metal  rods  approxi¬ 
mately  X/2  long  and  X/8  to  3X/8  in  diameter.  Three  typical  single -loaded  rods 
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are  shown  in  Figure  5,  and  cross  sections  with  dimensions  are  shown  in  Figure  6. 
Note  the  gap  extending  circumferentially  around  the  rod  and  the  symmetrical 
coaxial  loading  cavities.  Figure  7  shows  a  typical  double -loaded  rod. 


Figure  5.  Typical  Single  Loaded  Target  Models 

Curves  of  a  vs.  target  angle  Q  for  the  target  shown  in  Figure  6a  under 
horizontal  polarization  conditions  as  a  function  of  the  total  loading  cavity  length 
C  are  shown  in  Figure  8.  The  curve  labeled  C  =  0  is  for  the  unloaded  or  solid 
rod.  As  C  increases,  a  drops  significantly  until,  at  C  =  0.281  in.,  peak  crH 
has  dropped  approximately  25  db  below  overall  peak  <j.  Note  that  in  the  broadside 
aspect  cr  has  been  reduced  by  approximately  40  db.  The  value  of  C  for  minimum 
peak  Ojj,  designated  as  Cmin,  is  most  easily  determined  from  a  plot  of  normalized 
peak  Opj  vs.  C.  Curves  of  this  type  are  shown  in  Figure  9  for  two  different  rod 
diameters  (see  Figures  6a  and  6b).  For  D  -  0..5  in.,  Cm-n  =  0.281  in.;  for 
D  ~  1.0  in.,  C  .n  =  0.105  in.  Apparently  the  tuning  for  minimum  peak  is  less 
critical  for  D  =  1.0  in.  than  for  D  =  0.5  in.,  as  evidenced  by  the  broader  minimum. 

The  dual  cavities  used  in  the  loading  configurations  (see  Figure  10)  are  effec¬ 
tively  in  series.  This  was  experimentally  verified  by  using  only  one  cavity  and 
tuning  for  minimum  peak  c^.  The  resultant  inductive  load  was  equal  to  the  sum  of 
the  loads  of  the  two  symmetric  cavities  of  the  usual  configuration  under  minimum 
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peak  Ojj  conditions.  The  circumferential  gap  that  is  present  in  the  loading  con¬ 
figuration  acts  as  a  shunt  capacitance  across  the  cavity  loads.  To  determine  the 
effect  of  the  gap  width  G  on  the  load  impedance  i  squired  for  minimum  peak  <7^, 
loaded  rods  wit!  gaps  of  0.040  in.  and  0.030  in.  were  tested.  Figure  10  shows 
that  for  the  two  different  gap  sizes  the  individual  cavity  depths  C*  required  for 
minimum  peak  conditions  were  nearly  equal.  It  is  therefore  concluded  that 
the  effect  of  the  shunt  capacitance  associated  with  the  gap  is  small. 


Figure  7.  Typical  Double  Loaded  Target  Model 


The  effect  of  rod  diameter  D  and  rod  length  L  on  C  is  shown  in  Figure  11. 

From  Figure  11a  it  is  clear  that  C  .  decreases  as  D  increases,  and  from 

.  min 

Figure  lib  it  is  apparent  that  C  ^  changes  little  with  L  over  the  range  con¬ 
sidered.  The  latter  indicates  that  there  is  some  hops  for  bandwidth. 


6.  M0N0STATIC,  BISTATIC,  AND  BANDWIDTH  MEASUREMENTS  ON  V  ARIOUS  TARGETS 


6.1  Monostatic  Comparison  of  Cylindrical  Shape  »  ith  Double  Conical  Shape 

A  comparison  of  the  backscattering  from  a  cylindrical  rod  loaded  and  unloaded 
with  the  backscattering  from  a  double  conical  shape  loaded  and  unloaded  is  shown 


Figure  8.  Hackscattering  Cross  Section  a  vs  0  Target  Angle  for  Various  Loadings  of 
Loaded  Hod.  Polarization  Horizontal.  Frequency  =  3100.8  Mcps. 
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Figure  11.  Total  Cavity  Length  for  Minimum  Backscattering  Cmin  vs  Rod 
Diameter  D  and  Rod  Length  L.  Polarization  Horizontal.  Frequency  =  3100.  8 
Mcps. 
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in  Figure  12.  Dimensions  of  the  objects  are  given  in  Figures  6a  and  13a.  The 
solid  curves  of  Figure  12  are  cr  vs.  9  for  the  unloaded  configurations,  and  the 
dashed  curves  are  cr  vs.  9  for  the  loaded  configurations  tuned  for  minimum  peak 
ajj-  Comparison  of  the  solid  curves  of  Figures  12a  ana  12b  shows  that  peak  cr^, 
the  maximum  a  under  horizontal  polarization  condition' .  is  the  same  for  both 
target  shapes.  The  dashed  curves  in  Figures  12a  and  12b  show  the  large  reduc¬ 
tion  in  peak  to  be  gained  for  both  shapes  by  reactive  loading. 

For  both  shapes  with  polarization  vertical,  it  is  strikingly  evident  from  com¬ 
parison  of  the  solid  and  dashed  curves  in  Figures  12c  and  12d  that  the  reactive 
loading  has  practically  no  effect  on  v,  that  is,  there  is  very  little  difference 
between  the  solid  and  dashed  curves. 

The  conclusion  here  is  that  for  vertical  polarization  the  loading  cavities  as 
presently  designed  do  not  couple  with  the  induced  currents  on  the  surfaces  of  the 
targets.  Comparison  of  Figures  12c  and  12d  indicates  that  peak  Oy  is  about  6db 
less  for  the  double  conical  shape  than  for  the  cylindrical  rod  shape.  This  result 
is  to  be  expected  since  the  average  diameter  of  the  double  cone  is  considerably 
lees  than  that  of  the  rod. 

Clearly,  then,  changing  the  shape  from  cylindrical  to  double  conical  did  not 
reduce  the  overall  peak  cr,  whereas  reactive  loading  reduced  the  overall  peak  cr 
by  at  least  15  db.  The  condition  for  which  shaping  had  an  appreciable  effect,  that 
is,  vertical  polarization,  is  just  the  condition  for  which  peak  Oy  for  the  rod  is 
already  down  16db  from  peak  o^.  Note  the  similarities  of  these  conclusions  to 
those  obtained  from  using  absorbing  material. 

6.2  Bandwidth  Characteristics 

Three  differert  reactively  loaded  target  configurations  were  examined  over 
the  frequency  range  2500  to  3000  Mcps.  In  each  case  the  target  was  tuned  for 
minimum  peak  at  a  frequency  near  the  center  of  the  band.  The  configurations 
were  a  single  loaded  rod,  a  single  loaded  double  conical  shape,  and  a  double  loaded 
rod.  The  dimensions  of  the  targets  are  shown  in  Figures  6b,  13a,  and  13b.  Curves 
of  normalized  peak  cr^  vs.  frequency  are  shown  in  Figure  14,  and  the  approximate 
-15  db  bandwidth  is  indicated  for  each  target. 

Apparently  double  loading  produces  considerably  wider  bandwidth  than  single 
loading,  while  double  conical  shaping  with  a  single  loading  results  in  even  wider 
bandwidth.  It  is  interesting  to  note  here  that  although  shaping  does  not  decrease 
the  overall  peak  a  at  a  single  frequency,  it  does  apparently  increase  the  bandwidth. 

6.3  Bistatic  Characteristics 

Bistatic  reflection  measurements  for  horizontal  polarization  were  made  on  a 
single  loaded  rod  (see  Figure  13c)  and  on  a  double  loaded  rod  (see  Figure  13b). 
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Figure  12.  Backscattering  Cross  Section  a  vs  Target  Angle  Q  for  Single  Loaded 
Rod  (Figure  6a)  and  Single  Loaded  Double  Conical  Shape  (Figure  13a). 
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Figure  13.  Dimensions  of  Target  Models 


For  each  bistatic  angle  a  ,  each  target  was  rotated  260’  and  the  peak  bistatic  cross 
section  determined.  These  measured  peak  bistatic  cross  sections  were  then  nor¬ 
malized  to  the  overall  peak  backscattering  cross  section  oi  the  unloaded  configura¬ 
tion.  Curves  of  normalized  peak  bistatic  cross  section  ”S.  bistatic  angle  a  for 
both  targets  are  shown  in  Figure  15.  For  each  target  the  peal:  bistatic  cross 
section  increases  with  btstatie  angle;  in  the  neighborhood  of  a  -  90’  it  approaches 
the  reference  value,  that  is,  the  overall  peak  a  for  an  unloaded  or  solid  rod.  For 
bistatic  angles  in  the  range  a  =  0®  to  a=  30®,  the  single  loaded  ,-od  has  lower 
peak  return  than  the  double  loaded  rod  though  both  are  down  at  least  18  db  from 
the  reference  value.  The  peak  return  from  the  double  loaded  rod  is  15  db  or  more 
down  from  the  reference  value  o/er  the  bistatic  angular  range  a  =  0®  to  a  =  0.48®. 

7.  CONCLUSIONS 

For  elongated  objects  of  length  approximately  A/2,  coating  them  with  absorb¬ 
ing  material  or  shaping  them  into  the  double  conical  form  does  not  significantly 
reduce  overall  peak  a.  The  condition  under  which  these  two  techniques  have  their 
greatest  effect  in  reducing  a,  that  is,  vertical  polarization,  is  just,  the  condition 
under  which  the  a  for  an  elongated  object  is  already  well  below  its  overall  peak 
value.  Fortunately,  the  condition  under  which  overall  peak  <j  occurs  and  under 
which  the  two  above  techniques  are  least  effective  in  backscatter  reduction,  that 
is,  horizontal  polarization,  is  just  the  condition  under  which  reactive  loading  has 
its  greatest  effect.  Therefore,  the  three  techniques  tend  to  complement  each  other, 
and  it  is  reasonable  to  expect  that  a  combination  of  all  three  techniques  should  pro¬ 
duce  a  good  wideband  minimum  backscattering  object. 


8.  CURRENT  AND  FUTURE  R0RK 

The  following  investigations  are  either  currently  in  progress  or  will  be  under¬ 
taker  in  the  near  future: 

a.  Use  of  slightly  lossy  rather  than  purely  reactive  loading  to  improve  band¬ 
width  characteristics. 

b.  Redesign  of  the  slot  configuration  so  Uiat  it  will  couple  into  the  surface 
currents  under  vertical  polarization  conditions. 

c.  Determination  of  the  existence  of  higher  harmonic  resonances  of  the  load¬ 
ing  cavity  that  may  undesirably  enhance  the  scattering. 

d.  Design  of  loading  configurations  to  successfully  reduce  the  scattering 
from  spheres  and  jacks  in  the  resonance  region. 
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V!.  Analysis  of  Loaded  Terminal  Scatterers* 

Eduard  M.  Kennaugh 
Antenna  Laboratory 
Department  of  Electrical  Engineering 
The  Ohio  Slate  University 
Columbus.  Ohio 


I.  INTRODUCTION 

The  echoing  properties  of  antennas,  or  scatterers  with  one  or  more  terminal 
pairs,  have  been  the  subject  of  continuing  studies  at  our  Laboratory.  In  addition 
to  the  new  scattering  parameters  introduced  by  arbitrary  terminal  loadings,  we 
are  interested  in  the  subtle  difference  between  scattering  by  "  good"  antennas  and 
by  arbitrary  objects.  Since  a  "good"  antenna  acts  as  an  efficient  device  for  trans¬ 
fer  of  energy  from  a  wave  field  to  a  terminal  pair  when  receiving,  and  in  the  re¬ 
verse  direction  when  transmitting,  it  appears  that  its  scattering  properties  should 
possess  certain  unique  features.  Among  t!  ese,  of  course,  is  the  strong  depend¬ 
ence  of  these  properties  upon  the  coupling  between  the  incident  wave  and  the 
antenna  upon  termination. 

Several  of  the  applications  of  this  study  can  be  described.  First,  the  deter¬ 
mination  of  the  parameters  that  relate  scattering  by  an  antenna  to  its  load,  and 
the  selection  of  loads  that  will  maximize  or  minimize  the  echo  area.  Next,  the 


*The  research  reported  has  been  supported  in  part  by  Aeronautical  Systems 
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design  of  more  efficient  modulated  scatterers  for  use  m  field  probing,  passive 
communication  links,  or  in  radar  target  simulation.  In  many  cases,  the  control 
of  the  echo  from  a  specified  object  is  desired,  although  neither  an  antenna  nor 
terminal  pairs  have  been  specified.  Use  of  antenna  concepts  can  still  be  made, 
however,  to  find  the  characteristic  modes  that  can  most  significantly  contribute 
to  echo  area  and  determine  what  loadings  are  required.  Finally,  a  unified  theo¬ 
retical  treatment  of  scattering  and  radiating  properties  may  best  be  achieved,  at 
least  for  objects  of  resonant  size,  by  considering  these  objects  to  be  multi-mode 
antennas. 


2.  EFFECT  OF  TERMINATION  UPON  ANTENNA  ECHO  AREA 


Initially,  a  fixed-source  frequency  and  aspect  were  assumed,,  and  the  scat¬ 
tering  properties  of  antennas  or  bodies  with  antennas  mounted  upon  them  were 
investigated  as  functions  of  the  terminating  impedance.  It  is  easily  shown  that 
the  echo  signal  is  a  function  of  three  complex  parameters  and  the  load  impedance: 


Es  ZT  +ES  Z„ 
-j-,s  oc  L  sc  A 

E  =  — Z7+YZ - 


(1) 


g 

where  E  is  the  phasor  echo  signal  received  when  the  antenna  is  terminated  with 

s  o 

ZT ,  E  is  that  received  when  terminated  with  an  open  circuit,  E  is  that  re- 

L»  OC  SC 

ceived  with  a  short  circuit  termination,  and  Z^  is  the  antenna  impedance  at  the 
load  terminals.  A  more  useful  form  of  this  equation  is  readily  derived  by  use  of 

g 

the  phasor  E..,  the  phasor  echo  signal  received  when  the  antenna  is  terminated 
* 

with  Z^ ,  the  complex  conjugate  of  the  antenna  impedance: 


ES  =  EM-r(EM-Eoc>- 


(2) 


is  the  modified  voltage  reflection  coefficient  corresponding  to  the  load  Z^j 


T  = 


ZL  -  ZA 
ZL  +_iZ,A 


(3) 


Note  that  the  definition  of  r  differs  from  that  normally  used,  whenever  Z .  is 
not  a  real  quantity,  and  is  such  that  r  will  always  lie  on  the  unit  circle  for  re¬ 
active  loads  Zr  and  inside  the  unit  circle  for  passive  lossy  loads. 

A  graphical  interpretation  of  Eq.  (2)  is  that  the  phasor  signal  received  as  a 
function  of  load  impedance  ZT  is  proportional  to  the  difference  between  a  fixed 

*  Li 
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phasor  and  the  phasor r-  Thus,  if  one  plots  the  modified  reflection  coefficient 
on  a  pseudo-Smith  chart,  the  vector  from  r  to  the  tip  of  the  fixed  phasor  has  a 
length  proportional  to  the  echo  signal  received  as  r  (or  Zj)  is  varied.  Such  a 
geometrical  interpretation  is  shown  in  Figure  1  for  several  values  of  the  fixed 
phasor  component. 


Figure  1.  Dependence  of  Phaser  Echo  Signal  on  Load  Reflection  Coefficient 

We  might  assume  that  three  different  aspects  of  a  single  antenna  correspond 
to  the  phasor  values  OP^  ,  OPg,  OP^.  For  aspect  1,  the  fixed  component  exceeds 
the  load -dependent  component  and  it  is  not  possible  to  eliminate  the  echo  signal  by 
use  of  a  passive  load,  but  it  can  be  varied  between  amplitudes  proportional  to 
OP,  +  1  and  OP,  -  1  by  use  of  reactive  loads.  For  aspect  2,  the  fixed  component 
equals  the  load-dependent  component  and  the  echo  signal  can  be  eliminated  or 
maximized  by  use  of  reactive  loading.  For  aspect  3,  the  fixed  component  is  less 
than  the  load-dependent  component  and  the  echo  signal  can  be  eliminated  by  using 
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a  lossy  load  corresponding  to  the  point  P  ^  or  maximized  by  a  reactive  load. 

One  must  bear  in  mind  that  the  geometrical  diagram  omits  a  factor  of  pro¬ 
portionality,  so  that  only  relative  echo  signal  amplitude  and  phase  can  be  deter¬ 
mined;  the  factor  omitted  is  proportional  to  the  geometric  mean  of  the  antenna 
power  gains  in  the  directions  of  source  and  receiver  as  well  as  the  usual  radar 
range  parameters. 

Since  the  variation  of  echo  signal  with  load  is  uniquely  related  to  antenna  and 
scatterer  parameters,  we  have  considered  the  use  of  variable  reactive  loads  to 
determine  these  parameters  solely  through  scattering  measurements.  Combining 
graphical  and  analytical  techniques,  it  is  possible  to  quickly  reduce  echo  signal 
amplitude  data  for  several  calibrated  reactive  loads  to  determine  antenna  imped¬ 
ance,  gain,  and  scattering  parameters  at  a  single  aspect  and  frequency.  These 

12  3 

techniques  and  their  application  are  fully  described  by  Garbacz.  *  ' 


3.  AVERAGE  ECHO  AREA  AND  TOTAL  SCATTERING  CROSS  SECTION 

The  graphical  interpretations  of  Figure  1  are  limited  to  individual  aspects 
of  the  scattering  antenna,  but  it  is  frequently  necessary  to  consider  how  average 
echo  area  over  a  range  of  aspects  or  polarizations  varies  with  antenna  termina¬ 
tion.  Green*  has  shown  how  a  slight  modification  to  Figure  1  can  be  made  to 
permit  its  use  in  such  a  case.  We  merely  add  a  positive  constant  to  the  square  of 
the  distance  previously  derived  for  echo  signal  amplitude  in  Figure  1  to  obtain  a 
quantity  proportional  to  the  average  echo  area  and  to  determine  its  dependence 
upon  T.  As  shown  in  Figure  2,  the  average  echo  area  of  an  antenna  is  propor¬ 
tional  to  a  constant  C  plus  the  square  of  the  distance  PL,  where  L  is  the  point 
corresponding  to  the  modified  reflection  coefficient  of  the  load  defined  by  Eq.  (3). 

Depending  upon  the  position  of  the  fixed  point  P,  one  would  achieve  minimum 
average  echo  area  with  a  reactive  (P  outside  or  on  unit  circle)  or  lossy  load  (P 
inside  unit  circle).  In  every  case,  maximum  average  echo  area  can  be  obtained 
with  a  reactive  load.  The  echo  power  can  be  averaged  in  this  manner  over  orien¬ 
tations  of  source  and  of  receiver,  including  variations  in  polarization,  but  it  is 
assumed  that  the  source  frequency  remains  constant.  From  averaged  cross  sec¬ 
tion  values  obtained  with  as  few  as  four  loads,  it  is  poss  Ae  to  construct  a  dia¬ 
gram  such  as  Figure  2  that  will  yield  the  average  cross  section  for  any  load 
impedance. 

Of  special  interest  is  the  case  where  the  source  remains  fixed  and  the  receiv¬ 
ing  antenna  varies  over  all  bistatic  angles  to  obtam  the  average  scattering  cross 
section  of  an  antenna.  One  may  then  consider  how  the  total  energy  scattered  by 
an  antenna  compares  with  that  absorbed  as  a  function  of  termination.  In  the  early 
literature  on  antenna  scattering,  it  was  often  erroneously  assumed  that  an  antenna 
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must  back-scatter  at  least  as  much  energy  as  it  absorbs.  Although  no  such  re¬ 
striction  exists,  a  mere  fundamental  relation  possibly  applies  between  total  scat¬ 
tering  and  absorption  cross  sections  of  an  antenna.  If  we  postulate  that  an  antenna 
may  not  absorb  more  than  it  scatters,  this  would  imply  that  the  value  of  C  and  the 
location  of  P  in  Figure  2  would  be  restricted  in  the  case  of  averaged  bi static 
scattering. 


Figure  2.  Dependence  of  Average  Echo  Power  on  Loa'’  Reflection  Coefficient 

Green ^  used  the  concept  of  Figure  2  to  present  calculated  average  scattering 
cross  section  for  a  cylindr  xal  antenna  as  a  function  of  load,  and  included  contours 
of  constant  average  scattering  and  constant  absorption  cross  section  on  the  chart, 
as  shown  in  Figure  3. 

It  is  possible  to  fit’ j  a  load  impedance  that  will  maximize  the  ratio  of  absorp¬ 
tion  to  scattering  cross  section;  the  maximum  ratio  obtained  in  this  case  is  4 ■  ", 
obtained  at  point  B  where  the  power  absorbed  is  approximately  50  percent  of  that 
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for  a  matched  termination.  Note  that  the  scattering  and  absoi  ption  cross  sections 
are  approximately  equal  when  the  absorption  coefficient  is  a  maximum. 


isr ;t SCATTERING  CROSS-SECTION 
o-A  =  ABSORPTION  CROSS-SECTION 


Figure  3.  Contours  of  Constant  Total  Scattering  and  Absorption  Cross 
Sections  as  Function  of  Load  Reflection  Coefficient,  Cylindrical  Antenna 

The  relation  between  scattering  and  absorption  cross  sections  of  any  object 
is  implied  by  its  albedo,  which  is  the  ratio  of  scattering  cross  section  to  extinc¬ 
tion  (total)  cross  section.  Since  the  extinction  cross  section  is  the  sum  of  scat¬ 
tering  and  absorption  cross  sections,  an  object  with  an  albedo  less  than  1/2 
absorbs  more  than  it  scatters. 

5 

Results  obtained  by  Deirmendjian  for  the  scattering  by  lossy  dielectric 
spheres  show  that  the  albedo  is  below  1/2  for  lossy  dielectric  spheres  less  than 
2  wavelengths  in  diameter,  for  a  complex  index  of  refraction  1.  315  -  i.  4298.  For 
the  same  index,  the  ratio  of  absorption  to  scattering  cross  sections  is  approxi¬ 
mately  5  for  a  sphere  circumference  of  1  wavelength,  where  the  absorption  cross 
section  is  slightly  less  than  the  geometrical  cross  section. 

Tc  what  extent  can  a  lossy  dielectric  sphere  be  considered  an  antenna  with  a 
lossy  termination?  From  a  practical  point  of  view  there  does  not  appear  to  be  any 
way  to  extract  useful  power  from  such  a  structure,  but  a  more  fundamental  differ¬ 
ence  is  the  distributed  loss  rather  than  the  lumped  loss  corresponding  to  a  definite 
tertninal  pair. 
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4.  APPLICATION  OF  ANTENNA  CONCEPTS  TO  ARBITRARY  LOADINGS 

It  would  be  useful  to  fit  more  general  types  of  leading  or  scatterer  perturba¬ 
tion  to  the  model  of  Figures  1  and  2.  For  example,  the  effect  of  a  slot  or  gap  of 
variable  depth  in  a  conducting  object  might  be  considered  as  a  variable  reactive 
load  on  a  hypothetical  terminal  pair.  From  Figure  1,  if  this  model  may  be  used, 
we  would  expect  the  phasor  echo  signal  to  trace  out  a  circular  locus  as  the  depth 
is  varied.  Through  experimental  or  theoretical  studies  of  the  effect  of  gap  depth, 
we  should  then  be  able  to  determine  the  applicability  of  Figure  1  to  more  general 
types  of  loss-less  perturbations.  Similarly,  it  may  be  shown  that  a  purely  resis¬ 
tive  load  variation  should  lead  to  a  circular  locus  for  the  phasor  scattered  signal, 
and  it  might  be  possible  to  assess  the  effect  of  changes  in  conductivity  of  a  part  of 
the  scatterer  by  the  single  terminal  pair  concept  with  an  effective  lumped  load  of 
variable  resistance. 

In  short,  when  the  single -mode  concept  applies,  we  should  expect  a  circular 
locus  for  the  phasor  echo  signal  under  loss-less  perturbations  of  the  scatterer. 
When  this  is  not  obtained,  it  may  be  possible  to  treat  the  scatterer  as  a  multi¬ 
terminal  antenna  with  variable  loading  at  more  than  one  terminal  pair.  The  analy¬ 
sis  of  such  configurations  is  currently  under  study. 


5.  ANALYSIS  OF  ARBITRARY  SCATTER  AS  MULTI-MODE  ANTENNA 


When  one  considers  an  arbitrary  scatterer  as  a  multi-mode  antenna,  it 
follows  that  the  echo  area  as  a  function  of  source  and  receiver  directions  can  be 
expressed  as 


*<es.  *s;er,  V  -£  C.FjW,.  ^(e,.  #r) 


(4) 


where  es>  <ps  are  the  spherical  coordinates  of  the  source  and  0r>  <f>r  those  of 
the  receiver  (both  assumed  infinitely  remote),  and 

F.(0,  <p )  =  radiation  patterns  of  the  individual  modes  and 
C.  =  constants  determined  by  "terminal"  conditions. 

It  is  assumes  that  the  mode  fields  are  mutually  orthogonal  and  that  there  are  a 
finite  number  with  jcj  greater  than  some  arbitrary  minimum  value.  The  number 
required  wo. '.id  necessarily  increase  with  the  size  of  the  scatterer,  but  in  the 
resonance  region  the  number  would  be  less  than  ID  in  many  cases. 


These  orthogonal  or  "normal"  modes  as  we  have  called  them  are  simply  re¬ 
lated  to  the  eigenvectors  of  the  target  scattering  matrix,  and  the  are  determined 
from  the  corresponding  eigenvalues.  The  associated  patterns  are  useful  for  analy¬ 
sis  of  the  transmitting,  receiving,  or  scattering  properties  of  the  object.  We 
would  expect  the  effect  of  perturbations  to  the  object  shape,  such  as  those  pro¬ 
duced  by  slots,  cavities,  or  changes  in  surface  impedance,  to  be  represented  by 
changes  in  the  constants  CL,  as  reflecting  the  change  in  "  terminal"  conditions. 
Further,  the  relative  importance  of  the  various  modes  is  indicated  by  the  magni¬ 
tude  of  the  associated  C.,  so  that  one  may  determine  what  radiating  or  scattering 
patterns  will  be  most  readily  excited  on  the  object. 

Our  present  goals  are  to  extract,  by  a  combination  of  theoretical  and  experi¬ 
mental  methods,  the  "normal"  mode  patterns  of  a  scattering  object  (at  a  single 
frequency)  and  to  determine  the  associated  constants  Ci  as  influenced  by  various 
perturbations  of  the  scatterer  surface  and  by  lumped  loadings. 


6.  EXPERIMENTAL  INVESTIGATIONS 

Experimental  techniques  have  been  developed  for  measuring  antenna  parame¬ 
ters  with  conventional  echo-measuring  systems,  using  modulated  or  variable  re¬ 
active  loads.  To  permit  rapid  measurement  of  echo-signal  amplitude  as  a  function 
of  load,  a  self-contained  programmed  load  unit  was  devised.  This  battery-driven 
device  can  be  mounted  within  a  model  of  a  scattering  object  containing  an  antenna 
and  will  automatically  vary  the  position  of  a  sliding  coaxial  or  waveguide  short 
through  a  full  cycle,  giving  characteristic  pips  on  the  echo  pattern  corresponding 
to  each  1/72  in.  of  travel.  Typical  10  kMc  echo  patterns  obtained  from  a  stub 
antenna  on  a  7.  0  -  in.  -  diameter  sphere  and  from  c  1  1/2  in.  x  2  7/8  in.  compound 
rectangular  horn  are  shown  in  Figure  4. 

A  second  means  of  varying  antenna  termination  under  consideration  employs 
an  electrically  switched  load.  As  a  part  of  this  study,  the  echo  area  variation  with 
bias  of  a  germanium  switching  diode  was  investigated.  Figure  5  presents  the 
variation  in  broadside  echo  area  obtained  with  bias  over  a  range  of  X-band  frequen¬ 
cies.  It  is  clear  that  large  changes  in  cross  section  of  an  object  in  the  resonance 
region  can  be  obtained  through  electrical  tuning. 

The  effect  of  a  rectangular  slot  upon  the  echo  area  of  a  small  flat  plate  was 
also  studied.  An  analysis  by  Green  using  the  superposition  of  plate  scattering 
and  slot  radiation  patterns  to  predict  the  effect  of  the  slot  was  verified.  By  use  of 
a  small  shorting  bar  at  the  center  of  the  slot,  the  echo  area  of  the  slot-plate  com¬ 
bination  can  be  varied  by  approximately  15  db.  Some  of  the  experimental  results 
are  shown  in  Figure  6. 
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Figure  4.  Measured  Antenna  Echo  With  Programmed  Load  Variations 


7.  CONCLUSION 

Various  concepts  of  interest  in  the  analysis  ot  scattering  by  objects  that 
possess  one  or  more  terminal  pairs  (real  or  virtual)  have  been  described.  We 
have  briefly  summarized  the  contribution  of  our  co-workers,  R.  Green,  R.  Garbacz, 
D.  Moffatt.  and  S.  Skarote.  Additional  information  can  be  obtained  from  the 
References  and  Bibliography  listed  below. 
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Figure  6.  Echo  Patterns  of  Rectangular  Plate  With  Slot 
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VII.  Some  Bounds  to  the  Behavior  of  Small  Resonant 

Scatterers 

Roger  F.  Harrington 
Electrical  Engineering  Department 
Syracuse  University 
Syracuse  10,  New  lork 


A  loaded  scatterer  is  an  object  to  which  one  or  more  lumped-parameter  im¬ 
pedance  elements  are  connected.  It  has  been  shown  that  a  small  scatterer  can  be 

resonated  to  greatly  enhance  its  scattering  cross  section.1  In  the  vicinity  of  reso 

2 

nance,  the  echo  area  of  a  small  singly  loaded  scatterer  is  given  by 


u  ~ 


iT 


GRin 


Z.  +  ZT 
in  L 


2 


(1) 


where  \  is  the  wavelength,  G  is  the  directive  gain  of  the  scatterer  when  used  as 

a  transmitting  antenna,  Z.  =  R.  +  jX.  is  the  input  impedance  to  the  scatterer 

in  in  m 

when  used  as  an  antenna,  and  Z^  is  the  loading  impedance. 

According  to  Eq.  vl),  the  echo  area  of  a  small  resonant  scatterer  is  completely 
determined  by  its  properties  when  used  as  an  antenna.  Hence,  known  limitations 
to  the  behavior  of  antennas  can  be  used  to  obtain  corresponding  limitations  to  the 
behavior  of  scatterers.  If  the  general  formula  for  echo  area1  is  used  instead  of 
Eq.  (1),  bounds  to  the  behavior  of  scatterers  of  arbitrary  sizes  and  loads  can  be 
obtained. 


102 
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The  following  are  results  abstracted  from  the  more  general  antenna  theory. 

Let  a  denote  the  radius  of  the  smallest  sphere  that  can  contain  an  antenna.  Then 
for  small  (ka  <  1)  loss-free  antennas,  the  minimum  quality  factor  obtainable  is 

Qn,in  ■  5  I  (k*r‘  +  (ka)'3l  (2) 

where  k  =  2-/\  is  the  wavenumber.  This  minimum  Q  is  obtained  when  the  electric- 
dipole  mode  and  the  magnetic -dipole  mode  are  equally  excited.  If  a  high  Q  antenna 
is  resonated  by  an  input  reactance,  the  bandwidth  of  the  antenna  is  related  to  the 
Q  by 


(3) 


where  o?r  is  the  resonant  frequency  and  Aw  is  the  frequency  difference  between 
points  for  which  |z|  has  fallen  to  its  value  at  resonance.  Under  the  condi¬ 

tion  of  minimum  Q,  the  maximum  antenna  gain  is 


G 


max 


=  3 


(4) 


for  small  loss-free  antennas.  If  only  the  electric -dipole  mode  or  only  the  magnetic- 
dipole  mode  is  excited;,  instead  of  both,  then  the  minimum  Q  is  twice  that  of  Eq.  (1) 
and  the  gain  is  one -half  that  of  Eq.  (4). 

Applying  the  above  results  to  small  resonant  scatterers,  from  Eqs.  (1)  and  (4) 

one  finds  that  the  maximum  echo  area  obtainable  from  a  small  (ka  <  1),  loss -free, 

resonant  (Zr  =  -iX.  )  scatterer  is 
j_i  in 


a  =  -  X2  *  2.86  A.2 
max  tr 


(5) 


In  the  case  of  an  actual  scatterer,  losses  can  reduce  ct, 
tion  (5)  can  oe  compared  to  an  electric -dipole  scatterer, 
scatterer,  for  which  the  resonant  echo  area  is* 


as  discussed  later.  Equa- 
or  a  magnetic -dipole 


adipole 


9  i2  -  i 
4x  "4  ffmax  * 


(6) 
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Hence,  smart,  loss -free,  dipole  scatterers  can  give  one -quarter  the  maximum 
possible  echo  area. 

2 

If  has  been  noted  tha..  because  a  small  resonant  scatterer  behaves  as  a  series 
or  parallel  resonant  circuit  it  is  characterized  by  a  gain -bandwidth  product 

*J~o  P  =  constant.  (7) 


Using  Eqs.  (2),  (3),  and  (5)  in  Eq.  (7),  one  finds  the  maximum  gain-bandwidth 
product  for  small,  loss-free,  resonant  scatterers 


('/a’  p) 


max 


(3X_)  -  2 - 3.38 X(ka)3 

nTF"  (ka)  1  +  (ka) 


(8) 


where  a  is  the  radius  of  the  smallest  sphere  that  can  contain  the  scatterer.  For 
resonant  dipole  scatterers  (electric  or  magnetic),  both  \HT and  P  are  at  most  one- 
half  the  m.ximum.  Hence, 

^  ^dipole  -I  (^}max  (9) 


showing  that  gain-band"idth  products  for  dipoles  cannot  be  greater  than  one- 
quarter  the  optimum. 

When  practical  scatterers  are  made  small  compared  to  wavelength  the  radia 
tion  resistance  becomes  small,  and  the  conductor  losses  become  significant.  A 
careful  consideration  of  the  general  formula*  for  echo  area  shows  that,  if  the 
current  distribution  on  the  scatterer  is  not  changed  much,  instead  of  Eq.  (1) 


a 


3T 


GR 


rad 


Z.  +  ZT 

m  Jl, 


2 


(10) 


where  Rra{j  is  the  radiation  resistance.  The  maximum  echo  area  of  a  lossy 

scatterer  is  obtained  when  ZT  =  -jX.  ,  in  which  case  Eq.  (10)  gives 

L  m 


Rrad 

R  . 
rad 

12  9 

r  >■  G‘ 

Rrad  +  Rloss 

=  °0 

Rrad  +  Rloss 

a  « 
max 


(11) 
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where  ^oss  is  the  less  resistance  in  both  the  scatterer  and  the  load,  and  yQ  is 

the  loss-free  echo  area.  Similarly,  the  Q  of  the  scatterer,  which  is  determined 

by  the  resonance  curve  of  Z.  +  ZT ,  becomes 

m  Li 


R 


*  =  Qo  <R 


rad 


rad  +  Rloss 


(12) 


where  Qq  is  the  Q  of  the  loss -free  scatterer.  Forming  the  gain-bandwidth  prod¬ 
uct  of  Eq.  (7)  and  then  using  Eq.  (3),  one  finds 


4a  p=  nT^q  • 


(13) 


Hence,  the  gain -bandwidth  product  is  unchanged  by  losses,  assuming  that  the  current 
distribution  on  the  scatterer  is  unchanged.  The  primary  effect  of  losses  therefore 
is  to  reduce  the  echo  area  according  to  Eq.  (11)  and  to  reduce  the  Q  according  to 
Eq.  (12),  wiinout  changing  the  gain-bandwidth  product. 
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VIII.  A  Determination  of  the  Scattering  From  a 

Cavity- Backed  Plane  Surface 

J.  X.  Schindler.  1/Lt,  USAF* 
Microwave  Physics  Laboratory 
Air  Force  Cambridge  Research  Laboratories 
Bedford,  Massachusetts 

F.  V.  Schultz 
School  of  Electrical  Enginee.ing 
Purdue  University 
Lafayette,  Indiana 


1.  INTRODUCTION 

One  possible  technique  for  evaluating  the  scattered  return  from  surface -loaded, 
conducting,  scattering  bodies  consists  of  separating  the  scattered  field  into  two 
terms.  One  term  consists  of  the  scattered  field  from  the  body  with  all  loading 
structures  short-circuited  at  the  surface  of  the  body,  that  is,  the  unloaded  body. 

The  additional  term  represents  what  is  essentially  the  radiation  effects  of  the  aper¬ 
ture  that  couples  the  loading  structure  and  the  exterior  domain.  Of  course  the  inci¬ 
dent  field  plus  the  unloaded  return  gives  a  zero  tangential  electric  field  over  the 
whole  surface  of  the  assumed  perfectly  conducting  scattering  body.  The  additional 
scattered  field  term  is  then  necessary  to  account  for  the  non-zero  tangential  electric 
fields  induced  by  the  excitation  in  the  coupling  aperture. 


^Formerly  with  the  School  of  Electrical  Engineering,  Purdue  University,  Lafayette, 
Indiana. 


106 


With  this  technique  the  problem  of  analyzing  the  scattering  from  a  surface - 
loaded  body  is  broken  into  two  problems,  the  usual  unloaded  scattering  problem  for 
the  body  under  consideration  plus  the  problem  of  determining  the  radiation  from 
structures  equivalent  to  aperture  antennas  located  on  the  surface  of  the  scattering 
body.  The  aperture  antennas  are  of  course  excited  in  a  very  special  way  by  the 
external  radiation. 

A  key  problem  then  is  the  determination  of  how  the  aperture  fields  are  excited 
by  the  incident  radiation  and,  further,  how  the  magnitude  and  phase  of  the  aperture 
fields  depend  upon  the  cavity  parameters.  A  knowledge  of  this  latter  dependence 
allows  one  to  determine  how  the  cavity  may  be  made  to  control  the  reradiated  energy 
in  such  a  fashion  as  to  appropriately  modify  the  radar  re  turn  from  the  body. 

In  this  work  we  briefly  consider  the  analysis  of  the  scattering  from  a  highly 
idealized  cavity -loaded  structure.  The  structure  consists  of  an  infinite  plane  backed 
by  an  infinitely  long  rectangular  cavity  and  coupled  through  the  plane  by  a  narrow, 
infinitely  long  aperture  parallel  to  the  cavity  axis  We  do  not  approach  this  oroblem 
in  an  effort  to  minimize  the  return  from  an  infinite  screen;  rather,  it  is  believed 
that  the  results  obtained  for  this  idealized  structure  will  be  of  assistance  in  the 
analysis  of  more  practical  three-dimensional  bodies. 

For  example,  the  aperture  fields  obtained  may  approximate  those  excited  in  a 
cavity -loaded  finite  body,  if  the  principal  radii  of  curvature  of  the  scattering  surface 
are  large  at  the  point  of  loading  and  if  resonant  aperture  and  cavity  lengths  are 
avoided.  Ax  worst,  one  would  expect  the  same  qualitative  behavior  of  the  aperture 
fields  in  the  finite  body  as  is  exhibited  by  the  aperture  fields  in  the  present  problem. 
In  addition,  the  results  of  the  present  v/ork  will  provide  a  comparison  for  the  results 
of  simpler  though  more  approximate  schemes  for  the  calculation  of  the  aperture 
fields  in  finite  bodies. 


2.  THE  PROBLEM  AND  ITS  FORMULATION 

We  consider  the  scattering  of  a  plane  electromagnetic  wave  from  a  plane  screen 
of  infinite  extent  in  all  directions  and  backed  by  an  infinitely  long  rectangular  cavity. 
The  cavity  is  coupled  through  the  screen  by  an  infinitely  long  slit  parallel  to  the  axis 
of  the  cavity.  A  cutaway  view  of  the  structure  is  shown  in  Figure  1  and  a  cross 
sectional  view  is  shown  in  Figure  2.  The  structure  is  assumed  to  be  excited  by  a 
monochromatic  plane  wave  (with  time  variation  exp  (jwt)  )  having  an  arbitrary 
polarization  and  arbitrary  polar  and  azimuthal  angles  of  incidence,  6Q  and  4>o  , 
respectively.  (See  Figure  1. ) 

The  portion  of  space  where  electromagnetic  fields  may  exist  is  divided  naturally 
into  two  regions:  the  cavity,  given  by  -b  <  x  <  b  ,  -a  <  y  <  0  ,  and  the  semi -infinite 
half  space,  where  y  >  0  . 
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Figure  1.  Typical  Cross  Section  Taken  Figure  2.  Cross-Sectional  View 

between  Two  Constant  z  Planes  Showing  of  the  Cavity  Showing  the  Dimen- 

Cavity -Backed  Infinite  Plane  Screen  and  sion  Parameters  of  the  Structure 

Direction  of  Incidence  of  Excitation 

In  each  of  these  regions,  the  electromagnetic  field  may  be  expressed  as  a  linear 
superposition  of  divergenceless  solutions  to  the  homogeneous  vector  Helmholtz 
equation. 

In  the  cavity  region,  the  fields  are  expressed  as  the  superposition  of  a  discrete, 
doubly  infinite  set  of  modes,  both  transverse  magnetic  to  z  and  transverse  electric 
to  z  modes  being  included.  To  these  expressions  we  apply  the  boundary  conditions 
which  require  zero  tangential  electric  field  on  the  cavity  surfaces,  y  =  -a  ,  -b  £  x 
<  b  and  x  =  ±b  ,  -a5y<  0,  The  resultant  expressions  for  the  cavity  electric 
field,  Ec  ,  and  the  cavity  magnetic  field,  Hc  ,  are 
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EC 


OO  — 

E  !x  (jfl  dncos0o-mncn  )cosf§(x  +  b) 
1=0  L  \  ' 

+  V(imndncos8o  +  -lcn)  si 


sin  rn  (y  +  a) 
nv 


cos  0Q  +  -§g  cn )  sin  -fg  (x  *  b)  cos  mn(y  +  a) 


+  i  8  sin  0  d  sin  (x  +  b)  sin  m  (y  +  a) 


(1) 


Hl 


[!>i(mndn-i|icncose<.)sin2b 
+  V  ("5rdn  •  1  mn  cn  003  8o)  cos  2b 


(x  +  b)  cos  m  (y  +  a) 
n 


(x  +  b)  sin  m  (y  +  a) 

R 


+  i  0  sin  8  c  cos  (x  +  b)  cos  m  (y  +  a) 

*7  J  r\  n  * 


] 


(2) 


Here  cn  and  dn  ,  n=0, 1, . . . ,  are  unknown  complex  constants,  Z  =  y H0/eQ  , 

1/2 

.  The  assumed  harmonic  time 

and  z  variation,  exp  (ut  +  kz  cos  0Q)J  >  has  been  suppressed  from  these  ex¬ 
pressions. 

In  the  exterior  half  space  the  total  fields  are  written  as  the  sum  of  three  terms 
in  accordance  with  the  discussion  given  in  the  introduction.  The  first  term  is  the 
incident  field,  E1  ,  which  is  known.  The  second  term  is  the  unloaded  scattered 

g 

field,  E  1  ,  which  is  the  field  which  would  exist  when  the  aperture  i~d:£x:£f+d 
is  replaced  by  a  perfectly  conducting  sheet.  The  unloaded  scattered  field  is  simply 
the  field  specularly  reflected  from  an  infinite  screen  and  is  thus  known  from  ele¬ 
mentary  calculations.  The  third  term  in  the  field  expression  is  the  loaded  scattered 

field.  This  field  is  that  which  is  produced  by  the  existence  of  non-zero  aperture 

— Sg  -So 

fields.  It  is  these  field  quantities,  denoted  by  E  and  H  ,  that  we  wish  to  com¬ 
pute. 

Let  us  concentrate  on  expanding  these  latter  field  quantities.  It  is  known  that 
if  is  a  solution  to  the  homogeneous  scalar  Helmholtz  equation  and  a  is  any 
constant  vector,  the  vector  functions  M  and  N  ,  defined  by 


0  J  k  sin  0  ,  (k  =  2jt/X),  and  m 


n  ’  f 2  ‘(I)2" 


M  =  Vx(a$) 


(3) 


N  =  1/k  (V  X  M )  , 


(4) 
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are  purely  solenoiaal  solutions  to  the  vector  Helmholtz  equation.  Moreover,  there 
exist  collections  of  the  functions  of  Eqs.  (3)  and  (4)  which  are  a  basis  for  the  expan¬ 
sion  for  the  most  general  solenoidal  solution  to  the  vector  Helmholtz  equation. 

With  0  of  the  form  (suppressing  the  harmonic  time  dependence) 


<Mx,y,  z)  =  *(x,y)  exp  (jkz  cos  0Q)  , 
4'  may  be  written  as 


(5> 


*fx,y) 


oo 


exp  {j  [X(a)y  -  ox]  }  da  , 


(6) 


where 


c(a) 

is  an  unknown  to  be  determined,  and 

(  j(a2  -p2)112 

\a\  >p 

X(a) 

) 

1  (f2  -  a2)1'2 

\o  \  —  P  . 

It  may  be  readily  veriiied  that  Eq.  (6)  substituted  into  Eq.  (5)  gives  a  result  that 

obeys  the  scalar  Helmholtz  equation.  Further,  as  the  cylindrical  radius 
2  2  1/2 

p  =  (x  +  y  )  '  “  becomes  large,  Eq.  (5)  appears  in  the  form  of  a  diverging  cylin¬ 
drical  wave  and  thus  obeys  the  radiation  condition. 

We  now  use  in  Eqs.  (3)  and  (4)  the  result  of  substituting  Eq.  (6)  into  Eq.  (5). 

Using  both  transverse  electric  to  z  and  transverse  magnetic  to  z  field  expres  - 

-So  -So 

sxons,  the  quantities  E  and  H  take  the  form 


ES2 


=  L.  \  I  J  c(a)  +  a  D(a)  cos  0  ~j  exp  (  j  [X(a)y  -  ax]  }  da 
x  42*  L  °-J 

-oo 

CO 

- —  jj  Qja  c(a)  -  X(a)  D(a)  cos  0QJ  exp  f  j[X(a)y  -  ax]}  da 


y  ^/2F 
P  sin  9 


+  i_ 


4W 


—  Co(a)  exp  (j(X(a)y  -  ax]  }  da 


O 

-CC. 


(7) 
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H** 


7. 


—  SO 

c(a)  cos  6  +  jX(a)  D(a)]exp  {  j  (X(a)y  -  ax]}  da 
V  2*  ~  <-  -  J 


-ao 

CO 


+  1  — - — \  [-X(a)  c(a)  cos  0  +  ja  D(a)J  exp  (j[X(a)y  -  ax]}  da 

y  '^27  J  ° 


+  i 


0  sin  9C 

*/2 IT 


c(a)  exp  {j[X(a)y  -  ax]}  da 


(8) 


Here  D(a)  is  also  unknown  and  is  to  be  determined. 

Our  objective  at  this  point  is  to  use  the  results  of  Eqs.  (1),  (2),  (7),  and  (8)  to 
form  integral  equations  for  the  x  and  z  components  of  the  aperture  electric  field. 
The  technique  used  is  similar  to  that  proposed  by  Lewin*  for  a  class  of  waveguide 
discontinuity  problems. 

If  Ez(x)  denotes  the  unknown  z  component  of  the  electric  field  in  the  aperture, 
then  the  boundary  condition  requiring  continuity  of  the  tangential  electric  fields  at 
y  =  0  interface  gives 


5  (x)  i-d  =£  x  <  f+d 

— z 

0  otherwise 

when  -b  £  x  £  b  .  Substituting  for  Ec  and  using  the  orthogonality  of  the  sine 
functions  gives 


/ 


(9) 


flb  sin  6  sin  m  ad 
o  P  P 


f+d 

=  ^  Ez(x)  sin  -|£-  (x+b)  dx 
i-d 

f+d 

=  ~  ^  F(x)  cos  -gg-  (x+b)  dx  , 
f-d 


where 


?(x)  5  -t  £z{x)  ' 


(10) 


The  last  equality  follows  from  an  integration  by  parts  a3  suggested  by  the  work  of 
Lewin. 


4 


1 


Ill 


In  a  similar  fashion  the  continuity  of  E  relation  gives 


r-i  -S1 
E1  +  E  1 


\  Iz(x) 

I  0 


i-d  f  x  £  i+d 


otherwise 


(11) 


for  all  x  .  Note  that  (e*+E  /  '  *zjy=0  =  ®  from  the  definition  of  these  quanti¬ 

ties.  Substituting  for  E  2  from  Eq.  (7)  into  this  result  and  computing  the  Fourier 
transform  of  the  equation  gives 


I  +  d 

p  sin  0  D(a)  =  — - —  \  E  (x)  exp  (jar>.)  dx 

i+d 


<v/2 T 


jjj  J  F(x)  exp  (jax)  dxj  . 


(12) 


The  last  equality  follows  from  an  integration  by  parts.  Also,  P  denotes  that  the 
Cauchy  principal  value  of  the  expression  is  to  be  taken  upon  an  integration  with 
respect  to  c  . 

In  addition  to  the  conditions  above,  let  us  require  the  continuity  of  the  tangential 
components  of  the  magnetic  field  in  the  aperture,  i-d  £  x  £  f+d  .  That  is,  we 
require 


(i?+HSi+HS2) 

*  1z|y=o 

=  Hc  • 

lz 

II 

o 

i-d  <  x  <  l+d 

(13) 

(h‘  +  HS1  +  HS2) 

1  •  I  1  n 

x|y=0 

O 

IS) 

ii 

5x 

|y=o 

i-d  <  x  <  i+d  . 

(14) 

Note  that  due  to  the  definition  of  HS1  ,  (H^+  HS1)  •  I  I  n  =  2H*  •  I  I  ..  and 
s  z  |y=u  zjy=u 

(H*  +  H  1)  •  Ix|y_Q  =  2 H*  -  ixjy=Q  •  After  substituting  for  HS2  from  Eq.  (8) 

and  for  H  from  Eq.  (2),  we  note  that  the  second  equation  added  to 
-j  cos  0Q/p  sin  0Q  times  the  x-derivaiive  of  the  first  equation  gives  a  result  in¬ 
volving  only  the  unknown  quantities  D(a)  and  dn,  n -0,  1,  2, . . .  .  Substituting  for 
these  unknowns  from  Eqs.  (12)  and  (10)  respectively  gives  an  equation  in  the  single 
unknown  F  .  A  rearrangement  of  this  result  leads  to  the  integral  equation 


112 


+  1 

§  F(Xd+ i)  [G(X-t)  + H(X,r)]  dX  =  Eq  (t d  +  i )  -1<t<+1,  (15) 

-1 

where 


G(X-t)  =  j  P  y  exP  [  jord(X-r)  ]  da 


Em 

— —  cot  mn  a  sin  ^  (rd+  i  +  b)  cos  -g£.  (;v(j  +  i  +  b) 


EQ(rd+l)  =  -j  2np  sin  Qq  Zq 


iHi  •  jxiy=o  - 


j  cot  9  _  . 

O  d  r;i 


'  *z |y=0 


The  solution  to  this  integral  equation  ,  F ,  along  with  the  relation 


(x-i)/d 

Ez(x)  =  d  j*  F(Xd  +  i)dX  , 
-1 


makes  possible  the  determination  of  the  z  component  of  the  aperture  electric 
field.  Remember  in  connection  with  Eq.  (16)  that  F  =  d/dx  Ez(x)  and  that 
Ez(l-d)  =  0. 

Let  us  consider  now  the  derivation  of  a  second  integral  equation,  the  solution 
of  which  is  related  to  the  x  component  of  the  aperture  field.  Let  Ex(x)  denote 
the  x  component  of  the  electric  field  in  the  aperture.  Continuity  of  the  x  compo¬ 
nent  of  the  electric  field  at  y  =  0  then  gives 


f - d  <  x  s  Jt+d 

otherwise 


(17) 


for  -b  ^  x  ^  b  .Substituting  for  Ec  from  Eq.  (1)  in  this  expression  and  adding 
-jd  cot  0Q  times  the  x  derivative  of  Eq.  (9)  [with  Ec  substituted  from  Eq,  (l)[j 
gives 
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-0d  )  in  c  sin  in 
/  /  n  n  n 


a  cos  5r-  (x+b) 


^  2(x) 
/  0 


/-d<x5 /+d 


otherwise 


for  -b  <  x  <  b  .  Here 


G(x)  =  /3d  E  (x)  -  jd  cot  0  f(x)  . 


Using  the  orthogonality  properties  of  the  Fourier  cosine  series  gives 


(en  -  3)  m 


r 

b  0d  c  sin  in  a  =\  G(x)  cos  -Sr-  (x+b)  dx  , 
n  n  nj—  2b 


where 


(  1  n  =  0 

I  2  ntC 


In  a  similar  fashion,  continuity  of  E  ai  the  y=0  plane  requires 


((E1  +  ES1  +  ES2)  • 


,(x)  I  -  d  :=  x  £  f+d 


otherwise 


_i  —Si  _  —So 

for  all  x  .  Again,  noting  that  (E  +  E  )•  ix|y_0  =  0,  substituting  for  E  “  from 

Eq.  (7),  and  adding  -jd  cot  6q  times  the  x  derivative  of  Eq.  (11)  with  ES2  sub¬ 
stituted  from  Eq.  (7)  gives 


°°  l  G(x)  l- d  <  x  <  l+d 

d  C  )  ~ 

—  \  X(a)  c(a)  exp  (-jorx)  da  =  - 

Trr  J  *  ) 


otherwise  . 


Upon  taking  the  Fourier  transform  of  this  result,  there  results 


ITU 

j0d\(a)  c(a)  =  — - — \  G(x)  exp  (jax)  dx  . 

kI  ~ 
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~.So 

At  this  point  let  us  consider  Eq.  (13).  Substituting  for  H  from  Eq.  (8)  and 
for  Hc  from  Eq.  (2)  and  using  the  fact  that  (H*  +  H  1)  •  iz  |y_Q  =  2  K1  ■  iz  ,  we 

find  a  relation  in  the  unknowns  c  ,  n=0,  1,2, _ _  and  c(or)  .  Substituting  for  the 

c  's  from  Eq.  (20)  and  for  c(o;)  from  Eq.  (23)  and  rearranging  the  result  gives 
the  integral  equation 

+  1 

y  G(\d+i)  CI(X-t)  +  J(\,  t)]  d\  =  E1(rd+i)  ,  -1  =£  r  <  +1  ,  (24) 

-1 

where 


I(X“T)  =  t  I  x§)  exP  D0fd{X'7)D  d“ 


J(X,r)  = 


-  j  4  Z  ,  j 

E1(rd+l )  =  ~r~—  “  (h1  ■  Iz(rd+i)|y=0  j 


sin  0Q  9t 


Due  to  the  relation  defining  G  ,  Eq.  (19),  the  solution  to  Eq.  (24)  along  with 
the  solution  F  of  Eq.  (15)  makes  possible  the  calculation  of  E  . 


3.  PROBLEM  SOLUTION 

The  complexity  of  the  integral  Eqs.  (15)  and  (24)  seems  to  preclude  any  exact 
solution.  Instead,  we  consider  a  standard  perturbation  technique  in  order  to  obtain 
an  approximate  solution  under  the  assumption  of  a  small  aperture  width  in  wave  - 
lengths. 

Consider  the  integral  equation 

Jf(x)  A(x,y)  dx  =  g(y) 

L 


(25) 
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where  A  and  g  are  known  functions  and  f  is  an  unknown  function  to  be  de¬ 
termined.  Expanding  A  and  g  in  a  series  in  the  dimensionless  parameter  e 
and  assuming  a  similar  expansion  exists  for  the  unknown  f  results  in 


A(x,y)  =  T,  ertAn(x,y)  . 
n=0 

OO 

g(y)  =  Y  * 

n=0 

OO 

f(x>  =  Y/ gn  * 

n=0 


Substituting  these  expansions  into  Eq.  (25)  and  assuming  the  equality  of  the  coeffi¬ 
cients  of  like  powers  of  e  yields  the  following  sequence  of  integral  equations: 


§  f°(x)  A°(x,  y)  dx  =  g°(y) 

L 

§  f1(x)  A°(x,y)  dx  =  gJ(y)  -  ^  f°(x)  A*(x,y)  dx 
L  L 

and  so  forth.  Thus,  if  an  integral  equation  having  a  kernel  of  the  form  of  A^  is 

solvable,  then  each  integral  equation  in  the  above  sequence  is  solvable  for  the  suc- 

0  12 

cessive  approximations  f  ,  f  ,  f  ,  .  .  .  .  It  should  be  noted  that  the  right-hand 
side  of  each  equation  is  a  known  function  provided  th »  previous  equations  in  the 
sequence  have  been  solved. 

Let  us  apply  this  perturbation  technique  to  the  integral  Eq,  (15)  for  F  .  The 
small  parameter  e  in  this  development  will  be  taken  to  be  fid  .  To  apply  the 
technique,  we  require  series  expansions  in  fid  for  the  (unctions  G  ,  H  ,  and  Eq 
defined  following  Eq.  (15).  A  detailed  evaluation  of  G(\-t)  reveals  that 


G(\-t)  =  *p[x3f]+  0(fid)2  . 
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Also, 


H(A,r)  =  -  P  0(/3d> 


and 


E0(tcH-j£ )  =  j  sin  4>0\’  e  exp  (j pi  cos  $Q)  +  0(pd)  . 


Here  e  is  a  vector  in  the  direction  of  the  E  polarization  of  the  incident  wave  and 
of  magnitude  equal  to  the  incident  electric  field. 

Substituting  these  results,  along  with  the  assumed  pd  expansion  for  F  , 


F(Ad+X)  =  F°(Ad4-i)  +  0d  F1(Xd-5-i)+  .  .  .  , 


(26) 


into  Eq.  (15)  and  equating  the  coefficients  of  (pd/  gives 


+1  n , 

P  r  F°(Ad+i>  dA 

J  - A~f -  =  sin  *z  *  6  exp  ^  COS  -1  —  7  — +  1  '  (27) 


-i 


This  is  a  singular  integral  equation  of  the  form  considered  by  Tricomi.  The 
solution  to  Eq.  (27)  is 


Fu(Ad+t )  = 


-jp\  sin  <pQ  Iz  •  e  exp  (j pi  cos  4>n) 


(1  -  A2)l/2 


where  we  have  used  the  fact  that 


+1 


F  (Ad+Jt)  dA  =  0 


This  latter  result  follows  since 


0  = 


E  (i  +  d)  -  E  (i  -d) 
* 


+  1  +1 
=  j*  F(Ad+t)  dA  =  ^  Fn(Ad+i  )  dA  , 

-1  n=0  -1 
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or 


+  1 

jFn(AdH)d\  =  0  n  =  0,  ,2,...  . 

-1 

Here  the  first  equality  follows  since  the  z  component  of  the  aperture  field  must 
be  zero  at  the  edges  of  the  aperture.  The  second  equality  follows  from  the  defini¬ 
tion  of  F  given  following  Eq.  (10),  and  the  third  eqi  av.:y  follows  from  a  substitu¬ 
tion  of  the  assumed  /Jd  expansion  for  F . 

A  continuation  of  the  perturbation  scheme  yields  integral  equations  of  the  form 
of  Eq.  (27)  in  the  unknowns  Fn,  n=l,  2,  3, . . .  .  These  results  substituted  into 
Eq.  (26)  along  with  the  condition  Eq.  (16)  gives  for  the  aperture  field 


E  (A d+l)  =  /Jd(l-A^)1^2  sin*  i  ■  e  exo  (j fit  cos  *  )  j  j 
— z  O  Z  *  O’ 


D 


■  1  <  A  <  +1 


cos  0(3d)“J 
(28) 


Applying  this  same  perturbation  technique  to  Eq.  (24)  we  arrive  at  a  sequence 
of  integral  equations  of  the  same  form  as  Eq.  (27).  Without  going  into  detail* 
there  results  for  the  x  component  of  the  aperture  electric  field 


Ex(Ad-H ) 


2  i<£o  •  e  exp  (j0l  cos  <f>Q) 
/IdjrU-A2)1/2  Ao 


+  0  (0d)  . 


(29) 


Here  I^o  is  the  unit  vector  i^  evaluated  in  the  direction  from  which  the  incident 
wave  arrives.  Also, 


where  y  =  Euler -Mascheroni  constant  =  C.5772157  and  yn  =  £l  -  (2/?b/n7r)  j  . 

The  results  of  Eqs.  (28)  and  (29)  for  the  z  and  x  components  of  the  aperture 
electric  field  are  the  main  results  to  be  obtained  here.  It  should  be  noted,  however, 

A  detailed  analysis  of  this  and  other  aspects  of  the  problem  will  be  included  in  an 
AFCR  L  report  in  preparation. 
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that  by  applying  the  method  of  steepest  descent  to  the  Fourier  integrals  in  the  ex- 
— s<>  —so 

pressions  for  E  and  H  ,  Eqs.  (7)  and  (8)  respectively,  it  is  possible  to 
evaluate  the  far  scattered  fields  produced  by  the  non-zero  aperture  fields  Ez  and 
Ex  .  We  will  not  consider  these  results  here. 

4.  DISCUSSION  OF  RESULTS 

In  this  section  we  consider  Eqs.  (28)  and  (29)  for  the  z  and  x  components  of 
the  electric  field  in  th,.-  aperture. 

Consider  the  result  for  Ez  first.  We  note  that  this  field  component  is  propor¬ 
tional  to  the  z  component  of  the  incident  electric  field  and  is  independent  of  all 
other  components.  Further,  this  field  is  0(£d)  and  hence  for  a  narrow  aperture 
will  be  small.  A  more  complete  analysis  reveals  that  the  first  term  in  the  solution 
to  involve  the  cavity  parameters  a  ,  b  ,  and  i  is  O(0d)  .  Thus  it  would  seem 
that  for  a  narrow  aperture  the  cavity  will  do  little  to  change  or  control  the  z  com¬ 
ponent  of  the  aperture  field. 

Finally,  it  is  interesting  to  note  that  near  the  aperture  edge,  X  —  =1  ,  |ezI 
varies  according  to  the  square  root  of  the  distance  from  the  edge.  This  is  the  known 

3 

edge  condition  for  electric  field  components  parallel  to  a  sharp  edge. 

Let  us  consider  also  the  more  interesting  result  of  Eq.  (29)  for  E^  .  We  note 

that  the  dominant  term  for  £d  small  is  the  first  term,  it  being  0(l/£d)  .  This 

dominant  term  has  the  same  spatial  variation  as  the  static  electric  field  which 

would  be  produced  in  an  infinitely  long  narrow  slit  by  a  potential  difference  across 

the  slit.  Near  the  edges,  X  —  ±1  ,  this  component  of  the  field  is  proportional  to 
-1/2 

p  where  p  is  the  distance  to  the  edge.  Again,  this  is  the  known  edge  condi - 

3 

tion  for  electric  field  components  perpendicular  to  a  sharp  edge. 

It  is  of  interest  to  consider  Ex  at  the  center  of  the  aperture  (X=0)  when 
l  =0  ,  that  is,  Ex(0)  ,  as  a  function  of  the  cavity  parameters  a  and  b  .  A  detailed 
analysis  of  Eq.  (29)  for  0Q  =  90°  and  4>Q  =  45°  reveals  that  Ex(0)  =  0  for  ka  =  mr , 

n=0,  1,  2, _  as  long  as  kb*  pz  ,  p=i,  2, _  and  kd  *  0  .  Thus  the  cavity  appears 

to  short  circuit  the  aperture  when  the  cavity  depth  is  adjusted  to  integer  multiples 
of  a  half  wavelength. 

There  also  occur  combinations  of  the  cavity  dimensions  a  and  b  where  what 
might  be  called  "anti-resonances"  occur.  Here  the  dominant  first  term  in  the 
expression  for  Ex  is  a  maximum  in  magnitude. 

These  "anti -resonance"  dimensions  are  not  odd  integer  multiples  of  a  quarter 
wavelength  as  might  be  expected.  Rather,  these  dimensions  are  given  by  the  condi¬ 
tion  that 
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In?  £0  -  0  . 

Further,  it  can  be  shown  that  this  condition  occurs  in  the  range  0  <  a/\  <  1/2 

when  q  <  2b/A  <  q  +  1/2  ,  q=0,  1, 2, _  and  in  th.=  range  u/A  >  1  j  '  for  all  values 

2d 

of  2b/ A  .  A  numerical  analysis  for  the  conditions  -r-  =  0.  1  ,  *=0,0.  =  90°  , 

A 

o  =  45°  reveals  that  under  these  conditions  the  "anti -resonance'-  values  of  a/A 
o 

tend  to  n/2  ,  n*l,  2.  3, _  as  2b /a  becomes  large.  Also,  the  "anti -resonances" 

become  sharper  as  2b/ A  becomes  larger  since  the  magnitude  of  the  field  varies 
rapidly  to  zero  at  a/A  =  n/2  . 
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IX.  Some  Concepts  for  Reducing  Reflectivity 

From  Antenna  Apertures 

Edwin  M.  Tufnrr 
Electromagnetic  Environment  Branch 
Air  Force  Systems  Command 
Aeronautical  Systems  Division 
ttrighl-Palterson  AFB,  Ohio 


An  elementary  analysis  of  radar  reflectivity  will  generally  show  that  for  high¬ 
speed  aerospace  vehicles  the  most  difficult  surfaces  to  conceal  are  those  occupied 
by  the  antenna  apertures.  This  is  true  because  the  antennas  are  generally  focused, 
reciprocal,  and  linear  devices,  and  they  must  be  strategically  placed  to  perform 
their  assigned  functions.  Within  the  confines  of  these  assumptions,  one  has  a 
number  of  readily  available  alternatives  to  reduce  the  backscattering  from  antenna 
apertures;  the  alternatives  include  terminal  loading,  inductive  loading,  absorptive 
loading,  and  choosing  antenna  designs  which  are  inherently  low  in  reflectivity. 

Often  a  more  productive  solution  for  the  backscattering  problem  may  be  found 
from  a  study  of  the  flight  trajectory  and  the  operational  requirements  of  the  sys¬ 
tem.  For  instance,  an  antenna  configuration  may  be  used  over  only  a  small  part 
of  the  flight  path,  and  it  may  be  hidden  from  radar  view  throughout  the  remainder 
of  the  trajectory  without  being  operationally  objectionable.  This  may  be  accom¬ 
plished  by  either  changing  the  antenna  patterns  to  produce  a  beam  pointing  in  some 
innocuous  direction,  or  by  breaking  *he  aperture  up  into  small  segments,  each  of 
which  will  resonate  at  frequencies  above  or  below  the  operating  frequencies  of 
the  tracking  radars.  These  objectives  may  be  achieved  by  the  use  of  biased  diodes, 
non-reciprocal  ferrites  in  waveguides,  switches,  or  by  relay-operated  shorts  that 
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will  change  the  resonance  frequencies  or  the  patterns  of  the  apertures. 

Examples  of  the  techniques  suggested  herein  are  shown  b-_.ow  in  Figures  1 

to  4. 


Figure  1.  A  Horn  Aperture  Shorted  by  Means  of  a  Diode  Switch  Reduces 
Backscattering  to  Approximately  That  of  a  Flat  Plate 


Figure  2.  A  Conical  Spiral  Antenna  Showing  Fatterns  for  Both  In-f<‘ase  and 
Anti -Phase  Connections.  A  change  in  patterns  can  be  achieved  by  means  of 
placing  a  short  across  the  feed  terminals. 
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Figure  3.  A  Dipole  Antenna  With  Forward  Biased  Diode  Switches.  A  change  in 
dn  bias  across  diodes  will  segment  the  aperture. 


Figure  4.  An  Open-Ended  Wave -Guide  Radiator  Showing  Non-Reciprocal  Ferrite 
Phasing  Element.  The  angle  of  divergence  between  the  transmitted  and  received 
beams  is  proportional  to  the  magnetic  field  across  ferrite  phase  shifter. 
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X.  Radar  Cross  Section  of  Perfectly  Conducting, 
Dielectric  and  Dielectrically  Clad  Infinite 
Cylinders  at  Arbitrary  Incidence* 

Abdelnour  S.  Thomas 
A.  S.  Thomas,  Inc. 
Restv.ood,  Massachusetts 


1.  INTRODUCTION 

1  —  7 

The  work  reported  in  the  literature*  on  dielectrically  clad  cylinders  prima¬ 
rily  treats  normal  incidence  and  in  most  cases  is  restricted  to  the  geometrical 
optics  region.  As  far  as  could  be  ascertained  by  the  author,  the  dielectrically  clad 
infinite  cylinder  at  normal  (broadside)  incidence  for  ?  TM  incident  field  was  first 

7 

treated  by  C.  C.  Tang  with  computed  values  of  normalized  radar  cross  section 
experimentally  verified.  The  more  general  '.ase  of  arbitrary  incidence  was  treated 

g 

by  J.  R.  Wait  for  a  homogeneous  non-metallic  cylinder  for  an  incident  TM  mode. 
He,  however,  did  not  present  any  numerical  data.  The  dielectrically  clad  cylinder, 
following  the  work  of  Wait  and  Tang,  is  extended  here  to  the  more  general  case  of 
oblique  incidence  of  either  TM  or  TE  modes  for  N  concentric  layers.  The  normal¬ 
ized  radar  cross  section  in  the  plane  of  incidence  has  been  computed  for  normal 
and  oblique  incidence. 

It  is  believed  that  a  comparison  of  the  normalized  radar  cross  section  of  the 
dielectrically  clad  and  solid  dielectric  cylinders  for  small  radii  with  that  of  the 


*This  work  was  sponsored  by  Lincoln  Laboratory,  Subcontract  No.  265,  under 
Prime  Contract  No.  AF19(604)-7400. 
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perfectly  conducting  cylinder  would  provide  some  insight  into  the  problem  of  re¬ 
ducing  the  radar  cross  section  of  objects  of  small  geometrical  cross  section. 


2.  FORMULATION 


2. 1  Genera! 

The  equations  of  the  incident,  internal,  and  scattered  fields  at  arbitrary  inci¬ 
dence  are  presented  here  fcr  a  cylinder  consisting  of  N  homogeneous,  concentric, 
non-metallic  cylindrical  layers.  The  expressions  for  the  scattering  coefficients 
for  a  single  layer  about  a  perfectly  conducting  core  with  an  incident  TM  field 
(Hz  =  0)  are  derived,  and  then  the  expressions  for  incident  TE  mode  are  given. 
These  general  expi’essions  reduce  to  those  for  the  solid  homogeneous  dielectric 
cylinders  as  well  as  to  tae  solid  metallic  cylinder. 

Consider  a  plane-wave  incident  on  an  infinite  cylindrical  obstacle  at  an  angle 
4'  with  respect  to  the  normal  to  the  axis  of  the  cylinder  to  be  a  pure  TM  wave. 
Using  the  coordinate  system  given  in  Figure  1,  the  incident,  internal,  and  scattered 
fields  may  be  written  as  follows: 

(1)  Incident  Fields 


Einc  s  Eo  cos  *„£  ittjAr,F, 


,  inc 


=  E  cos*  >  3 


A 


-nk  sin*' 
o  o 


;Jn(3r)Fn 
no  n 


4 


H 


inc 

<?> 


E  cos  * 
o  o 


( 2)  Internal  Fields 


(1) 

(2) 

(3) 

(4) 


E 


m 


Dn(*nHn 


(1) 


(B  R  )  +  BmH(2)(/3  R  ) 
nr  n  n  'in  m 


K 


(5) 


128 


where  it  is  understood  that  the  summation  is  from  n  =  -  “to  n=°° 

En  =  magnitude  of  the  incident  field 

=  angle  of  incidence  with  respect  to  the  axis  of  the  cylinder 

8  -  k  cos  if  ,  with  in  -  0,  1,  2, ... .  N  relating  to  the  mth  layer, 

with  m  -  0  free  space  external  the  cylinder 


f  =  operating  frequency 

-  normalized  permeability  of  the  mth  layer 

=  normalized  permittivity  of  the  mth  layer 

Fn  =  exp  [  -j(n<t>  -  km  z  sin  *m)  ] 

Am  bw,  Cm  -  arbitrary  constants  relating  to  the  TM  mode 
n  n  n  * 

Am,  B™,  C™  =  arbitrary  constants  relating  to  the  assumed  internal  and 
n  n  n 

scattered  TE  mode 

J  ,  Hn  Hn2^  =  -T3es;5el  and  Kankel  functions  with  the  primes  relating  to 
differentiation  with  respect  to  the  argument. 
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2.2  Single  Layer  about  a  Perfectly  Conducting  Cylinder 

Nov/,  by  equating  the  internal  fields  to  the  sum  of  the  inc.dent  and  scattered 
fields  at  the  air  interface: 


E  cos*  )+C  h‘2)(x  )  =  A  H(1)(x.)+B 

o  ono  nn  o  nn  l  nn  l 


(13) 


nk  sin* 
o 


n*  r 
— £  E 
o  L  ° 


cos*  J  (x )  +  C 
ono 


H{2)(x  )] 
n  n  o  J 


-j-^  C  H<2)'(x  ) 
J  p  n  n  o 
o 


uk  sin  * 
o  o 


P,  x 


1*1 


C  H^2)(x  )  =  A  }  (x. )  +  B  (x. ), 
nn  o  nn  l  nn  l 


(15) 


and 


•>  2 
Jko 


i- 


1  nk  sin*  _ 

E  cos  *  J '  (x  )  +  C  H ( x  •  )  +  C  (x  J 

o  ono  nn  o  j  n  n  o 


o  o 


nk  sin  * 
o 


■S^*-  [ a  H(1 }  (x  )  +  B  H( 2)  (x  )]  +  — -g-  [a  H<1>'(x1)  +  Bn H*2)  '(x  )]  (16) 
jXj  [_nn  1  nn  1J  w/i^/3^|_n  n  1  n  n  1J 


where  *j  =  angle  wave  makes  with  normal  in  coating 

x  =  p  R 
o  o  o 

R  =  outer  radius 
o 

X1  =^lRo 

and  setting  at  the  metallic  interlace 


Bn 

T~ 

n 


H(l)  (xJ 
n  l 

H(2)>(x) 
n  & 


(17) 


and 


n 


H(1)"(x2) 


(18) 


the  Cn  and  Cn  may  be  readily  obtained  for  the  incident  TM  mode,  giving: 
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,TM 

'r. 


PN-[qnH|f>(xo)]2 


and 


_  as  (h!“)  V)  J  (x  )  -  H{n2)  (x  )  J'  (x  j] 
«  .  ti  ol  n  ono  n  onoj 


C  = 
n 


PN-kHn'(lto>) 


(19) 


(20) 


where:  N  =  s.  H^^(x  )D-  s  (x  ) 

1  n  o  on  o 

Pai'lH(„2)<VL'soHn,'(V 

H(1)'(Xl)  H(2)  (xJ  -hJ.2)'(x.)  H<1}  (xJ 
j-.  _ _ n _ in  £ _ n  l  n  i 

3JP (x  )  H{2)  (xJ  -  H(2)  (x.)  H(l)  (x,) 
n  1  n  l  n  l  n  2 

H*l,*{jc, )  H< 2) '(xJ  - H*2) '(x.)  H*1 1  '(x2) 

L  ’  imr,,,)  5^5 -h^Ix',)  h^Tx,)' 


s„  = 


1 


'e  -  sin  'J' 


So  cos  4f' 


r,  = 


~\Je  -  sin2  *£ 


nX  sin  ^  / 

O  I 


%  2r  RQ  ^ 


e-sin^’I'  cos  f 
o  o 


In  the  plana  of  incidence  (<t>  =  0°),  cr  may  be  written 


2A.cos  ¥ 


T, 


C  (-!)"< 


(21) 


Now  note  that  the  Cn>  Eq.  (20),  may  be  written  as  follows: 


nX  sin  'i' 


Cn  = 


“27 TT"  Gn' 


(22) 
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where 


hence,  for 

n  =  0 

Co=0 

and 

C  +C  „  =  0  for  n  |  0. 

Therefore,  in  the  plane  of  incidence,  the  orthogonal  roode  is  identically  zero. 

However,  the  scattered  field  in  planes  other  than  the  plane  of  incidence  will 

contain  the  orthogonal  mode  which  will  increase  as  a  function  of  the  angle  '1'  and 

the  angle  6  of  the  plane  of  observation,  with  a  maximum  in  the  plane  6  =  90°. 

Also,  if  the  axis  of  the  cylinder  is  skewed  with  respect  to  the  plane  of  incidence, 

the  scattered  field  in  the  plane  of  incidence  will  contain  the  orthogonal  mode. 

In  the  plane  of  incidence,  the  C  may  be  obtained  for  the  TE  incident  mode 
TE  n 

(Cn  )  by  replacing  E  by  H  and  H  by  -E,  e  by  p  and  p  by  e  throughout,  giving: 


where 


MN-qV(x )H<2)(x ) 

_  n  n  o  n  o 


(23) 


M*rlJn(xo)L"soJn(xo)* 
and  all  other  symbols  are  as  previously  defined. 

For  normal  incidence,  V  =  0°,  qn  =  0.  0.  Eqs.  (19)  and  (23)  become  the  well-known 
equations  for  normal  incidence 


,TM  _  V 

n  ""IT' 


and 


~TE  =  .  M 
un  P  ' 


(24) 


(25) 


2.3  Solid  Dielectric  Cylinder 

If  Xg  is  taken  equal  to  zero,  that  is,  a  solid,  homogeneous,  non-metallic 
cylinder,  Eqs.  (19)  and  (23)  immediately  reduce  to  those  for  the  solid  dielectric 
cylinder  at  oblique  incidence,  with 
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PTM  = 

VP  -  q“J  (x  )  H(  2)  (at)  J“  (x, ) 
n  n  o  n  o  n  1 

n 

PN  -  [qXX’ W] 2 

p.TE  _ 

MN-^Jn(VHn,(lto>Jn'='l) 

n 

PN- 

where 


V  =  WVW’VWW 

p  -  rl  ^n^W'^n'^W 

N  '  slHn2i(:'o,JA(>Il)-soHn2l',VJn<’t.) 

M’WVW^oWV*!1' 


Again,  for  normal  incidence. 


,TM 

'n 


V 

TT 


,TE  _  M 
'n  ~  "  -F  * 


(26) 


(27) 


(28) 

(29) 


2.4  Perfectly  Conducting  Cylinder 

Setting  €  =  03 ,  Eqs.  (19)  and  (23)  for  the  dielectrically  clad  cylinder  reduce 
to  the  following  expressions  for  the  perfectly  conducting  cylinder  at  oblique 
incidence: 


C™  - 


Jn(5£o) 


w 

ifTTT 


(30) 


(31) 


Equations  (30)  and  (31)  are  of  the  same  form  as  the  well-known  equations  for 

the  perfectly  conducting  cylinder  at  normal  incidence,  except  that  at  oblique 

incidence  x  =  k  R  cos  $  .  In  all  of  the  above,  the  factors  E^  cos  ❖  and  H  cos 

oooo  oooo 

have  been  surpressed  for  all  of  the  Cn's.  It  is  interesting  to  note  that,  ignoring 
the  factors  Eocos^0  or  H0cos^,  the  a  for  the  perfectly  conducting  cylinder  at 
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oblique  incidence  is  identical  to  that  of  normal  incidence  w ith  the  radius  reduced 
by  the  factor  cos  *  . 


3.  COMPLIED  NORMALIZED  SIGMA 

The  normalized  radar  cross  section 


S  = 


a 

zR  cos  $ 
o  o 


(32) 


was  computed  for  perfectly  conducting  solid,  homogeneous,  dielectric,  and  dielec¬ 
trically  clad  cylinders, 

3.1  Metallic  Cylinder 

For  the  metallic  cylinder,  a  single  computation  of  a  versus  Rq  for  both  TE  and 
TM  incident  modes  is  sufficient  since,  for  other  than  normal  incident,  multiplying 
Rq/A  by  cos  'I'Q  will  convert  normal  incidence  to  another  angle  of  incidence  'i'o- 

Figure  2  is  the  well-known  o7rRQ  of  normal  incidence  for  a  perfectly  con¬ 
ducting  infinite  cylinder  with  the  TM  mode  decreasing  monotonically  with  increas¬ 
ing  Rq/A  and  the  TE  mode  increasing  monotonically  to  unity  and  then  oscillating 
with  damped  oscillations  as  Rq/ A  is  increased  and  finally  settling  to  the  limiting 

value  of  air, R^  =  1. 

o 

3*2  Solid  Dielectric  Cylinder 

In  Figure  3,  the  normalized  sigmas  at  normal  incidence  for  TE  and  TM 
incident  modes  are  given  for  0.0002£  R  />  £1.0  and  e  =  4.  0.  Here  the  normalized 
sigmas  for  both  TE  and  TM  modes  are  monotonically  increasing  for  values  of 
R  / A  <  0. 10,  and  then  both  commence  to  oscillate  with  the  oscillations  increasing 
in  amplitude  and  rising  above  unity  with  increased  Rq/A.  It  is  interesting  that  the 
TM  mode  for  small  R0/A  is  approximately  at  the  same  level  as  the  TE  mode  for 
the  perfectly  conducting  cylinder.  The  level  of  course  depends  on  the  value  of  c 
and  will  be  higher  for  larger  values  of  e  and  lower  for  the  lower  values  of  e.  In 
Figure  4,  the  Rq/A  versus  S  is  given  for  e  =  4  and  =  84°.  In  comparison  with 
the  metallic  cylinder,  the  TE  mode  is  considerably  higher  up  to  Rq/A  =  0.  2  with 
the  TE  and  TM  modes  very  close  in  amplitude. 

It  has  been  found  that  the  relative  number  of  fluctuations  (min  to  min)  is  related 
to  the  dielectric  constants  of  the  materials  as  "Vcj  -  1  /  ,  where  Cj  is  the 

permittivity  oi  one  cylinder  and  Cg  the  permittivity  of  another  cylinder. 
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Figure  2.  Perfectly  Conducting  Infinite  Cylindei .  Normalized  Sigma  vs  R  /\,  y  =  0%  R  /X  =  0.  0002 
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Figure  3.  Solid  Dielectric  Infinite  Cylinder.  Normalized  Sigma  vs  I50/\. 
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Figure  4.  Solid  Dielectric  Cylinder  vs  Metallic  Cylinder.  Normalized 
Sigrnas  vs  Rq/A.  =  84°,  e  =  4.  0,  R^/A.  =  .  0002  —  1.  0. 

3.3  Dielectrically  Clad  Cylinder 

The  normalized  sigma  versus  R  l'\  for  the  range  0.  0002  £  R Q/\  ^  1  for 
shell  thicknesses  as  determined  by  R^  =  0.  9R0  and  dielectric  constant  of  the 
shell  equal  to  4.  0  at  normal  incidence  is  given  in  Figure  5.  Here  again, 
unlike  the  perfectly  conducting  cylinder,  the  oscillations  of  the  TE  mode  do  not 
dampen  out  in  the  same  manner  as  the  TE  mode  of  the  perfectly  conducting 
cylinder;  rather,  they  increase  in  amplitude  after  R Q/X  =  0.  5.  For  small  R o/.\, 
the  TE  mode  of  the  dielectrically  clad  cylinder  is  significantly  higher  than  that  of 
the  perfectly  conducting  cylinder,  and  the  TM  mode  is  significantly  lower  than  that 
of  the  perfectly  conducting  cylinder. 
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Figure  5,  Dielectrically  Clad  Infinite  Cylinder.  Normalized  Sigma  vs  R  /\. 
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It  is  interesting  to  compare  the  dielectrically  clad  cylinder  at  normal  incidence 
with  the  dielectrically  clad  sphere*  ,  Figure  6.  Note  the  considerable  similarity 
between  the  sphere  and  the  cylinder  at  normal  incidence.  Both  the  sphere  and 
cylinder  have  burst  of  oscillations  that  dampen  out  and  then  resume.  The  first 
burst  of  oscillations  from  the  sphere  corresponds  with  the  TE  mode  of  the  cylinder, 
and  the  second  burst  with  the  TM  mode  of  the  cylinder;  when  the  oscillations  of 
both  modes  of  the  cylinder  are  damped,  so  also  are  those  of  tue  sphere.  Both  sphere 
and  cylinder  have  equal  numbers  of  oscillations  occuring  at  the  same  Rq/X;  however, 
the  sphere  differs  from  the  cylinder  in  that  the  amplitude  of  the  oscillations  is  larger 
for  the  sphere  than  fox'  the  cylinder. 

Figure  7  shows  normalized  sigma,  S,  versus  Rq/X  for  a  shell  thickness  of 
0. 1  Rq,  e  =  4.  0,  and  the  angle  of  incidence  equal  to  84°,  for  0.  0002  ^  Rq/X:£  1.  0 
in  comparison  with  the  perfectly  conducting  cylinder.  Here  the  S  for  the  TE  mode 
of  the  dielectrically  clad  cylinder  is  significantly  larger  than  the  corresponding 
mode  of  the  perfectly  conducting  cylinder  for  Rq/X^  0.  2. 

Figures  8  and  9  give  scan  in  angle  for  Rj  =  0.  9Rq  and  Rq  =  0.  78  X,  respec¬ 
tively,  for  e  =  4.  0  and  e  =  2.  56.  Here  there  is  considerable  similarity  between 
the  two  values  of  e;  however,  for  the  large  e,  the  S  rises  higher  above  unity 
than  it  does  for  the  smaller  e. 

The  normalized  or,  S,  versus  e  is  given.  i,t  Figure  10  for  if  =  45°,  R Q/X  =  0.  78 
and  Rx  =  .  702X. 


*  The  sphere  was  computed  by  J.  Rheinstein  of  MIT,  Lincoln  Laboratories 
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Figure  7.  Dielectrically  Clad  Infinite  Cylinder.  Normalized  Sigma  vs  RQ/\,  *0=84°,  e  =  4.0,  0.9Ro,  RQ/X  =  0.  0002' 
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Figure  8,  Dielectrically  Clad  Infinite  Cylinder.  Normalized  Sigma  vs  Angle  of  Incidence.  e=4.0,  R.=0.  9R  ,  R  ”  78X 
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XL  Effect  of  Surface  Diffusivity  Upon  the  Scattering 

Characteristics  of  a  Plasma  Sphere* 

Philip  J.  Hyati*  * 
Plasmadyne  Corporation 


Abstract 


The  effect  of  varying  the  suriace  diffuseness  of  a  spherical  plasma  region  is 
examined  in  the  resonance  region.  It  is  found  that,  large  changes  in  electron  den¬ 
sity  may  be  readily  compensated  by  modest  changes  in  surface  diffuseness  as 
regards  the  scattering  characteristics  of  such  regions. 


Since  the  re-entry  characteristics  of  a  great  many  objects  often  include  large 
regions  of  highly  ionized  material,  it  is  of  particular  interest  to  consider  how 
slight  structural  modifications  of  these  plasma  volumes  can  appreciably  affect  the 
scattering  of  electromagnetic  waves  in  the  resonance  region.  For  purposes  of 
illustration,  a  hypothetical  plasma  will  be  examined  whose  structure  is  assumed 
to  be  spherically  symmetric.1  In  this  event,  an  earlier  formalism  developed  by 


This  work  was  completed  while  the  author  was  at  Plasmadyne  Corporation  and 
was  performed  entirely  at  company  expense.  Much  of  this  material  may  be 
found  in  J.  Appl.  Phys.  34:2078/1963. 

Now  a  member  of  the  Technical  Staff  of  the  Defense  Research  Corporation, 
Santa  Barbara,  California. 
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*)  O  4 

the  author  and  others  '  will  become  immediately  applicable  as  soon  as  an  asso¬ 
ciation  between  the  conventional  plasma  variables  and  the  so-called  index  of 
refraction  is  made. 

For  a  plasma  in  local  thermodynamic  r  j  dlibrium,  the  electron  density,  N, 

and  collision  frequency,  v,  are  sufficient  parameters  io  characterize  its  gross 

structural  features.  The  minor  contributions  from  the  heavier  ions  are  usually 

neglected  in  this  approximation.  The  electron  collision  frequency  itself  is  in 

5 

general  a  rather  complicated  function  of  temperature  and  electron  density. 
Defining  the  plasma  frequency, &p  ,  by 

o2  =  4 it  Ne2/m  ,  (1) 

P 

one  may  introduce  the  local  dielectric  constant 


whicn  depends  upon  the  incident  angular  frequency  as  indicated.  The  electron 

charge  and  mass  are  indicated  by  e  and  m  respectively. 

If  now  the  plasma  frequency  and  collision  frequency  are  functions  of  position, 

then  e  also  will  be  a  function  of  position.  The  complex  refractive  index  is  given 

by  the  square  root  of  e  and  thus  also  varies  with  position.  Since  all  variations 

are  for  the  present  example  assumed  radial,  one  has  the  situation  described  by 
2 

Wyatt  ,  namely  a  medium  whose  refractive  index  is  a  function  of  r  only. 

Consider  now  the  following  situation:  Radiation  of  frequency  w/2tt  is  incident 

upon  a  spherically  symmetric  plasma  with  a  uniform  interior.  In  terms  of  the 

vacuum  wavelength,  A.  the  dimensionless  parameter  p  =  2jrrA  may  be 

VoC  vac 

introduced.  Thus 

n(p)  =  ^efp)  =  Vl  -  «  (p)2  [j.  -  iv{p)/ofj  /  Q*j2  +  v(p)2J  •  (3) 


For  p  less  than  some  value  a,  say,  both  w  and  v  are  assumed  constant  and 
therefore  one  may  write 

n(p)  =  nQR  +  inQI,  o<  p  <  a  .  (4) 

Surface  diffusivity  may  now  be  introduced  by  letting  n(p)  vary  from  the  value 
given  by  Eq.  (4)  to  the  value  1.0  at  some  point  b(>  a),  say,  by  letting 
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n{p)  =  nR(p)  <■  in j(p),  a<  p  <b  (5) 

nip)  =  1.0,  p  >  b 

where  nR(p)  and  nj{p)  are  real  functions  of  P  .  This  functional  dependence  of  the 
complex  refractive  index  may  be  deduced  by  assuming  some  realistic  variation 
with  P  of  a.p  and  v  or,  alternatively,  by  assuming  a  variation  of  n(p)  directly. 
Although  either  approach  will  yield  a  diffuse  surface  with  similar  characteristics, 
the  latter  assumption  is  more  convenient  for  the  present  discussion. 

In  order  that  n(p)  and  its  first  two  derivatives  be  continuous  throughout  the 
entire  range  of  interest,  a  convenient  ioi  m  factor  of  the  Green-Wyatt  type  may  be 
introduced  as  follows: 

n(p)  =  nR(p)  +  inj(p) 

where 

nR(p)  =  1  +  (n0R-l)  |  (p)  (6) 

nl(p)  =  nQI  |(p) 

and  the  form  factor  ?(p)  is  given  by 
s(p)  =  1,  o  <  p  <a 

=  -  -  15Z/16  +  10Z3/16  -  3Z5/lt),  a<  p  <  b  (7) 

=  0,  p>  b 

Z  =  Q>  -i(b  +  a)  /  ^-{b-a)J  . 


The  parameters  a  and  b  may  be  related  to  the  half  fall-off  distance 
p  (=2-irr  A,,_  =  kr  )  and  the  surface  thickness  T  (diffuseness),  if  one  defines 

v  U  V  dv_.  U 

the  term  “thickness"  to  mean  the  0.9  to  0.1  fall-off  distance.  Thus  for  the  form 
chosen  in  Eq.  (7) 

a  =  PQ  ~  0.99868T  (8) 

b  =  pQ  +  0.99868T. 

The  variables  of  interest  are  illustrated  in  Figure  1  for  the  case  where 
nOR  =  1,37  an<*  nQI  = 
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Figure  1.  The  complex  index  of  refraction  for  an  object 
with  a  diffuse  surface.  The  form  factor  is  of  the  Green- 
Wyatt  type  and  the  radial  variation  is  in  units  of  1/k  ; 

i  p  ~  kvac  7 

For  the  present  example  the  incident  frequency  has  been  fixed  at  3kMc, 

v  =  10*°sec-1,  and  pQ  =  5.  Thus  the  problem  concerns  an  S-band  measurement 

of  a  plasma  whose  mean  radius  is  about  8  cm. 

Consider  first  the  effect  upon  the  radar  cross  section  of  such  a  region  when 

the  surface  thickness  is  varied.  The  curves  shown  in  Figure  2  illustrate  this 

1/2 

result  for  a  variety  of  values  of  '  )  at  the  uniform  interior  region  of  the 

plasma  scatterer.  Rather  than  express  the  results  in  terms  of  the  radar  cross 

-2 

section,  the  abscissa  refers  to  the  backscattered  intensity  in  units  of  (kr) 
normalized  to  unit  incident  flux. 


Figure  2.  The  differential  scattered  intensity  of  unpolarized  radiation  in  units  of 
(kr)"2  at  180°  shown  for  various  surface  thicknesses  as  a  function  of  plasma  fre¬ 
quency  (pQ=5,  a?  =  6tt  X  H)S  3ec_i,  v  =  lO^O  sec-*) 
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Several  important  results  are  clearlv  illustrated  by  Figure  2.  At  the  smaller 
electron  densities,  that  is,  when  the  interior  plasma  frequency  is  less  than  the 
incident  angular  frequency,  small  variations  in  surface  thickness  cause  appreciable 
variations  in  the  backscattered  intensity.  This  would  be  expected  in  general  since 
an  increased  surface  diffuseness  usually  results  in  a  diminished  reflection.  This 
trend  is  continued  throughout  the  region  where  the  plasma  frequency  is  comparable 
to  the  incident  angular  frequency,  except  for  the  particularly  sharp  surfaces.  For 
these  latter  exceptions,  the  xeilection  coefficients  are  critically  affected  for  even 
slight  variations  in  the  surface  region  a,,d  thus  produce  the  resonance  type  phenomena 
illustrated.  It  is  interestingto  note  that  from  a  given  measurement  of  backscattered 
intensity,  it  is  in  general  impossible  to  deduce  anything  concerning  the  structure 
of  the  scatterer.  This  result  is  vividly  illustrated  in  Figure  2  since  a  given  back- 
scatter  may  be  produced  from  an  infinite  variety  of  sp  and  T  combinations. 

Not  only  is  the  radar  cross  section  appreciably  affected  by  variations  in  sur¬ 
face  diffuseness,  but  the  scattering  at  all  other  angles  as  well.  Assuming  for  the 

$ 

moment  that  the  source  of  radiation  is  unpolarized,  then  Figures  3  and  4  illus¬ 
trate  the  differential  scattered  intensities  for  various  interior  plasma  frequencies 
at  two  distinct  surface  thicknesses  (0.0  and  1.0,  respectively).  Thus  for  the  case 
of  larger  surface  thickness,  the  scattering  characteristics  at  a  particular  angle 
behave  in  an  essentially  monotonic  manner.  For  the  sharp  surface  (Figure  3)  on 
the  other  hand,  this  simple  behavior  is  no  observed.  Indeed,  the  sharp-surfaced 
structure  may  be  said  to  have  a  high  Q  (at  a  variety  of  dielectric  constants)  which 
is  easily  "spoiled",  thus  resulting  in  the  resonance  behavior  illustrated.  The 
diffuse  structure  may  be  said  to  always  have  a  poor  Q. 

If  the  interior  plasma,  frequency  is  fixed,  then  the  diffusivity  effects  upon  the 
differential  scattered  intensity  may  be  examined  directly.  Figure  5  illustrates  this 
situation  for  w  p  ~  lO^sec  *.  Note  the  general  shifting  of  the  diffraction  pattern 
toward  smaller  angles  with  increasing  surface  diffuseness.  For  less  than 
this  value  a  significant  drop  in  the  backscattered  intensity  is  also  observed. 

The  above  example  illustrates  that  marked  variations  in  the  scattering  charac¬ 
teristics  of  plasma  volumes  may  in  general  be  expected  subsequent  to  small  varia¬ 
tions  in  surface  diffusivity.  In  addition,  large  changes  in  ionization  density  may  be 
easily  "disguised,"  insofar  as  conventional  radar  measurements  are  concerned,  by 
relatively  small  changes  in  surface  diffuseness.  Applied  to  larger  systems  at  other 
frequencies,  the  above  indicated  results  could  have  interesting  practical  applications. 

Although  this  is  completely  unrealistic  insofar  as  microwave  sources  are  con¬ 
cerned,  this  assumption  permits  an  average  differential  scattering  pattern  to  be 
presented.  As  was  so  vividly  brought  out  at  the  Symposium,  i*  is  not  enough  to  be 
just  concerned  with  the  effects  upon  the  scattering  of  certain  modifications  or  load¬ 
ing  techniques  for  a  given  polarization  of  the  incident  radiation.  The  effects  upon 
two  mutually  perpendicular  polarizations  must  always  be  considered  simultane¬ 
ously.  Thus,  the  assumption  of  an  unpolarized  source  (which,  incidentally,  renders 
the  computations  somewhat  more  difficult)  results  in  a  clear  picture  of  the  net  effects. 
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1.  INTRODUCTION 

It  is  the  purpose  of  this  paper  to  report  on  theoretical  work  dealing  with  the 
absorption  of  electromagnetic  power  by  plasma  spheres.  The  electromagnetic 
energy  is  propagated  to  the  plasma  sphere  in  the  form  of  plane-wave  radiation. 

The  sphere  is  assumed  to  be  homogeneous  and  to  be  surrounded  by  a  medium 
whose  refractive  index  is  near  unity.  Under  these  conditions,  Mie  scattering 
theory  may  be  used  to  compute  the  absorption  cross  section  of  the  sphere  as  a 
function  of  its  complex  dielectric  constant,  its  radius,  and  the  frequency  of  the 
incident  electromagnetic  wave.  In  particular,  the  conditions  for  maximum  power 
absorption  efficiency  will  be  discussed.  It  will  be  shown  that  strong  absorption 
resonances  are  to  be  expected  only  for  cases  where  the  ratio  of  sphere  radius  to 
free  space  electromagnetic  wavelength  is  less  than  unity.  Consequently,  models 
which  assume  infinite  or  semi-infinite  media  are  not  applicable  and  would  not  yield 
the  absorption  resonances  that  are  predicted  by  the  application  of  electromagn 
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scattering  theory.  On  the  other  hand,  treatments  based  on  electrostatic  or  mag¬ 
netostatic  solutions  o'  >.he  wa»»e  equation  /.  ill  yield  a  discrete  set  of  characteristic 
" resonance"  conditions  in  the  dielectric  constant  or  m  the  permeability  but  will 
not  lead  to  a  determination  of  absorption  cross  section  or  to  the  conditions  for 
absorption  cross  section  maxima. 


2.  APPROACH 

Since  it  is  the  objective  of  this  work  to  optimize  the  efficiency  for  power 
absorption  in  dielectric  spheres,  the  results  of  formal  Mie  theory  are  used  to 
write  down  the  exact  expression  for  absorption  cross  section.  It  is  first  deter- 
minted  that  strong  absorption  resonances  are  to  be  expected  only  in  the  low  fre¬ 
quency  limit.  A  low  frequency  approximation  is  then  applied  which  reduces  the 
problem  to  that  of  optimizing  the  absorption  cross  section  of  lossy  spheres  having 
small  radius  to  wavelength  ratio.  This  is  done  by  first  determining  the  optimum 
values  of  the  real  and  imaginary  parts  of  the  appropriate  partial  wave  coefficient 
and  then  computing  the  maximum  possible  value  for  the  power  absorption  efficiency 
as  a  function  of  ka,  where  k=2  sr/X  and  a  =  sphere  radius.  Next,  the  optimum 
values  of  the  real  and  imaginary  parts  of  the  dielectric  constant  are  found  as 
functions  of  ka  by  means  of  performing  a  power  series  expansion  of  the  partial 
wave  coefficient  in  terms  of  ka.  In  this  step,  the  condition  that  ka  <<  1  is  again 

used.  Since  the  magnitude  of  the  optimized  absorption  efficiency  is  proportional 
_2 

i  y  (ka)  ,  and  since  the  approximations  used  become  increasingly  accurate  as  ka 
decreases,  this  method  should  be  quite  valid  in  regions  of  very  high  absorption 
resonance. 

Finally,  the  effect  of  an  applied  magnetic  field  is  investigated  and  it  is  found 
that  the  existence  of  the  magnetic  field  serves  to  decouple  the  parametric  equations 
for  absorption  resonance. 


3.  RESULTS  OF  FORMAL  MIE  THEORY 

The  absorption  efficiency  factor  Qabs  for  a  scatterer  is  defined  as  the  ratio 

of  its  absorption  cross  section  to  its  projected  area.  For  a  sphere  this  is 
2 

^abr  =  °abs^ra  *  w^ere  ^abs  =  al5SorPtion  cross  section  and  a  =  sphere  radius. 

Likewise,  the  extinction  and  scattering  efficiency  factors  are  defined  as 
2  2 

Q  =  <r  ./jra  and  Q  =  <r  /it a  .  From  conservation  of  energy,  it  follows 

ext  ext  sea  sea  ' 

that 


i  ' 


^abs  "  ^ext  ”  ^sca  '  ^ 

The  extinction  and  scattering  efficiency  factors  are  usually  given  by. the  following 
series  representations: 


where  x  =  ka  =  2«a/\,  and  an  and  bfi  are,  respectively,  the  coefficients  for  the 
n'th  electric  and  magnetic  partial  waves.  These  can  be  shown  to  be  the  contribu¬ 
tions  due  to  the  induced  electric  and  magnetic  dipole,  quadrupole,  and  higher-order 
moments  of  the  sphere.  They  may  be  written  as  follows: 


\  nn(x)  (yjn(y)]  '-  <rjn(y)  [xnn(x)T  j 
I  jn(x)  ^jn(y)l  ejn(y)  txjn(x))  '  \ 


b  = - - - - -  ,  (5) 

(  «rnn(x)  [yjn(y)J '  -  jp(y>  [xnn(x)] '  j 

j  cjn(x)  (yjn(y)] '  -  jn(y)  [x  jn(x)] '  j 

where  the  Jn(y)  and  nn(x)  are  spherical  Bessel  functions  of  the  first  and  second 
kinds  respectively,  x  =  ka,  y  =  x  -\JT ,  and  e  =  complex  dielectric  constant.  The 
prime  in  Eqs.  (4)  and  (5)  indicate  derivatives  with  respect  to  the  argument.  Since 
the  Bessel  functions  will  in  general  be  complex,  the  an  and  bn  will  also  in  general 
be  complex.  Substituting  Eqs.  (2)  and  (J)  into  Eq.  (1),  we  have 

«  2  2 
E  !«•«*„> -KJ  +Re<bn>- |b/J.  w 

ra  x  n»l  1  1 

Equation  (6)  indicates  that  it  is  reasonable  to  expect  resonances  in  Qa^s  when 
resonance  exists  for  one  or  more  of  the  partial  wave  coefficients.  Such  resonances 
can  only  occur  for  plasma-like  media  with  e'<  1,  as  shown  by  approximations  of 
the  first  few  partial  wave  coefficients. 
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A  computer  code  for  the  computation  of  the  Q's  has  been  developed  by  J.  C. 
Pedersen  of  the  Avco  Corporation.  Some  results  of  this  program  have  been 
checked  against  published*  results  of  a  similar  program  at  RAND  Corporation 
and  good  agreement  was  obtained  for  ka  values  as  low  as  0.  2  (the  lower  limit  on 
the  published  RAND  results).  Both  programs  were  operable  for  dielectric  con¬ 
stants  having  a  negative  real  part  as  well  as  an  imaginary  part.  Such  cases 
represent  plasma  spheres.  In  the  case  of  a  plasma,  the  complex  dielectric 
constant  is  given  by 


ie"  =  1  - 


i"ee  1 

17~M - 2 

o  e  u) 

<x) 

-JUi 


(7) 


where  ng  =  electron  density,  e  =  electronic  charge,  «  =  electron  collision  fre- 
quency,  and  e  ~  permitivity  of  free  space.  It  is  convenient  to  define  v  =  (u>  /o>) 

2  °  P 

ne  e 

= - -  and  y  =  (a?  /«),  where  u>  is  the  plasma  frequency.  The  real  and 

<oMe* 

imaginary  parts  of  the  dielectric  constant  are  then  given  by 


_  1  +  y"  - 
c - 2 - 


(8a) 


1  +y 


=  ^  y 

1  +  y 


2  ' 


(8b) 


Some  general  statements  can  be  made  about  the  limiting  values  of  rj  and  y  as 
related  to  absorption.  First,  it  is  obvious  that  no  absorption  can  occur  for  e"  =  0. 
Thus,  the  cases  of  tj  =0  (no  electrons)  or  y  =  0  (no  collisions)  will  yield  Qa^s  -  0. 
The  case  of  t)  =  °°  with  y  finite  represents  a  perfect  conductor  (imaginary  refractive 
index).  No  absorption  can  occur  in  this  case.  Finally,  for  y  =°°  and  rj  finite,  we 
have  e"  =  0  and  again  no  absorption  can  occur.  It  has  therefore  been  shown  that 
any  resonances  in  the  absorption  cross  section  will  occur  for  finite  values  of  both 
t j  and  y. 

A  number  of  sets  of  computer  runs  were  made  for  plasma-like  spheres,  and 
the  results  are  shown  in  Figures  1  through  5.  Each  set  of  curves  is  a  plot  of  Qabs 
vs  y  for  a  particular  fixed  value  of  ka,  with  r\  the  running  parameter.  The  max¬ 
imum  value  of  Qaks  for  a  given  ka  is  seen  to  be  characterized  by  a  discrete  value 

*D.  Deirmendjian  and  R.  J.  (Jlasen,  Light  Scattering  on  Partially  Absorbing 
Homogeneous  Spheres  of  Finite  Size,  HAND  Report  R- 3^3 -PR,  Feb.  1962T 
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Figure  3.  Qabg  as  a  Function  of  i?  and  y  for  ka  =  0.  8 


r— *■ 


abs 


Figure  5.  Q 


as  a  Function  of  r?  and  y  for  ka  -  2.  5 
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of  »}  and  a  discrete  value  of  y  .  Thus,  for  a  given  ka,  this  set  of  values  of  r?  and 
y  represents  the  resonance  condition  for  Qabs.  Of  particular  interest  in  these 
results  is  the  fact  that  the  maximum  value  of  Q  -  .  which  we  define  as  Q  .  (max), 

aOS  305 

increases  strongly  with  decreasing  ka.  The  various  Qa^s(max)  in  Figures  1  to  5 
have  been  identified  with  resonances  in  a^ ,  the  electric  dipole  partial  wave  coeffi¬ 
cient.  For  ka  values  below  0.  4,  only  the  a^  coefficient  need  be  considered. 
Although  resonances  due  to  the  higher  order  coefficients  have  been  observed  in 
some  computer  runs,  these  are  very  much  sharper  than  those  due  to  the  electric 
dipole  term  and  arise  from  values  of  y  which  are  too  small  to  be  of  practical 
interest.  Since  the  magnetic  permeability  was  set  equal  to  uq,  no  magnetic 
resonances  were  expected  in  these  runs,  nor  were  they  found.  Resonances  due  to 
*he  bR  can  cccur  only  in  the  case  of  negative  permeability  and  were  therefore  not 
investigated. 


4.  APPROXIMATIONS  TO  MIE  THEORY 


Since  small  ka  implies  large  Qa^s(max),  and  since  it  has  been  demonstrated 
that  the  use  of  alone  v.ill  suffice  under  these  conditions,  it  is  next  meaningful 
to  make  two  approximations. 

(a)  Assume  that  Qabs  results  from  only  the  electric  dipole  term  a^. 

Find  Qabs<max)  as  a  function  of  ka.  In  doing  this,  it  will  be  necessary 
to  find  the  optimum  values  of  the  real  and  imaginary  parts  of  a5  for  the 
occurrence  of  Qabg(max). 

(b)  Write  a^  as  a  power  series  expansion  in  ka  and  make  the  appropriate 
approximation  for  ka  <<  1.  Separate  the  approximate  expression  of  a, 
into  its  real  and  imaginary  parts.  The  optimum  value  of  these,  found 
from  (a)  above,  will  then  be  used  to  arrive  at  two  equations  which  define 
the  necessary  values  of  c*  and  e"  as  a  function  of  ka  for  the  occurrence 
of  Qabs  (max).  Finally,  Eq.  (8)  will  be  used  to  obtain  a  set  of  equations 
in  t)  and  y  giving  the  resonance  values  of  these  as  functions  of  ka. 

For  the  determination  of  Qabs^maxK  assuming  a  particular  afi  to  be  the  only 
contributor,  Eq.  (6)  reduces  to 


<5abs  =  4!2'I+1)[Re(an)  -  k|2+1VV  -  I6,/} 


(9) 


If  we  now  let  a  =  (c  +  id) 
n 


we  obtain 
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£W'Jr<2n+1> 


To  find  Qabs  =  Qabs{rrjax),  we  must  have 


^  ,  a  _  ( 2n+l) 

Qahs(max) - 7— 


2x 


(10) 


(li) 


(12) 


or 


Q 


abs 


(max)  = 


(electric  dipole  resonance). 


(13) 


Exact  machine  calculations  verify  Eq.  (13)  for  values  of  x  below  0.  5. 

The  next  step  is  to  calculate  the  optimum  properties  of  the  medium  as  functions 

of  x.  It  will  be  most  convenient  to  first  find  the  absorption  resonance  conditions: 

e'  =  fj(x)  and  c"  =  gj(x).  Equation  (8)  will  then  be  used  to  find  ij  =  fgfx)  and 

Y  =  g0(x).  For  this  purpose,  it  is  convenient  to  write  as  a  pov/er  series 

expansion  in  x  and  to  then  make  a  low  frequency  approximation,  making  sure 

that  sufficient  terms  are  kept  in  the  denominators  to  prevent  the  occurrence  of 

infinities  which  would  otherwise  result.  Due  to  the  algebraic  complexity  of  this 

operation,  the  series  expansion  for  the  general  a  term  is  given  in  the  Appendix. 

n  -1 

The  result  is  given  by  Eq.  (A -12).  Since  Eq.  (11)  led  to  the  result  a  =  (2  +  iO) 
for  Qabs  =  Qabs(max),  we  obtain  the  absorption  resonance  conditions  by  equating 
the  real  and  imaginary  parts  of  the  denominator  of  Eq.  (A- 1 2)  to  2  and  0  respec¬ 
tively,  which  yields  Eqs.  (A-13)  and  (A-14).  We  now  let  x  <<1.  The  zero'th 

order  approximation  of  the  absorption  resonance  conditions  in  Eqs.  (13)  and  (14) 
n  |  ^ 

is  a  = - - — ,  /)  =  0.  This  is  obviously  not  of  sufficiently  high  order  in  x  since 

the  x  dependence  on  e"  and  e"  does  not  appear.  In  addition,  this  result  gives 
zero  collision  frequency.  The  next  higher -order  approximation  gives 


(n+1)  (2n+ 1) 


[(2  n+ 1) !  l]  1  (n+ 1)  n 


7T  x 


2n  +  1 


a  =  - 


n+1 


1  +  x 


2(2  n+  1) 
n(2n-  1)  (2  n+  3) 


(14) 


(15) 


in  the  case  of  electric  dipole  resonance  (n-=-  1),  we  have  the  resonance  conditions: 
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jS  *  2x3 


(16) 


a  =  - 


(17) 


Substituting  Eqs.  (16)  and  (17)  into  Eqs.  (8a)  and  (8b),  we  finally  get 


=  ^  =  2 
a  1 


3 


x  . 


(18) 


cv 


3  +  i|x2. 


(19) 


Equations  (16)  and  (17)  constitute  the  necessary  conditions  for  electric  dipole 
resonance  in  the  absorption  efficiency  factor  for  homogeneous  plasma  spheres 
having  small  radius-to-wavelength  ratios.  Since  Eq.  (13)  shows  that  x  must  be 
small  in  order  for  large  absorption  efficiencies  to  occur,  Eqs.  (18)  and  (19)  are 
useful  in  the  region  of  interest. 


5.  DISCUSSION 

It  was  demonstrated  in  the  last  section  that  large  absorption  resonances  can 
occur  in  plasma  spheres  (Eq.  13).  The  conditions  for  electric  dipole  absorption 
resonance  (Eqs.  18  and  19)  show  that  in  the  region  of  strong  resonance  (small  x), 

3 

the  value  of  tj  goes  to  a  constant  value  of  3,  while  y  is  proportional  to  x  .  Thus, 

_3 

the  ratio  of  plasma  frequency  u>p  to  collision  frequency  is  proportional  to  x  . 
Thus,  for  Qabs  =  150  (Eq.  13),  we  would  require  Up/a^  =  VrT / y  ss  2600.  As  x  is 
further  decreased,  this  ratio  will  greatly  increase.  It  is  probably  this  fact  more 
than  any  other  which  will  limit  the  highest  physically  realizable  values  of  Qabg. 
Another  consideration  is  the  mean  free  path  of  electrons  within  the  plasma  sphere. 
This  distance  is  really  the  mean  distance  a  "hot"  electron  must  travel  in  order  for 
it  to  equilibrate  with  the  surrounding  ion-neutral  atom  medium.  Obviously,  the 
sphere  radius  should  not  be  smaller  than  the  mean  free  path.  On  the  other  hand, 
the  radial  homogeniety  of  the  plasma  sphere  is  not  believed  to  be  too  important  in 
this  case.  The  reason  is  that,  since  a/X  is  small,  the  wave  will  "see"  a  sphere 
with  some  effective  dielectric  constant  and  some  effective  radius.  These  assump¬ 
tions  have  been  at  least  partially  confirmed  by  Dr.  ?.  Wyatt*  ,  Defense  Research 
Corporation. 


*  Private  communication. 
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The  only  real  way  to  predict  the  degree  of  physical  realizability  of  this  absorp¬ 
tion  x*esonance  phenomenon  is  to  combine  the  electro-magnetic  resonance  conditions 
(Eqs.  18  and  IS)  with  appropriate  physical  expressions  for  the  equation  of  state  and 
the  collision  frequency,  and  to  derive  a  set  of  equivalent  thermodynamic  resonance 
conditions.  A  preliminary  calculation  of  this  nature  was  performed  by  one  of  the 
writers  (N.  P. )  who  used  the  weakly  ionized  approximation  of  the  Saha  equation. 
Difficulties  were  encountered  in  that,  in  the  region  of  interest  (small  x),  the 
assumption  of  weak  ionization  was  not  valid.  Mr.  J.  Davis  of  the  Avco  Corporation 
is  presently  extending  this  work  to  include  strong  ionization.  These  results,  when 
they  are  available,  will  be  investigated  for  violations  of  the  assumption  of  equilib¬ 
rium,  etc.  If  the  assumption's  implicit  in  the  calculation  prove  to  be  valid  in  the 
resonance  region,  we  will  conclude  that  this  is  the  proper  thermodynamic  model. 

If  not,  we  will  proceed  to  the  more  complicated  non-equilibrium  case.  Present 
indications  are  that  this  may  be  necessary. 

One  of  the  v/riters  (N.  P. )  has  done  some  preliminary  work  on  the  effects  of  a 
strong  magnetic  field  on  the  absorption  resonance  conditions.  The  anticipated 
results  are: 

2 

W  12  2 

*3+-±£x  ,  (20) 


OJ 


2 

I 


(2i) 


where 

g/3  « 

o)  „  = - —  =  electron  cyclotron  frequency. 

p  e 

Although  Eqs.  (21)  and  (18)  are  identical,  Eqs.  (20)  and  (19)  differ  in  that 
replaces  in  Eq.  (19).  Thus,  the  application  of  a  strong  (w^>>w)  magnetic 
field  serves  to  decouple  the  resonance  conditions,  and  leads  to  an  extra  degree  of 
freedom  in  physically  realizing  strong  absorption  resonances.  If  severe  limita¬ 
tions  should  be  shown  as  a  result  of  the  zero  magnetic  field  calculations,  the 
strong  field  case  will  be  investigated  in  detail  regarding  its  relaxation  of  thermo¬ 
dynamic  requirements. 

When  the  correct  thermodynamic  description  of  the  plasma  has  been  found, 
the  power  input  and  the  power  output  will  both  be  written  in  analytic  form.  Thus, 
an  analytic  equation  for  steady  state  can  be  written.  This  equation  will  be  set  up 
ard  the  required  incident  electric  field  strength  will  be  computed.  A  more  interest 
ing  study  will  be  an  examination  of  the  steady  state  equation  for  stability.  Is  the 
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plasma  sphere,  under  the  (initial)  steady  state  conditions  stable  against  fluctuations 
of  the  thermodynamic  parameters?  This  question  will  hopefully  be  answered  when 
the  thermodynamic  resonance  conditions  are  expressed  in  analytic  form. 
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